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ABSTRACT 


The  objective  of  this  thesis  is  to  investigate  the  use  of  the  Auxiliary  Resonant 
Commutated  Pole  Converter  (ARCP)  as  a  DC-to-AC  power  converter.  The  advantages 
and  disadvantages  of  resonant  converters  over  non-resonant,  or  hard-switched  converters, 
are  investigated.  Basic  ARCP  circuit  operation  is  modeled,  with  emphasis  placed  on 
examining  the  commutation  between  high  and  low  voltage  states.  Detailed  ARCP 
converter  operation  is  modeled  in  software  and  compared  to  a  software  model  of  a  hard- 
switched  converter.  Comparisons  are  made  using  total  harmonic  distortion  calculations, 
to  establish  the  reliability  of  using  the  hard-switched  model  to  perform  control  synthesis 
for  the  ARCP. 

Several  control  algorithms  are  tested  through  simulation  and  the  results  analyzed. 
The  advantages  of  performing  control  in  the  synchronous  vice  stationary  reference  frame 
are  shown.  Testing  on  a  reduced-scale  circuit  model  using  a  digital  signal  processing 
system  (dSPACE)  to  implement  control  algorithms  is  used  to  validate  the  control 
algorithm  simulations.  A  new  method  of  waveform  modulation,  Space  Vector  Control,  is 
introduced  and  compared  with  conventional  methods.  Finally  the  operation  of  a 
prototype  ARCP  unit  is  discussed,  and  recommendations  for  improvements  in  future 
designs  are  presented. 
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/.  INTRODUCTION 


Current  naval  shipboard  power  requirements  are  satisfied  using  classical  AC 
power  distribution  employing  transformers  and  three-phase  transmission  cables.  As  the 
Navy  moves  toward  the  "more  electric"  initiative  for  propulsion,  electromechanical 
actuation  and  high  power  weapons  systems  on  ships,  submarines  and  aircraft,  there  is  a 
need  to  capitalize  on  electric  power  developments  [Ref.  1].  Improvements  in 
semiconductor  technology,  specifically  increased  current  density  capabilities,  have  made 
advanced  shipboard  power  distribution  schemes  using  semiconductor  devices  realizable. 
Critical  requirements  for  Naval  electrical  systems  include  reduction  in  size  and  cost, 
increase  in  power  density,  and  the  inclusion  of  intelligence  within  power  devices  that  can 
provide  electric  power  control,  conditioning,  and  very  high-speed  switching  [Ref.  1], 

A.  PROPOSED  SHIPBOARD  DC  POWER  SYSTEM 

In  order  to  best  take  advantage  of  the  power  electronics  revolution,  a  new 
architecture  for  shipboard  power  distribution  has  been  proposed  [Ref.  2].  In  this  new 
distribution  scheme  the  main  feeders  are  DC.  The  ship  is  divided  into  zones,  each  zone 
containing  common  energy  conversion  devices  fed  from  the  DC  busses.  As  shown  in 
Figure  1-1,  DC  is  distributed  via  port  and  starboard  busses  form  the  source(s)  into  the 
separate  zones  [Ref.  3].  Each  zone  contains  a  number  of  Ship  Service  Converter  (SSCM) 
and  Inverter  Modules  (SSIM).  The  SSCM  is  used  in  each  zone  to  step-down  the 
distribution  bus  voltage  to  a  regulated  level  for  use  in  the  zone.  In  this  way  the  SSCM 
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inserts  intelligence  into  the  system  by  acting  as  a  buffer,  preregulator  and  fault  protection 
for  each  zone.  Electric  loads  within  the  zone  are  fed  by  either  SSCM’s  or  SSIM’s 
depending  on  the  load’s  requirement  for  DC  or  AC  power. 

The  focus  of  this  thesis  is  on  the  SSIM.  Traditionally  DC-to-AC  converters  have 
been  uni-directional  power  flow,  non-resonant  type  devices.  Allowing  bi-directional 
(DC-to-AC  and  AC-to-DC)  power  transfer  improves  the  efficiency  of  the  system  since 
energy  that  would  be  unrecoverable  in  a  uni-directional  system,  such  as  during  motor 
braking,  can  now  be  captured  and  returned  to  the  system.  Resonant  converters  have 
lower  losses  and  emit  less  electromagnetic  interference.  It  is  proposed  to  use  a  resonant 
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converter,  specifically  an  Auxiliary  Resonant  Commutated  Pole  (ARCP)  converter,  as  a 
SSIM. 

The  remainder  of  this  chapter  discusses  the  inherent  advantages  and  disadvantages 
of  resonant  converters.  In  Chapter  II  the  basic  operation  of  the  ARCP  is  presented  and 
the  equations  governing  circuit  behavior  are  developed.  Chapter  III  deals  with  simulation 
of  single-phase  and  three-phase  converters  and  the  comparison  of  hard-switched  (non¬ 
resonant)  and  ARCP  converter  simulations.  Control  of  converters  with  various  loads  and 
with  various  control  designs  is  investigated  in  Chapter  IV.  A  new  pulse-width- 
modulation  technique  (Space  Vector  Control)  which  provides  advantages  over  sine- 
triangle  control  is  presented  in  Chapter  V.  Testing  of  control  algorithms  using  a  hard- 
switched  converter  and  digital-signal-processing  hardware  is  covered  in  Chapter  VI.  The 
operation  and  physical  characteristics  of  a  prototype  ARCP  converter  are  discussed  in 
Chapter  VII,  together  with  suggestions  for  improving  the  circuit  implementation.  The 
final  chapter,  contains  a  summary  of  the  research  work,  notable  conclusions  and 
recommendations  for  future  work. 

B.  ADVANTAGES  OF  RESONANT  CONVERTERS 

Examining  a  typical  switch  in  a  power  converter  provides  insight  into  the 
advantages  of  resonant  converters.  In  the  ideal  situation  with  a  switch  closed  and 
conducting  a  current  ILoad,  the  voltage  across  the  switch  is  zero  volts.  When  the  switch  is 
gated  off  the  current  falls  off  linearly  to  zero  and  the  voltage  rises  linearly  to  some  final 
value  which  is  based  on  the  circuit  topology.  For  purposes  of  illustration,  the  final  value 
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of  the  voltage  across  the  switch  will  be  called  V^..  The  power  loss  during  switching  will 
be  the  voltage  times  the  current.  Figure  1-2  shows  idealized  waveforms  for  turn  on  and 


Figure  1-2  -  Power  Loss  During  Switching 


turn  off  in  the  case  of  a  resistive  load.  The  current  rise  and  fall  times  are  controlled  by 
the  type  of  devices  and  the  drive  current  used.  For  a  typical  150  amp/1200  volt  Insulated 
Gate  Bipolar  Transistor  (IGBT)  (POWEREX  CM150DY-24E),  the  rise  and  fall  times  are 
400  nsec  [Ref.  4].  Most  IGBT  rise  and  fall  times  are  on  this  order  of  magnitude,  with 
larger  current  ratings  having  longer  rise  and  fall  times. 

Each  time  the  switch  changes  state  a  specified  amount  of  energy  is  lost.  Some  of 
the  energy  is  dissipated  as  heat  and  some  goes  to  the  generation  of  electro-magnetic- 
interference  (EMI).  The  short  rise  and  fall  times  cause  significant  energy  to  be  generated 
in  the  radio  frequency  (RF)  region  [Ref.  5].  The  more  frequently  the  switch  changes  state 
the  more  energy  is  lost,  thus  the  power  loss  is  proportional  to  the  switching  frequency. 
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For  high-frequency  pulse-width-modulation  techniques  switching  losses  can  be  quite 
significant.  This  not  only  lowers  the  efficiency  of  the  overall  system,  it  also  places  stress 
on  the  devices  used  for  switches.  The  ability  of  a  semiconductor  device  to  remove  heat  is 
limited.  As  the  heat  load  increases,  temperature  rises  which,  in  turn,  degrades 
performance. 

To  reduce  the  stress  on  the  devices  used  for  switches  and  limit  EMI,  snubbers  can 
be  used.  The  goal  of  a  snubber  during  device  tum-off  is  to  provide  zero  voltage  while  the 
current  dies  off.  During  device  tum-on  the  snubber  tries  to  maintain  a  zero  current  while 
the  voltage  drops  off  to  zero.  In  this  way  current  and  voltage  are  not  present  in  the  switch 
simultaneously  and  thus  ideally  no  power  is  lost  in  the  switch.  Additionally  snubbers 
limit  the  di/dt  and  dv/dt  of  the  switches  which  reduces  the  EMI  emissions. 

Figure  1-3  shows  a  typical  turn-off  snubber  attached  to  a  transistor  [Ref.  5].  To 
analyze  the  snubber  it  is  assumed  that  the  transistor  is  initially  gated,  the  voltage  across  it 
is  zero  and  it  is  conducting  a  current  ILoad.  When  the  transistor  is  turned  off  the  current 
will  drop  at  a  rate  determined  by  the  fall  time  of  the  device,  and  the  voltage  will  rise  as 
the  snubber  capacitor  begins  to  charge.  With  the  snubber  capacitor  in  the  circuit,  the 


Figure  1-3  -  Turn-Off  Snubber 
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voltage  rise  is  slower  than  the  current  fall  off  (depending  on  the  value  of  the  capacitor) 
and  the  energy  lost  in  the  switch  is  smaller  than  without  the  snubber.  Figure  1-4 
illustrates  the  waveforms  for  transistor  turn  off. 


The  energy  stored  in  the  capacitor  after  turn  off  is  dissipated  in  the  resistor  when 
the  switch  is  turned  on,  thus  the  snubber  does  not  improve  the  overall  efficiency  of  the 
circuit  but  simply  limits  the  stress  on  the  transistor. 

A  turn-on  snubber  [Ref.  5]  and  the  waveforms  associated  with  it  are  shown  in 
Figures  1-5  and  1-6.  The  snubber  inductor  limits  the  current  build  up  until  the  voltage 
across  the  switch  has  died  down.  This  limits  the  current  loss  in  the  switch.  However, 
during  the  off  time  of  the  switch  the  energy  stored  in  the  inductor  is  dissipated  in  the 
snubber  resistor,  so  the  overall  efficiency  of  the  circuit  is  not  improved.  Since  the 


6 


4 


inductor  is  required  to  carry  the  entire  load  current,  it  is  expensive  and  is  seldom  used 
alone. 


Figure  1-6  -  Turn-On  Snubber  Waveforms 


The  previous  discussion  has  shown  that  switch  stress  and  EMI  can  be  reduced  by 
the  use  of  snubbers  but  that  the  overall  efficiency  of  the  circuit  is  not  improved.  To 
improve  the  efficiency  of  the  overall  system  the  switching  losses  must  be  reduced  and  the 
energy  saved  in  reduced  switching  losses  must  be  somehow  returned  to  the  system. 
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Resonant  converters  attempt  to  accomplish  this  by  using  some  type  of  LC  resonance.  As 
the  voltage  and  current  resonate,  switches  are  turned  on  and  off  at  either  zero  voltage 
(zero  voltage  switching)  or  zero  current  (zero  current  switching),  thus  avoiding  power 
loss.  Since  the  circuits  resonate,  the  energy  used  to  drive  the  circuit  to  a  zero  state  is 
returned  to  the  system  avoiding  energy  loss.  Electro-magnetic  interference  is  also 
eliminated  since  no  power  is  lost  during  switching. 

Resonant  converters  are  generally  more  complex  than  non-resonant  converters. 
As  a  result,  resonant  converters  are  not  normally  used  in  low-power  or  low-frequency 
applications  where  the  added  cost  and  complexity  outweigh  the  benefits  of  improved 
efficiency,  lower  switch  stresses  and  reduced  EMI.  Higher  power  and  frequency 
converters  are  more  likely  to  be  of  the  resonant  type  since  the  efficiency  improvements 
begin  to  become  significant.  In  general  the  decision  to  use  a  resonant  converter  is  a 
tradeoff  between  complexity  and  efficiency  and  between  increased  initial  cost  and 
continued  savings. 
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II.  BASIC  OPERATION  OF  ARCP  CIRCUIT 


This  chapter  contains  a  detailed  analysis  of  the  basic  operation  of  the  ARCP 
converter.  A  single-phase  converter  is  used  to  explain  the  dynamics  involved  in 
commutation  and  how  ARCP  commutation  differs  from  that  occurring  in  a  hard-switched 
inverter.  Three  different  commutation  situations  are  explained  in  detail,  and  the 
extension  to  three-phase  power  applications  is  demonstrated.  The  majority  of  the 
material  presented  on  the  ARCP  is  derived  from  DeDoncker  [Ref.  6]. 

A.  BASIC  COMMUTATION 

Operation  of  the  Auxiliary  Resonant  Commutated  Pole  (ARCP)  inverter  is  similar 
to  that  of  a  hard-switched  inverter.  The  difference  lies  in  the  commutation  between 
states.  For  the  hard-switched  inverter,  commutation  involves  opening  one  switch  and 
closing  another.  Thus,  the  inverter  leg  output  voltage  changes  instantaneously,  in  the 
ideal  case,  from  VDC  to  zero,  or  vice-versa.  In  the  ARCP  in  order  to  reduce  switching 
losses,  commutation  is  accomplished  through  additional  auxiliary  circuitry  in  a  finite 
amount  of  time. 

The  basic  single-phase  ARCP  circuit  is  shown  in  Figure  2-1 .  It  consists  of  two 
main  switches,  SI  and  S2,  with  associated  free-wheeling  diodes,  D1  and  D2;  two 
auxiliary  switches,  A1  and  A2;  two  DC  supply  capacitors,  labeled  2Cdc;  two  resonant 
capacitors,  labeled  Cr/2;  a  resonant  inductor,  Lr;  and  a  filter  inductor  Lf.  The  main 
switches  operate  much  the  same  as  in  a  classical  hard-switched  inverter.  When  the  output 
voltage,  Vf,  is  desired  to  be  high  switch  SI  is  closed  and  when  the  output  voltage  is 
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Figure  2-1  -  Single-Phase  ARCP 


desired  to  be  zero  switch  S2  is  closed.  Switches  A1  and  A2  along  with  the  resonant 
inductor,  Lr,  make-up  the  auxiliary  circuit.  Auxiliary  switches  A1  and  A2  are  not 
required  to  turn  off  with  current  present  and  can  be  SCR’s,  GTO’s,  ZTO’s  or  MCT’s. 

The  main  switches  SI  and  S2  must  be  turn-off  devices  and  can  be  GTO’s,  MCT’s  or 
IGBT’s. 

The  auxiliary  circuit  is  involved  in  commutation  as  is  explained  in  the  following 
sections.  The  resonant  inductor  along  with  the  resonant  capacitors  determine  the  resonant 
frequency  of  the  converter  and  thus  define  the  speed  at  which  the  converter  can  be 
switched.  The  DC  supply  capacitors  are  large  capacitors  which  provide  a  voltage  equal  to 
half  of  the  source  voltage  at  the  center  of  the  two  capacitors.  If  the  negative  rail  of  the 
load  is  connected  to  this  midpoint,  an  AC  voltage  with  no  DC  offset  is  applied  to  the 
load.  The  filter  inductor  is  normally  several  orders  of  magnitude  larger  than  the  resonant 
inductor  and  helps  maintain  the  load  current  constant  during  commutation. 
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In  order  to  implement  an  ARCP  inverter  in  the  laboratory,  additional  circuitry  is 
required  to  monitor  signals,  gate-on  the  switches  and  provide  necessary  protection.  This 
additional  circuitry  is  not  considered  in  the  forthcoming  discussion  to  focus  attention  on 
the  commutation  process.  The  auxiliary  circuit  and  all  of  its  attendant  additional  circuitry 
substantially  increase  the  parts  count  of  an  ARCP  inverter  as  compared  to  a  hard- 
switched  inverter  and  is  one  of  the  drawbacks  that  must  be  considered. 

In  a  classical  hard-switched  inverter,  shown  in  Figure  2-2,  commutation  involves 
simply  turning  off  one  switch  and  turning  on  another.  In  an  ideal  circuit  the  turn-on  and 
turn-off  can  take  place  simultaneously;  however,  in  a  real  circuit  the  time  delays  involved 
can  result  it  a  direct  short  across  the  source.  As  a  consequence  a  time  delay  is  typically 
implemented  into  the  commutation  process.  Even  with  the  time  delay  issue,  the 
commutation  and  circuitry  involved  are  much  simpler  than  that  associated  with  the 
ARCP. 


Figure  2-2  -  Single-Phase  Hard-Switched  Converter 
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There  are  three  basic  cases  for  commutation  in  the  ARCP  inverter.  In  a  steady- 
state  switching  state  either  the  upper  or  lower  switch  (SI  or  S2)  will  be  closed  at  a  given 
instant,  and  either  the  diode  (D1  or  D2)  or  the  switch  will  be  carrying  the  current 
(depending  on  the  polarity  of  the  load  current).  When  it  is  desired  to  change  states  from 
the  upper  to  the  lower  switch,  the  direction  of  current  flow  and  hence  the  device  in 
conduction  must  be  determined.  For  example,  if  the  output  voltage,  Vf,  is  V0C  and  the 
current  flow  is  into  the  inverter,  then  diode  D1  is  conducting  the  load  current. 

The  three  commutation  cases  are  explained  in  detail  in  the  following  sections. 

The  first  of  the  three  cases  occurs  when  either  diode  D1  or  D2  is  conducting.  In  this  case 
commutation  is  accomplished  using  the  auxiliary  circuit  to  add  additional  energy  into  the 
circuit.  Additional  energy  is  required  to  overcome  the  losses  in  the  circuit.  For  example, 
if  no  additional  energy  were  added  to  the  circuit  during  commutation  from  Vf  =  0  to  Vf  = 
Voc,  the  output  voltage  would  never  quite  reach  the  upper  rail  voltage  while  resonating 
and  the  upper  switch  would  have  to  be  turned  on  with  a  non-zero  voltage  present.  The 
second  and  third  cases  occur  when  either  of  the  main  switches  are  conducting.  In  the 
second  case  the  magnitude  of  the  current  is  below  some  threshold  level,  and  the  auxiliary 
circuit  is  used  in  commutation.  This  threshold  level  is  chosen  to  limit  the  peak  current 
seen  in  the  main  switches.  The  third  case  occurs  when  the  magnitude  of  the  current  is 
above  a  threshold  level.  In  this  final  case  the  auxiliary  circuit  is  not  used  in  commutation; 
the  main  switch  is  simply  turned  off  and  after  the  voltage  swings  to  the  opposite  state  the 
other  main  switch  is  gated  on. 


12 


In  summary,  commutation  from  a  diode  always  involves  the  auxiliary  circuit  and 
commutation  from  a  switch  does  not  always  involve  the  auxiliary  circuit. 

B.  CASE  1  -  COMMUTATION  FROM  DIODE 

The  starting  point  for  Case  1  is  shown  in  Figure  2-3.  Initially  diode  D2  is 
assumed  to  be  conducting  the  full  load  current  and  the  output  voltage,  Vf,  is  zero.  Switch 
S2  remains  gated  on.  A  DC  source  voltage  is  continuously  applied  across  the  two 
capacitors,  2Cdc.  The  mid-point  of  these  capacitors  is  at  a  voltage  VDC/2.  It  is  assumed 
that  the  commutation  process  occurs  at  a  much  faster  rate  than  the  fundamental  of  the 
output  quantities  is  changing.  Thus,  it  is  reasonable  to  assume  that  the  load  current  is 
effectively  constant  during  the  commutation  interval. 


2Cdc  = 

2Cdc  = 

Vdc 

01 

IE 

'si  _ 

_ 

_  Cr/2 

iLoad 

0 

bi 

D2 

i 

S2  ” 

Lf  + 

=  Cr/2  Vf=0 

Figure  2-3  -  Initial  Condition,  D2  Conducting  S2  Gated  On 


To  start  the  commutation  process,  switch  A2  is  gated  on.  With  the  voltage  drop 
across  conducting  devices  assumed  to  be  zero,  this  introduces  a  forcing  potential  of  Vdc/2 
across  the  filter  inductor  Lr ,  as  shown  in  Figure  2-4.  The  ramp-up  phase  begins  as  the 
auxiliary  current,  ir,  increases  at  a  linear  rate  of  VDC/(2Lr),  dictated  by  the  relationship 
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Figure  2-4  -  Ramp-up  A2  -  Gated  On 


V  =  L  — .  Switch  S2  remains  gated  on  and  the  auxiliary  current  begins  to  displace  the 
dt 

load  current  flowing  through  D2. 

When  ir  exceeds  zload,  assumed  to  be  approximately  constant  throughout 
commutation,  the  boost  phase,  shown  in  Figure  2-5,  begins.  The  amount  of  auxiliary 

current  in  excess  of  the  load  current  flows  through  switch  S2,  and  is  defined  as  zboost  =  it  - 

Li2  _ 

zload.  This  boost  current  adds  energy  of  an  amount  equal  to  to  the  circuit.  This 

additional  energy  is  used  to  overcome  the  losses  in  the  circuit  during  the  resonant 
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Figure  2-5  -  Boost,  Ir  exceeds  IL  and  1^,  flows  through  S2 
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commutation  phase.  The  boost  current  adds  to  the  magnitude  of  the  voltage  resonance 
such  that  the  overall  magnitude  is  large  enough  to  reach  the  upper  rail  voltage  during 
commutation.  Since  iboosl  is  increasing  at  a  linear  rate,  the  boost  phase  can  be  controlled 
by  a  simple  time  delay.  At  this  point,  the  auxiliary  circuit  is  carrying  the  entire  load 
current  plus  the  boost  current.  Also  of  note  in  this  boost  interval  is  that  the  recovery 
current  in  D2  is  in  the  direction  of  the  boost  current,  and  thus  adds  additional  boost 
energy  which  assists  in  commutation. 

When  the  required  boost  energy  has  been  added  (as  specified  by  the  control),  S2 
is  gated  off  beginning  the  resonant  commutation  phase.  The  current  in  S2  at  shut-off  is 
diverted  into  the  resonant  capacitors,  Cr  /  2 ,  as  shown  in  Figure  2-6.  Some  switching 
losses  will  occur  in  switch  S2  at  turn-off  but  the  resonant  capacitors  will  act  as  snubbers 
to  limit  these  losses.  The  auxiliary  current  and  output  voltage  now  resonate  at  a 
frequency  determined  by  the  values  Lr  and  Cr .  The  auxiliary  current  will  rise  to  a  peak 
and  then  drop  while  the  output  voltage  swings  towards  the  upper  rail  voltage. 


Figure  2-6  -  Resonant,  S2  Gated  Off 
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In  analyzing  the  resonant  commutation  interval,  the  load  current  can  be 
considered  constant,  since  Lf  »  Lr,  and  the  front-end  capacitors  can  be  modeled  as 

voltage  sources  with  a  value  of  VDC/2,  since  2Cdc  »  Cr  /  2 .  Details  of  the  analysis  of 
the  resonant  phase  are  provided  in  Appendix  A.  Assuming  ideal  circuit  elements,  the 
maximum  current  in  the  auxiliary  circuit  can  be  approximated  by: 


h,max 


(Eq.  2-1) 


This  value  is  used  to  size  the  auxiliary  branch  devices. 

When  the  output  voltage  becomes  equal  to  the  DC  input  voltage,  zero  voltage 
exists  across  SI .  At  this  point  SI  is  gated  on  with  no  switching  losses,  clamping  the 
output  voltage  at  the  high  state  (Figure  2-7). 


Figure  2-7  -  Clamp,  SI  Gated  On 

When  SI  is  closed  in  order  to  clamp  the  output  voltage  at  the  upper  rail,  a  forcing 

potential  of  VDC/2  is  locked  across  the  resonant  inductor  Lr  as  shown  in  Figure  2-8.  This 

-V 

causes  the  auxiliary  current  to  drop  off  linearly  at  a  rate  of — — .  As  the  auxiliary 


current  drops  off,  load  current  shifts  to  switch  SI . 
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Figure  2-8  -  Ramp-Down,  SI  Picks-up  Load,  Auxiliary  Current  Falls  Off 

As  auxiliary  current  falls  to  zero,  A2  is  gated  off  with  no  current  flowing  through 
it;  therefore,  there  are  no  switching-off  losses.  In  actuality  A2  will  be  gated  off  in 
anticipation  of  the  current  crossing  zero  so  that  the  crossover  point  can  be  accurately 
predicted.  The  diode  shunting  A1  will  stop  the  current  from  resonating  back  negative, 
but  some  losses  will  occur  if  the  diode  blocks  the  current  due  to  reverse  recovery  losses. 
A  snubber  is  normally  added  across  the  auxiliary  devices  to  eliminate  the  possibility  of 
damaging  the  switches  due  to  a  small  timing  error. 


Figure  2-9  -  Completion,  A2  Gated  Off 
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Figure  2-10  shows  the  output  voltage  and  auxiliary  circuit  current  for  Case  1 
commutation.  Initially  the  auxiliary  current  is  zero,  the  output  voltage  is  zero  and  diode 
D2  is  conducting  the  load  current.  The  devices  in  conduction  are  indicated  below  the 


voltage  graph  in  Figure  2-10.  When  the  auxiliary  device  A2  is  switched  on,  the  ramp-up 
phase  begins  and  the  auxiliary  current  begins  to  rise,  and  both  D2  and  A2  are  conducting 
As  the  auxiliary  current  increases  above  the  level  of  the  load  current,  diode  D2  stops 
conducting  and  current  begins  flowing  through  the  lower  main  switch  S2  ,  initiating  the 
boost  phase  of  commutation.  After  the  boost  current  reaches  the  desired  level  S2  is 
switched  off  commencing  the  resonant  phase.  Throughout  the  resonant  phase  auxiliary 
switch  A2  is  conducting  the  entire  load  current  in  conjunction  with  the  auxiliary  current. 
The  output  voltage  swings  high  and  the  auxiliary  current  peaks  and  begins  to  swing  low. 


« 
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When  the  output  voltage  reaches  the  DC  voltage  the  upper  main  switch  (SI)  is  gated  on. 
In  the  diagram  this  occurs  just  as  the  auxiliary  circuit  current  drops  to  the  boost  plus  load 
current  level.  In  an  actual  circuit  the  voltage  would  not  reach  this  level  until  slightly  after 
the  auxiliary  current  passed  below  the  boost  plus  auxiliary  current  level  due  to  losses. 
Figure  2-10  is  generated  assuming  no  losses.  Once  SI  is  gated  on  the  auxiliary  current 
drops  off  at  a  linear  rate.  When  the  auxiliary  current  drops  below  the  level  of  the  load 
current,  SI  begins  to  conduct.  As  the  auxiliary  current  drops  to  zero,  A2  is  gated  off 
completing  the  commutation. 

C.  CASE  2  -  SWITCH  CONDUCTING,  CURRENT  BELOW  THRESHOLD 

Commutation  for  Case  2  is  similar  to  Case  1  except  that  there  are  no  ramp-up  or 
ramp-down  phases  involved. 

Initially  switch  SI  is  assumed  to  be  conducting  the  full  load  current,  and  the 
output  voltage  is  equal  to  VDC,  as  illustrated  in  Figure  2-11. 


Figure  2-11  -  Initial  Condition,  SI  Conducting 
Switch  A1  is  gated  on  to  start  the  commutation  process.  This  introduces  a  forcing 

potential  of  VDC/2  across  the  filter  inductor  Lr .  The  current  flow  through  the  auxiliary 
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circuit  is  in  a  direction  to  increase  the  current  flow  through  SI  as  shown  in  Figure  2-12. 
Thus  no  ramp-up  phase  is  needed  as  the  auxiliary  current  immediately  initiates  a  boost  of 
the  load  current. 


Figure  2-12  -  Boost,  Ir  Flows  through  SI 


When  the  boost  current  reaches  the  desired  level,  the  resonant  commutation  phase 
is  started  by  switching  SI  off.  The  output  voltage  begins  to  fall  while  the  current  through 
SI  is  diverted  to  the  resonant  capacitors.  This  action  is  depicted  in  Figure  2-13.  Ideally, 
the  auxiliary  current  consists  of  a  half-sinusoid  superimposed  upon  the  boost  current 
level,  and  the  output  voltage  falls  as  the  sum  of  cosine  and  ramp  waveforms. 
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Figure  2-13  -  Resonant,  SI  Gated  Off 
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As  the  output  voltage  continues  to  fall  it  will  attempt  to  drop  below  zero.  At  this 
point  diode  D2  becomes  forward  biased  and  begins  to  pick  up  load  current.  Switch  S2  is 
turned  on  at  this  point  to  clamp  the  voltage  to  zero.  Energy  remaining  in  the  resonant 
inductor  is  returned  to  the  DC  bus  capacitors  as  the  potential  across  the  inductor  is  now  of 
opposite  polarity  as  at  the  initiation  of  the  resonant  phase. 


Figure  2-14  -  Clamp,  S2  Gated  On 

When  the  auxiliary  current  falls  to  zero  A1  is  gated  off.  Diode  D2  is  now 
conducting  all  of  the  load  current  and  commutation  is  complete. 


Figure  2-15  -  Completion,  A1  Gated  Off 


Figure  2-16  illustrates  the  output  voltage  and  auxiliary  current  waveforms  for 
Case  2  commutation.  Initially  the  output  voltage  is  high,  auxiliary  current  is  zero  and  the 
upper  main  switch  SI  is  conducting  all  of  the  load  current.  Auxiliary  switch  A1  is  gated 
on  to  start  the  commutation  process;  this  immediately  initiates  the  boost  phase.  The 


boost  current  increases  to  a  specified  level,  then  the  upper  main  switch  SI  is  gated  off 
initiating  the  resonant  phase.  The  output  voltage  swings  low  and  the  auxiliary  current 
peaks  in  magnitude  and  begins  to  decrease.  When  the  output  voltage  reaches  zero,  diode 
D2  begins  conducting  and  S2  is  gated  on  with  a  zero  voltage  differential.  When  the 
auxiliary  current  drops  to  zero,  A1  is  gated  off  completing  the  commutation.  It  can  be 
noted  that  in  Figure  2-16  the  output  voltage  reaches  zero  just  as  the  auxiliary  current 
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drops  to  the  boost  level.  If  losses  are  considered,  the  time  at  which  the  output  voltage 
reaches  zero  would  be  shifted  past  this  point. 


D.  CASE  3  -  SWITCH  CONDUCTING,  CURRENT  ABOVE  THRESHOLD 

When  the  load  current  exceeds  a  pre-defined  threshold  level,  the  auxiliary  circuit 
is  not  used  when  commutating  from  a  conducting  switch.  To  determine  the  threshold 
current  level,  the  maximum  commutation  time,  Tmax ,  (i.e.  the  time  required  for  a  Case  1 
commutation  at  full-load  current)  must  first  be  computed.  It  is  desired  to  maintain  the 
Case  3  commutation  time  less  than  or  equal  to  the  maximum  commutation  time  for  Case 
1.  In  Case  3  the  change  in  output  voltage  is  governed  by  Equation  2-2. 


dVoul  _  j f_Load 
dt  ~  Cr 


(Eq.  2-2) 


The  larger  the  load  current,  the  faster  the  output  voltage  will  change;  therefore,  a  small 
output  current  will  cause  commutation  to  take  longer.  To  determine  the  smallest  current 
at  which  the  Case  3  commutation  time  will  be  less  than  or  equal  to  Tmm ,  Equation  2-2  is 
used  with  dVoul  =  FDC  and  dt  =  Tmax .  This  current,  defined  as  the  threshold  current,  takes 
on  the  following  value 

I,,  -  ^  (Eq.  2-3) 

1  max 

In  component  selection  the  threshold  current  should  be  selected  to  be  sufficiently 
below  the  peak  load  current.  Thus,  IA  +  I^oA  <  1^.^.  If  this  criteria  is  met,  then  the 
primary  switches  do  not  have  to  turn  off  any  higher  current  than  if  a  hard-switched 
converter  were  used. 
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Case  3  commutation  is  described  as  follows.  Initially  switch  SI  is  assumed  to  be 
conducting  the  full  load  current,  and  the  output  voltage  is  equal  to  VDC.  This  is  the  same 
initial  condition  as  given  for  Case  2. 


Figure  2-17  -  Initial  Condition,  SI  Conducting 


Switch  SI  is  gated  off  to  begin  commutation.  This  causes  the  load  current  to  be 

shunted  into  the  two  resonant  capacitors  as  shown  in  Figure  2-18.  The  output  voltage 

dV 

falls  at  a  rate  proportional  to  the  output  current  following  the  relationship  /'  =  C — . 

dt 


When  the  output  current  attempts  to  fall  below  zero,  diode  D2  becomes  forward  biased 


Figure  2-18  -  Commutation,  SI  Gated  Off 


24 


and  picks  up  the  load  current.  Switch  S2  is  gated  on  at  this  point,  with  no  losses, 
completing  the  commutation. 


Figure  2-19  -  Completion,  D2  Begins  Conducting,  S2  Gated  On 


Figure  2-20  illustrates  the  output  voltage  and  auxiliary  current  waveforms  for 
Case  3  commutation.  Initially  the  output  voltage  is  high,  auxiliary  current  is  zero  and  the 
upper  main  switch  SI  is  conducting  all  of  the  load  current.  SI  is  gated  off,  dumping  the 
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load  current  into  the  resonant  capacitors.  The  output  voltage  falls  linearly  following  the 

dV 

constitutive  capacitor  relationship  ihad  =  Cr  — — .  Once  the  voltage  falls  to 

dt 

zero,  diode  D2  begins  conducting  and  S2  is  gated  on,  completing  the  commutation. 

E.  THREE-PHASE  APPLICATIONS 

To  integrate  the  ARCP  into  a  three-phase  power  system  is  a  relatively  simple  task 
of  connecting  three  single-phase  ARCP  units  into  a  correct  configuration.  A  common 
DC  voltage  bus  is  used  for  all  three  units  and  the  AC  sides  are  connected  in  the  desired 
configuration,  either  Y  or  A.  A  Y-connected  configuration  is  shown  in  Figure  2-21. 


Figure  2-21  -  Three-Phase  ARCP,  Y-Connected  Load 


One  consequence  of  the  above  configuration  is  that  the  neutral  point  of  the  Y- 
connected  load  will  be  at  the  same  potential  as  the  average  DC  input  voltage.  To  avoid 
an  offset  at  the  neutral  point  the  DC  voltage  can  be  made  bi-polar  with  the  positive  and 
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negative  voltages  of  equal  value,  thus  the  average  value  will  be  zero.  This  problem  is 
avoided  in  a  A-configuration  since  no  neutral  point  exists. 

Power  flow  can  be  in  either  direction,  DC  to  AC  or  AC  to  DC.  No  additional 
modifications  are  necessary  in  the  ARCP  to  accommodate  the  bi-directional  power  flow. 
However,  the  load  and  source  must  be  capable  of  sourcing  and  sinking  power. 
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III.  SIMULATION 


This  chapter  documents  the  approach  and  issues  related  to  the  simulation  of  a 
number  of  structures  considered  in  this  work.  Initially  a  digital-based  simulation  of  a 
single-phase  hard-switched  inverter  is  developed  with  sine-triangle  voltage  control.  Next, 
the  details  required  to  model  a  single-phase  ARCP  converter  with  sine-triangle  control 
are  examined.  Finally  the  simulations  of  the  three-phase  versions  of  both  the  hard- 
switched  and  ARCP  converters  are  obtained.  Comparison  of  the  converters  using  total 
harmonic  distortion  calculations  shows  that  the  hard-switched  inverter  can  be  used  as  a 
practical  reduced-order  model  of  the  ARCP.  This  facilitates  faster  simulation  of 
controllers  and  the  convenient  investigation  of  system  load  transients.  In  addition,  the 
effect  of  a  minimum  switching  time  on  high-frequency  modulation  will  be  explored. 

A.  SIMULATION  LANGUAGE 

Simulations  described  in  this  chapter  are  performed  using  the  Advanced 
Continuous  Simulation  Language  (ACSL).  ACSL,  pronounced  ‘axle’,  is  a  Fortran  based 
language  which  uses  a  free  format  which  is  not  as  restrictive  as  Fortran.  Programs  are 
compiled  after  being  translated  into  Fortran  by  the  ACSL  translator.  Once  compiled,  the 
program  is  linked  with  the  ACSL  run-time  library  to  allow  user  interface. 

Models  are  constructed  from  state  space  descriptions  of  systems.  Discrete-event 
blocks  and  if-then-else  constructs  are  used  to  switch  from  one  set  of  equations  to  another. 
The  general  construction  of  a  program  is  shown  in  Figure  3-1  [Ref.  7].  The  derivative 


29 


PROGRAM 

INITIAL 

Statements  executed  before  the  run  begins. 

State  variables  do  not  contain  the  initial  conditions  yet. 

END 

DYNAMIC 

DERIVATIVE 

Statements  to  be  integrated  continuously. 

END 

DISCRETE 

Statements  executed  at  discrete  points  in  time. 

END 

Statements  executed  each  communication  interval 

END 

TERMINAL 

Statements  executed  after  the  runtime  terminates. 

END 

END 


Figure  3-1  -  ACSL  Program  Outline  (from  [Ref.  7]) 

section  is  sorted  so  that  as  long  as  the  inputs  and  outputs  of  various  subsystems  are 
compatible,  no  specific  sequential  ordering  of  statements  is  required  by  the  programmer. 

The  ACSL  language  was  selected  over  SIMULINK  for  several  reasons.  First, 
ACSL  programs  are  compiled  so  they  run  faster.  ACSL  also  handles  switching-type 
networks  far  more  easily  and  provides  more  options  for  solving  algebraic  loops.  Overall 
the  major  advantage  over  SIMULINK  is  the  speed  at  which  ACSL  provides  in  the 
simulation  of  complex  systems. 


B.  SINGLE-PHASE  HARD-SWITCHED  CONVERTER 

In  simulation  it  is  assumed  that  the  hard-switched  inverter  can  change  states 
instantaneously.  As  soon  as  one  switch  is  opened,  the  other  is  closed  with  no  delay.  In 
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practice  some  delay  must  be  introduced,  as  switches  do  not  turn  off  instantly,  to  avoid 
both  switches  being  shut  simultaneously  causing  a  short  across  the  source. 

A  sine-triangle  method  of  modulation  was  used  to  generate  the  control  signal 
[Ref.  8].  A  control  sinusoid  at  the  desired  output  frequency  is  generated  and 
superimposed  on  a  high-frequency  triangle  wave.  The  two  waves  are  compared  and 
when  the  sinusoid  is  greater  than  the  triangle  wave,  a  control  signal  to  switch  the 
particular  leg  of  the  inverter  to  the  high  state  is  generated.  When  the  control  sinusoid  is 
less  than  the  triangle  wave,  a  control  signal  for  the  low  state  is  generated.  Figure  3-2 
shows  graphically  how  sine-triangle  modulation  generates  an  output  voltage  scheme. 


The  output  voltage  waveform  is  a  square  wave  which  contains  many  harmonic 
components.  To  reduce  the  lower-order  harmonics  the  frequency  of  the  triangle  wave  is 
increased  in  relation  to  the  frequency  of  the  control  sinusoid.  In  the  simulations 
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considered,  the  frequency  of  the  control  sinusoid  is  60  Hz  while  the  frequency  of  the 
triangle  wave  is  6  kHz. 

Appendix  B  (pp  127-128)  lists  the  ACSL  code  for  the  single-phase  hard-switched 
converter.  A  resistor  with  an  LC  filter  is  the  assumed  load. 


C.  SINGLE-PHASE  ARCP  CONVERTER 

The  first  step  in  simulating  the  ARCP  was  to  generate  a  control  signal  for 
switching.  As  in  the  hard-switched  converter,  sine-triangle  modulation  was  considered. 

Establishing  the  switching  states  for  the  ARCP  is  much  more  complicated  than  in 
the  hard-switched  converter.  How  the  ARCP  changes  states  is  described  in  detail  in 
Chapter  II.  To  simplify  modeling  the  switching  dynamics  six  situations  encompassing  all 
possible  ways  of  changing  states  were  identified.  Three  situations  involved  a  positive 
load  current,  and  three  a  negative  load  current.  Each  situation  corresponds  to  a 
commutation  sequence  described  in  Chapter  II.  Table  3-1  shows  all  the  possible 
commutation  states  for  the  ARCP.  Figure  3-3  illustrates  the  decision  cycle  that  must  be 
implemented  in  the  ACSL  code.  The  load  current  sign  and  magnitude,  as  well  as  the 
initial  output  voltage,  are  the  quantities  which  dictate  the  commutation  case. 


Table  3-1  -  ARCP  Commutation  Possibilities 


State 

# 

Type  of 
Commutation 

Voltage 

Transition 

Load 

Current 

Above/Below 
Thres  Current 

1 

Case  1 

Lo  to  Hi 

Pos 

N/A 

2 

Case  2 

Hi  to  Lo 

Pos 

Below 

3 

Case  3 

Hi  to  Lo 

Pos 

Above 

4 

Case  2 

Lo  to  Hi 

Neg 

Below 

5 

Case  1 

Hi  to  Lo 

Neg 

N/A 

6 

Case  3 

Lo  to  Hi 

Neg 

Above 
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Figure  3-3  -  Flowchart  for  ARCP  Commutation 


Once  the  commutation  state  has  been  determined,  the  equations  governing  the 
circuit  behavior  must  be  found.  In  Chapter  II  the  Laplace  transform  solution  for  the 
resonant  portion  of  commutation  was  established.  For  ease  of  programming,  the  Laplace 
transform  solution  is  not  used.  Instead,  a  state  space  representation  is  used.  State 
variables  of  interest  are  the  auxiliary  circuit  current  (/,),  the  load  current  (iL ),  the  output 
voltage  (Vf),  and  the  load  voltage  ( VLoad ).  The  appropriate  state  space  representation 

depends  on  the  status  of  both  the  main  and  auxiliary  switches.  These  changes  are 
controlled  by  conditional  IF  statements  in  the  program.  The  derivatives  of  the  auxiliary 
current  and  output  voltage  for  various  switch  status  combinations  are  shown  in  Table  3-2. 
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Table  3-2  -  State  Variable  Equation  vs  Switch  Positions 


A1 

off 

off 

off 

off 

off 

on 

on 

on 

off 

A2 

off 

on 

on 

on 

off 

off 

off 

off 

off 

SI 

off 

off 

off 

on 

on 

on 

off 

off 

off 

S2 

on 

on 

off 

off 

off 

off 

off 

on 

off 

dL_ 

dt 

0 

2  L, 

V*-2  Vf 

2Lr 

0 

V*-2Vf 

2  Lr 

v«-iv, 

2  Lr 

0 

dVf 

dt 

0 

0 

*r-h 

Cr 

0 

0 

0 

m 

0 

-h 

Cr 

The  derivatives  of  load  current  and  load  voltage  are  independent  of  the  switch  status  and 
are  given  by: 


^ Load 

~dt  ~  Tf 


(Eq.  3-1) 


Load 


dV, 


Load 


R 


Load 


dt 


(Eq.  3-2) 


Appendix  B  (pp  129-133)  contains  the  ACSL  code  for  the  single-phase  ARCP. 


D.  THREE-PHASE  SIMULATIONS 

Expansion  to  three-phase  simulation  involved  combining  three  separate  single¬ 
phase  inverters  and  their  control  signals.  Load  dynamics  were  simulated  in  the  stationary 
reference  frame.  The  control  signals  consisted  of  three  sinusoids  of  equal  amplitude, 
phase  displaced  by  120  degrees.  This  resulted  in  a  standard  three-phase  output 
waveform.  Code  for  the  three-phase  hard-switched  and  ARCP  converters  is  listed  in 
Appendix  B  (ppl34  -  148). 
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E.  COMPARISON  OF  HARD-SWITCHED  AND  ARCP  CONVERTERS  USING 
TOTAL  HARMONIC  DISTORTION 


It  should  be  noted  that  less  simulation  code  is  required  for  implementing  the  hard- 
switched  converter  as  compared  to  the  ARCP.  This  is  due  to  the  ARCP  simulation 
modeling  each  commutation  in  detail,  while  the  hard-switched  converter  simulation  takes 
a  simplified  view  of  commutation.  Due  to  this  fact  the  hard-switched  simulation  runs 
much  faster  than  the  ARCP  program.  For  an  integration  time  step  of  O.lpsec  the  ARCP 
simulation  required  133.8  seconds  to  complete  a  .032  second  real-time  simulation  on  a 
SPARC  10  workstation.  The  hard-switched  simulation  took  63.4  seconds  to  run  an 
identically  structured  simulation.  Since  the  hard-switched  inverter  does  not  need  to 
model  the  commutations  in  detail,  the  integration  step  can  be  changed  to  a  larger  value 
without  loss  of  accuracy.  This  reduces  the  run  time  even  further.  Given  an  integration 
time  step  of  1  psec  the  same  hard-switched  simulation  run  took  only  6.6  seconds.  For 
this  reason,  it  is  desirable  to  substitute  the  hard-switched  model  for  the  ARCP  when 
evaluating  new  control  schemes  or  loads.  To  verify  that  this  type  of  substitution  is  valid  a 
comparison  of  output  waveforms  for  both  converters  was  undertaken. 

To  compare  the  two  models  the  total  harmonic  distortion  (THD)  of  each  model’s 
output  voltage  was  evaluated.  Total  harmonic  distortion  is  defined  as  the  square  root  of 
the  sum  of  the  squares  of  the  harmonic  components  divided  by  the  square  of  the 
fundamental,  or  in  equation  form: 


THD%  =  100  x . 


irz+r/s+KW 


s  1 


=  100x  2 

\  h*  1 


sh 


UsJ 


(Eq.  3-3) 


The  method  used  to  calculate  THD  is  described  in  detail  in  Appendix  C  [Ref.  9]. 
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A  simulation  was  run  at  a  switching  frequency  of  6  kHz  to  compare  the  models, 
harmonic  components  up  to  seventh-order  were  used  in  calculating  THD.  Results 
showed  the  THD  for  the  hard-switched  converter  was  0.8537%  and  for  the  ARCP  was 
0. 1745%.  These  are  on  the  same  order  of  magnitude  as  expected.  Total  harmonic 
distortion  for  the  ARCP  was  lower  than  for  the  hard-switched  converter.  This  also  was 
expected  since  the  waveforms  for  the  ARCP  were  less  sharp  at  the  transitions  from  high 
to  low  than  for  the  hard-switched  converter  which  had  square  transitions  and  thus  more 
harmonics.  For  this  reason  using  the  hard-switched  inverter  in  simulations  will  give  a 
conservative  result  with  respect  to  the  amount  of  distortion  created. 

F.  MINIMUM  SWITCHING  TIME 

Pulse-width-modulation  (PWM)  schemes  such  as  sine-triangle  modulation  can 
require  very  high-frequency  switching.  For  instance,  in  sine-triangle  modulation  with  a 
duty  cycle  of  95%  the  highest  frequency  component  will  be  20  times  the  fundamental.  At 
higher  switch  frequencies  this  can  become  a  problem  because  there  is  a  minimum  time 
required  for  the  ARCP  to  switch  states  which  may  be  longer  than  the  switching  time 
commanded  by  the  PWM  technique.  Minimum  switching  time  is  commonly  referred  to 
as  dead  time.  For  the  ARCP  dead  time  is  on  the  order  of  one  over  the  resonant  frequency. 
Hard-switched  converters  typically  have  a  short  dead  time  also,  to  allow  the  switches  to 
completely  turn  off. 

Dead  time  causes  additional  distortion  in  output  waveforms.  To  check  the  amount 
of  distortion,  a  simulation  using  a  hard-switched  converter  with  dead  time  was 
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implemented.  The  results  in  Table  3-3  show  a  dramatic  increase  in  THD  when  dead  time 


is  incorporated. 


Table  3-3  -  Simulation  Parameters 


Dead  Time 

25  psec 

Switching  Frequency 

6  kHz 

Maximum  Duty  Cycle 

95% 

Method  of  PWM 

sine-triangle 

Minimum  Switching 

Time  from  PWM 

8.3  psec 

Voltage  THD  -  with  dead  time 

1.010% 

Voltage  THD  -  w/o  dead  time 

0.039% 

Current  THD  -  with  dead  time 

3.584% 

Current  THD  -  w/o  dead  time 

0.080% 

In  an  ARCP  the  dead  time  can  be  approximated  by  Equation  3-4  [Ref.  6]. 

7’  O  T  \  * Load-max  +  ^ Boost)  ,  _  /  j  ^  /T-,  - 

Dead  =  2Lr - - - -  +  ^ LrCr  (Eq.  3-4) 

* DC 

Using  the  values  from  Table  3-4,  a  dead  time  value  of  3.44  psec  is  calculated. 


Table  3-4  -  ARCP  Parameters 


Resonant  Inductance  Lr 

2.9  pH 

Resonant  Capacitance  Cr 

0.3  pF 

Maximum  Load  Current 

20  Amps 

Boost  Current 

4  Amps 

So  for  a  maximum  duty  cycle  of  95  %  the  highest  frequency  triangle  wave  which  can  be 
used  in  sine-triangle  modulation  is  29  kHz.  If  the  dead  time  is  assumed  to  be  5  psec  for  a 
non-ideal  ARCP,  the  highest  triangle  wave  frequency  to  be  used  without  dead-time 
distortion  is  20  kHz.  Although  the  numbers  used  in  the  dead  time  simulation  vary  from 
those  of  the  ARCP,  the  simulation  provides  insight  to  the  analyst  into  the  amount  of 
distortion  which  dead  time  can  introduce. 
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Using  the  simulation  tools  presented  here  for  circuit  simulation,  the  issue  of 
control  can  be  addressed.  Various  control  algorithms  can  be  simulated  and  joined  with 
the  circuit  simulations  to  predict  system  behavior.  The  following  chapter  addresses  the 
control  issue  in  this  manner. 
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IV.  CONTROLLERS 


Open-loop  control  of  converters  can  produce  satisfactory  output  waveforms; 
however,  using  some  type  of  feedback  generally  improves  the  quality  of  the  output 
response.  Feedback  helps  make  the  system  response  less  sensitive  to  parameter 
variations  and  exogenous  disturbances.  In  addition,  feedback  is  introduced  to  improve 
the  transient  and  steady-state  response  or  make  reference  signal  tracking  automatic.  To 
control  the  output  of  a  three-phase  converter,  three  separate  control  schemes  are 
proposed.  The  first  uses  voltage  feedback  to  control  the  load  voltage,  the  second  uses 
current  feedback  to  control  the  load  current,  and  the  last  uses  a  combination  of  current 
and  voltage  feedback  to  control  the  load  voltage.  Reference  frame  theory  is  used  in  the 
control  algorithms,  and  the  advantages  of  implementation  in  different  reference  frames 
are  discussed. 

A.  VOLTAGE  FEEDBACK  CONTROLLER 

A  voltage  controller,  controlling  in  the  stationary  reference  frame,  is  presented 
and  analyzed  in  this  section.  The  controller  compares  the  actual  output  voltage  to  the 
desired  voltage,  generates  an  error  and  passes  this  error  in  parallel  through  proportional 
and  integral  amplifiers  to  generate  a  control  signal.  This  control  signal  is  used  to  derive 
the  gating  signals  for  the  switches  in  the  converter. 

For  the  purpose  of  simulation  a  circuit  shown  in  Figure  4-1,  consisting  of  a  three- 
phase  converter  driving  a  resistive  load  attached  to  an  LC  filter,  was  considered.  The 
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Figure  4-1  -  Converter  with  LRC  Load 


purpose  of  the  controller  is  to  regulate  the  voltage  across  the  filter  capacitor  and  load 
resistor,  Vout.  Possible  extensions  to  machine  loads  will  be  apparent. 

To  begin  analysis  of  the  circuit  the  three  phase  voltages  of  the  converter  are 
transformed  into  an  equivalent  wye-connected  source  as  shown  in  Figure  4-2.  For  a 
balanced  load  condition,  the  phase  voltages  are  determined  according  to  the  following 
relationships: 

CEq-4-l) 

(&!•«) 

where  Va.fj.c  s  represents  a  set  of  wye-connected  phase  voltages  and  Va.jj.c  g  represents 
phase-to-ground  voltages  at  the  converter  output.  The  voltages  Va.]j.c  g  are  known  once 
the  switch  status  of  the  inverter  is  specified.  Wye-connected  quantities  are  used  in  the 
remainder  of  the  analysis. 
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Figure  4-2  -  Equivalent  Wye-Connected  Source 


Equations  for  the  circuit  are  next  established  for  use  in  the  simulation.  The  state 
variables  include  the  inductor  current  and  the  capacitor  voltage.  Assuming  balanced 
conditions  the  center  of  the  wye-connected  source  and  the  common  point  of  the  resistive 
load  have  the  same  potential.  Writing  an  equation  from  one  common  point  to  the  other 
across  phase  A  gives: 


V  -  I  -2--V  =0 

y  as  x-y  ¥  Load -a  v 


(Eq.  4-4) 


where  ia  is  the  current  flowing  through  the  inductor  Lf.  Solving  for  the  derivative  of  the 
inductor  current  yields: 


di  V  -  V 

Ula  v  as  v  Load- a 


(Eq.  4-5). 


Summing  currents  at  the  common  node  between  the  inductor,  capacitor  and  load  resistor 
produces  the  following  equation: 


V  dV 

j  Load-a  ^  r  Load -a  q 

"  ’  <* 


(Eq.  4-6). 


Solving  for  the  derivative  of  capacitor  voltage  gives 
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(Eq.  4-7). 


l„  ~ 


Load-a 


Load- a  ^ Load 


dt  Cf 

Equations  for  phases  B  and  C  give  similar  results  with  the  subscript  ‘a’  replaced  by  ‘6’  or 
‘c’  as  appropriate. 

Next  the  state  equations  were  transformed  into  the  stationary  q-d-0  reference 
frame  using  the  arbitrary  reference  frame  transformation  matrix  [Ref.  10] 


K  =  — 
3 


cos#  cod 


0- 


sin#  sin  9  — 


2  n 

"3 
2  n 


\_ 

2 


2 


fa 

cos  0  +  — 
\  3  J 

sin^  +  yj 
2 


where 


(Eq.  4-8) 


6 =  jco  dt  +  90 


(Eq.  4-9) 


For  the  stationary  reference  frame,  a>  is  set  to  zero  and  90  is  selected  to  also  equal  zero 
for  convenience.  The  transformation  matrix  then  becomes 


K 


S 

s 


1  _I 

"2  2 

V3  V3 

2  2 
\  l 
2  2 


(Eq.  4-10). 


where  the  ‘5  ’  superscript  denotes  the  stationary  reference  frame. 

With  balanced  conditions  the  zero-sequence  variables  must  be  zero,  so  they  are 
not  considered  in  the  controller  analysis.  Zero-sequence  variables  may  be  used  to  check 
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for  balanced  operation,  but  for  the  simulation  it  is  assumed  that  balanced  conditions  are 
maintained. 


The  transformed  state  equations,  written  in  vector  form,  become 


<& 


qd  0 


Xs  _  Xs 

T  qdOs  v  qdOLoad 


dt 


(Eq.  4-11), 


dVs 

Uy  qdO  Load 

dt 


\S  - 
lqd  0 


r  qdOLoad 

h 


Load 


(Eq.  4-12). 


A  balanced  three-phase  set  of  desired  voltages  is  generated  and  transformed  into 
the  stationary  reference  frame  to  give  the  desired  q-d-0  voltages. 


V5  * 

q-Load 

>* 

a- Load 

ys  * 
v  d-Load 

=  K 

v * 

V  b- Load 

Vs  * 

\J  0- Load  J 

V* 

J  c-Load  _ 

(Eq.  4-13) 


Again  with  balanced  conditions  the  zero-sequence  voltage  will  be  zero,  so  it  is 
ignored.  Desired  q-d  voltages  are  compared  with  actual  q-d  voltages  to  provide  error 
signals,  one  error  signal  for  the  q-quantities  and  one  for  the  d-quantities.  These  error 
signals  are  next  sent  through  proportional-integral  controllers  which  produce  control 
signals  in  the  q-d  frame.  The  control  signals  are  transformed  back  into  the  a-b-c  frame  to 
give  a  commanded  control  signal  to  each  leg  of  the  converter.  Based  on  beating  these 
control  signals  up  against  a  high-frequency  carrier  triangle  waveform,  the  devices  in  each 
leg  of  the  inverter  are  switched  to  give  the  required  output  voltage  at  the  inverter 
terminals.  This  voltage  drives  the  load  to  give  the  desired  voltage  across  the  resistors  and 
capacitors.  Signal  flow  for  the  controller  is  illustrated  in  Figure  4-3.  The  gain  K 
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Figure  4-3  -  Per-Phase  Signal  Flow  of  PI  Controller 


represents  the  gain  inherent  to  the  sine-triangle  PWM  process.  This  gain  is  a  constant  in 
the  linear  region  of  operation  and  decreases  in  a  non-linear  fashion  in  the  overmodulation 
region. 

To  select  gains  for  the  controller  the  transfer  function  for  the  system  must  be 
derived.  Two  methods  can  be  used  to  find  the  transfer  function,  either  drawing  the  block 
diagram  and  using  Mason’s  gain  formula  or  by  using  state  space  methods.  Both  methods 
are  considered. 

First  the  classical  approach  of  drawing  the  block  diagram  is  shown.  This  is  a 
useful  method  in  visualizing  the  signal  flow  of  the  system.  The  flow  diagram  is  shown  in 
Figure  4-4.  It  is  assumed  that  the  converter  will  have  a  unity  gain;  that  is,  the  control 
signal  into  the  converter  times  the  converter  gain  is  the  voltage  which  is  given  at  the 
output.  This  simplification  ignores  nonlinear  effects  and  the  impact  of  harmonics; 
however,  it  is  mathematically  valid  as  long  as  the  converter  is  operated  in  the  linear 
range,  with  a  modulation  index  of  less  than  one. 


t 
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Figure  4-4  -  Signal  Flow  Diagram,  q-Quantities  Only 


In  the  following  discussion  only  the  q-quantities  are  considered  since  the  q-  and  d- 
quantities  are  independent.  Derivation  of  the  d-quantity  controller  mirror  those  shown 
here  for  the  q-quantities. 

Mason’s  rule  states  that  the  input-output  transfer  function  associated  with  the 
signal-flow  graph  of  Figure  4-4  is  given  by,  [Ref.  1 1] 

rt  \  _  Kl  is)  _  Forward  Path  Gain 

V*'L  (5)  1  -  ^  Individual  Loop  Gains  ^ 

Inserting  values  and  solving,  the  gain  or  closed-loop  transfer  function  is  given  by 
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(Eq.  4-15). 


The  characteristic  equation  is  given  by  the  denominator  of  the  closed-loop  transfer 


function 


*3+s2 
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Cf  R/oad 
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(Eq.  4-16) 


An  alternative  method  to  solve  for  the  characteristic  equation  is  to  use  the  state 
equations,  Equations  4-1 1  and  4-12.  Additional  states  are  defined  to  account  for  the 
integrators  employed  in  the  control  as  follows 
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(Eq.  4-17) 


Placing  the  q-quantities  into  state  space  normal  form  results  in  Equation  4-18 
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(Eq.  4-18) 
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where  the  3x3  system  matrix  is  typically  denoted  as  the  system  A  matrix.  To  find  the 
characteristic  equation  the  determinant  of  (si  -  A)  is  found 
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(Eq.  4-19) 


Note  that  Equation  4-19  is  identical  to  Equation  4-16,  as  expected. 

With  the  characteristic  equation  determined,  desired  pole  locations  can  be  selected 
and  implemented  by  choice  of  Kj  and  Kp  and  the  resultant  transient  response 
investigated.  Examining  the  effects  of  the  constants  KjK  and  KpK  on  transient  response 
is  most  easily  accomplished  through  simulation.  A  simple  MATLAB  routine  can  easily 
show  transient  response  and  pole  locations  for  specific  values  of  K\K  and  KpK.  It  is 
typically  desired  to  have  a  quick  response  with  a  limited  amount  of  overshoot.  In 
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addition,  the  gains  should  not  be  made  too  large  or  else  the  PWM  algorithm  is  pushed 
into  overmodulation  and  the  aforementioned  linear  analysis  becomes  invalid. 

To  facilitate  analysis,  the  system  is  transformed  into  the  synchronous  reference 
frame  where  the  desired  input  voltage  is  represented  by  a  steady-state  constant  value.  In 
the  synchronous  reference  frame  the  q-  and  d- variables  are  no  longer  uncoupled,  so  all  the 
states  must  be  considered  together.  The  state  equation  for  the  system  with  a  stationary 
controller,  represented  in  the  synchronous  reference  frame,  is  given  by  Equation  4-20. 
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(Eq.  4-20) 

Steady-state  values  for  load  voltage  in  terms  of  the  desired  voltages  are  given  by 
the  relationship: 

Vload.qi=  AV?oad_gd  (Eq.  4-21) 

where  the  elements  of  A  (not  to  be  confused  with  the  system  matrix)  are  given  by 
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(Eq.  4-22) 


(Eq.  4-23) 


where 
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(Eq.  4-24) 
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(Eq.  4-25) 


Examining  the  diagonal  elements  of  A,  it  can  be  seen  that  if  the  values  of  KpK 
and  KjK  are  made  large  enough  then  the  squared  terms  will  dominate  and  the  value  will 
tend  towards  one.  In  addition  the  off-diagonal  elements  will  tend  towards  zero. 
Reducing  the  steady-state  error  must  be  weighed  against  sending  the  controller  into 
overmodulation,  as  large  gains  will  drive  the  modulation  index  to  be  greater  than  one. 
Overmodulation  is  not  necessarily  an  undesirable  phenomenon,  but  the  previous  analysis 
does  not  hold  and  therefore  the  results  will  not  be  reliable. 

Using  the  values  in  Table  4-1  for  simulation  purposes,  a  settling  time  of  25  msec 
can  be  achieved  with  constant  values  of 


K,K  =  50 
KPK  =  7. 

The  controller  stayed  in  the  linear  modulation  range  for  a  50  Volt  step  change  in 

commanded  voltage  amplitude.  A  larger  step  change  with  the  gains  used  would  send  the 

controller  into  overmodulation.  The  steady-state  error  was  found  to  be  approximately 

1 1%.  If  a  smaller  step  change  were  expected  then  the  gains  could  be  increased  to  further 

reduce  the  steady-state  error. 

Table  4-1  -  Simulation  Parameters 


Lf 

lOmH 

Cf 

O.lmF 

■^ioad 

25  Ohms 
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Results  of  the  MATLAB  program  (Appendix  D,  pi  53)  are  shown  in  Figure  4-5. 
The  waveform  is  the  q-axis  load  voltage  represented  in  the  synchronous  reference  frame. 
The  STEP  function  was  used,  together  with  the  matrices  implicit  to  the  state-space 
description. 


To  verify  that  the  MATLAB  simulation  provides  representative  results,  an  ACSL 
program  (Appendix  D,  ppl54  -  157)  incorporating  a  three-phase  hard-switched  converter 
supplying  a  resistive  load  through  an  LC  filter  was  developed.  A  Pi-controller  was 
implemented  with  the  gains  from  the  MATLAB  simulation.  A  transient  consisting  of  a 
50  volt  step  change  in  the  commanded  q-axis  synchronous  load  voltage  from  the  ACSL 
simulation  (Figure  4-6)  shows  a  settling  time  of  around  25  msec  also.  Examining  the 
dynamics  of  the  step  change  shows  an  almost  identical  response  to  that  of  the  MATLAB 
simulation,  which  shows  that  the  controller  did  not  go  into  overmodulation.  A  significant 
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difference  between  the  ACSL  and  MATLAB  results  would  be  evidence  of 
overmodulation,  since  the  ACSL  program  incorporates  the  dynamics  of  the  sine-triangle 
PWM  scheme  and  the  MATLAB  code  simply  considers  the  PWM  a  constant  linear  gain 
regardless  of  the  input  magnitude. 

ACSL  Step  Response  I 
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Figure  4-6  -  Stationary  Voltage  Controller  Start-up  Transient 

It  is  obvious  that  the  steady-state  error  which  exists  with  the  stationary  controller 
can  be  large.  This  error  is  due  to  the  fact  that  the  control  is  performed  in  the  stationary 
reference  frame  and  the  control  quantities  are  time- varying  sinusoids  in  that  frame. 

One  way  to  correct  for  the  steady-state  error  would  be  to  use  the  A  matrix  from 
Equation  4-21,  which  relates  the  desired  voltages  to  the  steady-state  voltages  in  the 
synchronous  reference  frame.  The  desired  voltages  would  be  transformed  into  the 
synchronous  reference  frame,  and  then  be  pre-multiplied  by  A'1,  to  give  a  new  set  of 
desired  voltages.  This  new  set  of  voltages  would  be  transformed  back  into  the  stationary 
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reference  frame  to  be  used  by  the  controller.  In  this  way  the  steady-state  error  would  be 
eliminated.  The  practicality  of  this  method  would  make  it  difficult  to  implement  as  the 
delays  caused  in  performing  transformations  and  calculating  A'1  would  introduce 
additional  errors.  The  steady-state  error  of  stationary  controllers  will  be  investigated 
further  in  the  next  section. 


B.  CURRENT  FEEDBACK  CONTROLLER 

A  controller  using  current  feedback  together  with  a  commanded  output  current 
waveform  is  considered  in  this  section.  This  type  of  set-up  can  be  used  as  a  first-order 
approximation  of  an  induction  machine  drive,  where  the  induction  machine  is  modeled  as 
a  simple  RL  load.  The  circuit  diagram  is  shown  in  Figure  4-7. 


Figure  4-7  -  Converter  with  LR  Load 


Using  the  same  technique  as  the  voltage  feedback  controller,  the  state  equation  for 
the  current  controller  with  RL  load  is  found  and  listed  as  Equation  4-26.  The  signal  flow 


dt 


Vs  _  p  Is 

qdOL  lvloadlqdO 


(Eq.  4-26) 
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diagram  for  the  current  feedback  controller  is  illustrated  in  Figure  4-8  for  the  q-variables. 
A  similar  diagram  holds  for  the  d-variables.  In  the  stationary  reference  frame  for  the 
system  considered,  the  ‘q’  and  ‘d’  variables  are  not  coupled.  In  Figure  4-8,  K  is  the 
PWM  gain  and  isqs  is  the  desired  q-axis  current. 


Figure  4-8  -  Signal  Flow  Diagram,  q-Quantities  Only 

Using  the  same  techniques  as  before  the  transfer  function  for  the  current  feedback 
controller  is  uncovered  and  is  given  in  Equation  4-27. 


/;  KKs  +  K,K 

77T  =  —  7 - x -  (Eq.  4-27) 

V  L/S2+(RLoad+KpK)s  +  K,K 

Dynamic  response  can  be  investigated  using  simulation  as  in  the  voltage  feedback 
converter.  One  problem  encountered  is  that  the  control  signals  are  not  DC  values,  instead 
they  vary  at  the  output  frequency.  This  also  occurred  with  the  voltage  feedback  converter 
in  the  previous  section.  To  determine  the  steady-state  error,  Equation  4-27  is  used  with 
‘s’  set  equal  to  ja> . 


/;  KpKs  +  K,K 

C  L/s2+(Riml+KpK)s  +  K,K 
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(Eq.  4-28) 


r _ MW _ Tr  (Eq.. 

'  K,K  -  LjO1  +  jo^R^  +  KpK )  " 

It  can  be  demonstrated  from  Equation  4-28  that  unless  the  output  frequency  is  zero,  a 


steady-state  error  will  exist.  For  an  output  frequency  of  60Hz  the  commanded-to-actual 
current  relationship  predicted  by  Equation  4-28  is:  i*s  =  i**s  0.889Z  - 1.36  . 

As  the  previous  analysis  has  pointed  out,  one  drawback  of  using  stationary 
reference  frame  variables  in  the  control  is  that  the  steady-state  error  is  non-zero.  To 
facilitate  further  investigation  of  the  controller,  the  analysis  is  shifted  to  the  synchronous 
reference  frame.  In  this  frame  the  commanded  quantities  will  be  DC  values  and  analysis 
will  be  simpler.  The  synchronous  reference  frame  state  space  formulation  of  the  system 
with  a  stationary  reference  frame  controller  is  given  by  Equation  4-29.  [Ref.  12] 


KpK  +  RLoad 
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-Z—  0 
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J_  1 r*  +  o  ^  (Eq-  4-29> 


L,  LC 


0  0  -e>,  K,K  0 
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The  relationship  between  commanded  and  actual  steady-state  currents  represented  in  the 
synchronous  reference  frame  may  be  derived  from  Equation  4-29  by  setting  the 
derivatives  of  the  state  variables  equal  to  zero.  Synchronous  reference  frame  quantities 


j«  _  A  I** 

lqds  Alqds 


(Eq.  4-30) 


are  used  because  for  a  desired  sinusoidal  three-phase  input,  steady-state  DC  commanded 
quantities  result.  The  elements  of  A  (not  to  be  confused  with  the  system  matrix)  are 


given  as  [Ref.  12] 
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(Eq.  4-31) 
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(Eq.  4-32). 


In  Equations  4-3 1  and  4-32  the  values  for  K'  and  K\  incorporate  the  pulse- width- 

modulation  gain  K.  To  verify  the  accuracy  of  Equations  4-29  through  4-32,  an  ACSL 
simulation  (Appendix  D,  ppl58  - 161)  incorporating  a  controller  in  the  stationary 
Reference  frame  with  an  RL  load  was  developed  and  investigated.  Values  used  in  the 
simulation  are  given  in  Table  4-2. 


Table  4-2  -  Simulation  Parameters 
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10mH 

^•load 

25  Ohms 

Ki’ 

3000 

V 

200 

<ye 

377 

For  the  simulation  the  synchronous  reference  frame  commanded  current  values 
for  the  synchronous  quantities  were  /**  =  5.0  Amps  and  =  0.0  Amps.  Using 
Equations  4-31  and  4-32,  au  is  computed  as  0.889  and  an  as  -0.019.  Thus,  the  steady- 
state  values  for  ie  and  ieds  are  predicted  to  be 

igs  =  flu c  +  anC  =  0.889  *  5.0  -  0.019  *  0.0  =  4.445  Amps. 

&  =  a2Xie;  +  a22C  =  0.019  *  5.0  -  0.889  *  0.0  =  0.095  Amps. 
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This  corresponds  well  with  the  graph  of  commanded  and  actual  currents  during  start-up, 
shown  in  Figure  4-9. 


Figure  4-10  -  Stationary  Current  Controller  Start-up  Error  Transient 
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To  overcome  the  steady-state  error,  the  control  can  be  performed  in  the 
synchronous  reference  frame  where  the  signals  being  controlled  are  all  DC  quantities. 
The  analysis  of  such  a  system  results  in  a  state  space  description  given  by  [Ref.  12] 
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(Eq.  4-33) 


Note  that  the  only  difference  between  Equation  4-29  and  Equation  4-33  is  that  elements 
(3,4)  and  (4,3)  are  now  zero.  The  relationship  between  the  actual  and  commanded  qd- 
currents  is  now  given  by,  [Ref.  12] 

= Ai;i  (Eq- 4-34) 

where  the  elements  of  A  are 


-  a22 

(Eq.  4-35) 

=  <*2l  =  0 

(Eq.  4-36) 

An  ACSL  simulation  (Appendix  D,  pp  162  -  165)  with  a  controller  in  the  synchronous 
reference  frame  was  implemented  to  verify  the  results.  Values  for  the  simulation  are  as 
shown  in  Table  4-2.  Results  of  the  simulation  are  illustrated  in  Figure  4-11.  A  zero 
steady-state  error  was  achieved  as  desired. 
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Figure  4-11  -  Synchronous  Controller  Simulation  Results 


C.  A  VOLTAGE  CONTROLLER  WITH  INNER  CURRENT  CONTROL  LOOP 


A  controller  which  uses  current  and  voltage  feedback  as  well  as  current  and 
voltage  feedforward  has  been  proposed  in  the  literature  [Ref.  13].  The  input  voltage  is 
the  line-to-neutral  voltage  synthesized  from  a  three-phase  inverter  as  described 
previously.  An  LC-filter  is  used  to  reduce  the  voltage  harmonics  across  the  load.  The 
control  is  predicated  on  specifying  a  desired  output  voltage  and  using  various  feedback 
and  feedforward  signals  to  force  the  load  voltage  to  track  that  value.  A  per-phase 
representation  of  the  system  is  illustrated  in  Figure  4-12. 

The  state  equations  for  a  phase  of  the  system  are  given  by 


diL  ^ Load 

_____  — 


(Eq.  4-37) 
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^Load  _  h  hood  _  jc_ 

dt  Cf  Cj 


(Eq.  4-38) 


The  desired  a-phase  of  the  load  voltage  is  specified  as 

v;oad  =  v0  cos(<»./)  (Eq.  4-39) 

where  V0  is  the  desired  peak  amplitude  and  coe  the  desired  frequency  in  radians  per 
second.  Knowing  the  desired  output  voltage  the  desired  capacitor  current  can  be 
determined  from  Equation  4-38. 


i'c  =  C,  =  -Voffl.C,  sinfoj/)  (Eq.  440) 

The  desired  capacitor  current  is  fed  forward  by  adding  it  to  the  actual  load  current  to 
obtain  the  desired  inductor  current.  This  value  is  then  compared  with  the  actual  inductor 
current  to  obtain  a  current  error,  which  is  passed  to  a  PI  controller  to  generate  a  control 
signal  for  a  particular  phase  of  the  converter.  The  current-loop  controller  is  illustrated  in 
Figure  4-13. 

Voltage  feedback  is  incorporated  by  comparing  the  desired  and  actual  values  of 
load  voltage  to  obtain  a  voltage  error.  The  voltage  error  is  input  to  a  PI  controller  to 
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Figure  4-13  -  Per-Phase  Current-Loop  Controller 


generate  a  control  signal.  Feedforward  action  is  realized  by  using  the  desired  inductor 
current  and  desired  load  voltage  to  obtain  the  desired  converter  output  voltage  as 
specified  in  Equation  4-41 . 

+  (Eq.  4-41) 

Equation  4-41  can  be  simplified  by  omitting  the  time  derivative  of  inductor  current  and 
re-writing  the  desired  output  voltage  as: 

V;'=V^+M  il  (Eq.  4-42) 

Because  the  time  derivatives  are  ignored  in  Equation  4-42,  the  feedforward  voltage  is 
only  correct  for  sinusoidal  steady-state  conditions.  Feedback  loops  are  necessary  to 
make-up  for  this  shortcoming.  Feedback  will  also  correct  for  discrepancies  between  the 
value  of  inductance  used  in  the  controller  and  the  actual  value  of  the  circuit  filter 
inductance  which  may  change  due  to  temperature  and  saturation  effects.  Figure  4-14 
shows  the  signal  flow  for  the  per-phase  controller. 

Implementation  of  the  controller  is  accomplished  in  the  synchronous  reference 
frame  where  sinusoidal  quantities  become  DC  quantities.  This  simplifies  the  control  and 
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Figure  4-14  -  Per-Phase  Controller 


allows  the  integral  action  to  guarantee  zero  steady-state  error.  The  synchronous  reference 
frame  transformation  is  given  by  Equation  4-43  [Ref.  10]. 


2 

k:=j 


2k 


cos (#e)  cos  6e - cos 


3 

2  n 


0.+ 


2k 


sin(0e)  sin  0e - sin 


1 

2 


1_ 

2 


e’+f 

2 


(Eq.  4-43) 


where 

/ 

0e  =  \coedt  (Eq.  4-44) 

o 

The  ‘e’  superscript  denotes  the  synchronous  reference  frame.  The  transformation  of  a 
balanced  set  of  three-phase  voltages  is  illustrated  in  Equation  4-45. 

Zero-sequence  quantities  are  zero  as  expected,  and  note  that  for  the  particular  set 
of  voltages  selected,  the  commanded  d-voltage  is  also  zero.  Transformation  of  control 
quantities  into  the  synchronous  frame  is  straightforward  and  will  not  be  discussed  in  any 
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V0  cos(<ye  t) 


=  K *  V0cos ycQet-  — 

„  (  2  n 

V0cos^©,  t  +  — 


(Eq.  4-45) 


additional  detail.  The  composite  controller  in  the  synchronous  reference  frame  is  shown 
in  Figure  4-15. 
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Figure  4-15  -  Signal  Flow  for  Synchronous  Controller 
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To  choose  gains  for  the  PI  controllers  the  transfer  function  relating  the 
commanded  and  actual  output  voltages  must  be  determined.  Assuming  an  RL  load,  the 
system  state  equations  in  the  synchronous  reference  frame  are  given  by 


y  £  T  T  g  m  g  f  ^ 

' Load-q  Load-d  1  *  L-q  1  ^  Load-q 

Ve  to  Ve  Cf  7e  Cf  ,•« 

■Load-d  J  u/er  Load-q  J  J  \J  L-d  J  J  1  Load-d  _ 


(Eq.  4-46) 
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(Eq.  4-48) 


[}  Load-d]  L  Load-q  J  ~Load  ^  Load_d  ]  ~ Load  ^  Load_d  ] 

In  order  to  incorporate  the  dynamics  of  the  Pi-compensators,  the  following  additional 
state  equations  are  required 
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(Eq.  4-49) 
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(Eq.  4-50) 


Equations  4-5 1  and  4-52  show  the  synchronous  reference  frame  commanded  voltages 
represented  in  terms  of  state  variables, 


v;:q  =(KpcKpv -a]LfCf)w0+KpcKivXV;  +  a>eLfKivXV*d  +  Kicxi'q 

+  ^ pc  ^-pv)^ Load-q  ~  e^f  ^pv)^ Load-d  +  ^pc^  Load-q  +  ^e^/hoad-d  ~  & pch-i 

Vf-d  =  -(V»eCf  +coeLfKpv)wo -C0eLfKlvxv;+KpcKivxv‘+Kicxi‘ 

+  ooeLfKpvV [oad_q  +(l-KpcKpv)V[oad_d -coeLfilad_q+KpJload_d-KpJlc 


(Eq.  4-51) 


(Eq.  4-52) 
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where  V0  is  the  commanded  q-axis  voltage  V** .  The  actual  three-phase  source  voltages 
from  the  converter  are  the  control  voltages  times  the  PWM  gain,  in  the  synchronous 
reference  frame  this  is 


v;.,  -  W (Eq.  4-53) 
V’,.a  =  Kpmv;_a  (Eq.  4-54) 

Combining  Equations  4-46  through  4-54,  the  state  space  representation  of  the  system  is 
generated  and  shown  in  Equation  4-55. 

From  the  state  space  representation  of  Equation  4-55,  it  is  desired  to  obtain  a 
transfer  function  relating  the  commanded  voltage  to  the  actual  load  voltage.  With  the 
transfer  function  known,  controller  gains  can  be  selected  to  give  the  system  a  desired 
stable  response.  As  the  transfer  function  will  be  of  tenth-order,  an  analytical  solution  was 
not  attempted.  Instead,  a  MATLAB  command  file  was  used  to  obtain  transfer  functions 
for  various  gains.  The  poles  of  the  system  could  not  all  be  placed  by  choosing  gains, 
since  there  are  only  four  gains  to  choose  and  ten  system  poles.  Actual  values  examined 
in  the  analysis  are  documented  in  Table  4-3. 


Table  4-3  -  Simulation  Parameters 


Lf 

10.1  mH 

Cf 

1.0  mF 

©e 

377  rad/sec 

^“Load 

25  ohms 

^Load 

1.0  mH 

Kpv 

0.25 

Kpc 

10.0 

Kiv 

1.0 

Kic 

1.6 
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IIMirf  ..  [  Ad  3  tUA\d 


The  resultant  poles  and  zeros  of  the  transfer  function 


Load-q 


are  presented  in  Table  4-4. 


Table  4-4  -  Poles  and  Zeros  of  Transfer  Function 


Poles 

Zeros 

-249600+j376 

-419+jl815 

-969+jl060 

-3.92+j0.59 

-0.087+j0.033 

-249800+j433 
-907+j  1 500 
-3.92jj0.59 
-0.087+0.033 

The  response  of  the  system  to  a  50  volt  step  change  in  commanded  voltage,  V0, 
was  simulated  in  MATLAB  using  the  STEP  function  and  the  state  space  description  of 
the  system  given  in  Equation  4-55.  Figure  4-16  illustrates  the  MATLAB  output  of  the  50 
volt  step  input.  The  MATLAB  program  is  shown  in  Appendix  D,  pp  166  -  167. 
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When  considering  the  poles  of  the  transfer  function  shown  in  Table  4-4,  the 
MATLAB  step  response  does  not  seem  to  give  the  correct  response.  It  would  appear  that 
the  dominant  poles  would  be  the  set  closest  to  the  origin  (-0.087+j0.033).  If  this  were  the 
case  then  the  settling  time  for  the  step  response  would  be  much  longer  than  that  seen  by 
the  simulation.  To  explain  the  rapid  settling  time  the  zeros  of  the  transfer  function  must 
also  be  examined.  Taking  the  zeros  into  account,  pole-zero  cancellation  occurs  leaving 
the  dominant  set  of  poles  as  -419+j  1815.  The  settling  time  from  Figure  4-17  is 

4.6 

approximately  0.01  seconds.  Using  the  approximation  ts  =  —  ,  where  o  is  the  negative 

G 

46 

real  part  of  the  dominant  pole  [Ref.  1 1],  the  settling  time  would  be  ts  =  =  0.01 1  sec . 

Thus  the  poles  located  at  -419+j  1815  are  the  dominant  poles  for  the  system. 

To  verify  that  the  MATLAB  simulation  was  valid,  a  detailed  simulation  of  the 
system  and  controller  using  ACSL  was  performed.  The  simulation  (Appendix  D,  pp  168 
-171)  incorporated  the  parameter  values  listed  in  Table  4-4.  A  hard-switched  converter 
was  used  in  the  simulation  with  a  sine-triangle  control  scheme.  The  three  modulating 
sinusoids  for  the  control  scheme  are  the  Va,  Vfr,  and  Vc  signals  illustrated  in  Figure  4-15. 
A  50  volt  step  change  in  commanded  voltage  was  simulated  and  the  results  are  shown  in 
Figure  4-17. 

Comparing  Figures  4-16  and  4-17  shows  an  almost  identical  response  for  the  step 
change.  This  validates  the  MATLAB  simulation  as  an  accurate  model  of  the  system  and 
controller,  at  least  in  the  linear  operating  range  of  the  PWM  scheme. 
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While  investigating  the  effect  of  controller  gains  on  system  response  with  the 
MATLAB  model,  it  was  noted  that  larger  gains  tended  to  speed  up  the  response.  In  fact 
system  responses  in  the  order  of  a  few  microseconds  were  achieved  using  very  large 
gains.  When  these  same  gains  were  used  in  the  ACSL  simulation  the  results  showed  very 
erratic  and  unpredictable  outputs.  This  variation  is  due  to  the  fact  that  the  state  space 
description  used  in  the  MATLAB  code  assumed  that  no  matter  what  the  commanded 
signal  specified  by  the  controller,  the  inverter  could  theoretically  produce  it.  Thus,  if  the 
commanded  signal  was  50000  volts,  that  is  what  the  system  would  produce.  In  reality  the 
system  cannot  put  out  a  voltage  above  the  source  voltage  (which  was  500  V  for  the 
simulations).  Thus  if  the  control  signal  became  too  large  saturation  occurred  which  was 
not  represented  by  the  MATLAB  model. 


« 
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This  problem  could  be  overcome  by  using  a  source  with  a  very  large  value,  but 
this  is  impractical  in  most  realizable  systems.  Instead  the  gains  were  chosen  lower  so  that 
normal  transients  would  not  drive  the  controller  into  saturation.  This  gave  a  predictable 
response  to  transients. 

Now  that  the  dominant  set  of  poles  has  been  uncovered  for  the  set  of  gains 
chosen,  the  effect  of  different  gains  must  be  addressed.  Specifically,  does  changing  the 
gains  negate  the  pole-zero  cancellation  seen  earlier  and  thus  shift  the  dominant  poles 
closer  to  the  origin.  To  investigate  this  several  sets  of  gains  were  simulated  and  the  poles 
and  zeros  created  were  compared.  Results  for  two  different  sets  of  poles  are  listed  in 
Tables  4-5  and  4-6. 


Table  4-5  -  Poles  and  Zeros 
Gains:  Kpv=1.5  Kpc=15  Kiv=10  Kic=5 


Poles 

Zeros 

-249600+j376 

-350+j5043 

-1531+j4288 

-6.66+j0.167 

-0.313+j0.0079 

-249800+j380 
-1046+j4775 
-6. 66+j0.169 
-0.313±0.0079 

Table  4-6  -  Poles  and  Zeors 
Gains:  Kpv=0.1  Kpc=1.75  Kiv=0.01  Kic=0.01 


Poles 

Zeros 

-249600+j376 

-10.54+j78.66 

-563+J32.0 

-0.0876+j0.0330 

-.000997+j0.00217 

-249800+jl530 

-10700 

-590.8 

-0.0876+j0.0330 
-.000997+j 0.002 1 7 
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Inspection  of  the  results  shows  a  pole-zero  cancellation  for  the  four  poles  closest 
to  zero  for  both  sets  of  gains.  All  of  the  other  sets  of  gains  investigated  showed  the  same 
cancellation.  It  is  therefore  not  a  particular  set  of  gains  which  cause  the  four  poles  closest 
to  zero  to  be  canceled,  but  the  system  configuration  which  causes  it. 

With  the  knowledge  of  pole-zero  cancellation,  the  system  response  can  be  better 
optimized.  Again  using  simulation  as  a  tool,  several  sets  of  gains  were  programmed  and 
the  responses  to  a  step  input  analyzed.  An  optimum  set  of  gains  was  chosen  and  is 
illustrated  in  Table  4-7  along  with  the  poles  and  zeros  of  the  transfer  function. 


Table  4-7  -  Optimum  Poles  and  Zeros 
Gains  -  Kpv=0.05  Kpc=10  Kiv-0.05  Kic=0.05 


Poles 

Zeros 

-249600+j376 
-814.8+j220.1 
-576.8+j534.8 
-0.632+J0.482 
-.0043 7+j0.00 166 

-249800+jl530 

-1955 

-638.4 

-0.632+j0.482 

-,00437+j0.00166 

The  MATLAB  generated  system  response  to  a  50  volt  step  increase  is  depicted  in 
Figure  4-18. 

To  confirm  that  the  controller  did  not  saturate  and  that  the  MATLAB  code  gave 
realistic  results,  an  ACSL  simulation  was  performed  and  the  results  reported  in  Figure  4- 
19. 

Comparing  Figures  4-18  and  4-19  shows  an  almost  identical  response  with  a 
settling  time  of  about  0.009  seconds.  With  dominant  poles  of  -576.8+j534.8  the 

46  46 

predicted  settling  time  is  —  =  — : —  =  0.008  seconds,  which  corresponds  well  with  the 

<j  576.8 

step  response  plots. 
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Volts  |  Amplitude 


Figure  4-18  -  MATLAB  Step  Response  with  Optimized  Gains 
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If  a  faster  response  were  desired  a  bang-bang  type  of  control  could  be  derived 
where  the  system  would  be  driven  by  either  the  maximum  or  minimum  control  effort 
until  the  transient  was  stabilized. 

Normally  large  step  changes  in  commanded  values  are  not  encountered.  System 
start-up  most  likely  would  be  accomplished  with  manual  control  or  by  using  a  ramp  input 
in  commanded  voltage.  This  minimizes  the  importance  of  a  quick  response  to  step 
changes  and  emphasizes  the  steady-state  control  for  which  the  controller  was  designed. 

With  the  control  issue  addressed,  investigation  of  a  new  modulation  technique  is 
undertaken  with  a  focus  on  improving  the  classical  sine-triangle  modulation  method. 
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V.  PRINCIPLES  OF  SPACE  VECTOR  CONTROL 


A  new  method  of  pulse-width-modulation  has  been  introduced  recently  in  the 
literature  [Ref.  14,15].  This  chapter  explains  the  basics  of  this  new  technique  called 
Space  Vector  Control  (SVC).  In  addition,  a  simulation  showing  the  feasibility  of  the 
concept  is  performed,  and  finally  the  advantages  and  disadvantages  are  discussed. 


A.  DEVELOPMENT  OF  SPACE  VECTOR  CONTROL  ALGORITHM 


The  basic  ideas  of  Space  Vector  Control  have  been  introduced  and  discussed  in 
the  literature  [Ref.  16,17].  The  following  discussion  lays  out  the  steps  necessary  to 
implement  a  SVC-based  algorithm.  To  realize  Space  Vector  Control  for  three-phase 
applications,  a  set  of  three-phase  variables  is  represented  by  a  single  vector.  This  is 
accomplished  by  transforming  quantities  in  the  a-b-c  reference  frame  into  the  stationary 
q-d-0  frame.  The  three  phase  voltages  are  denoted  by  the  notation 


[Vcm  V bn  Kn]r  =  v abcn  •  A  transformation  matrix  is  then  applied  to  the  voltage  vector 
placing  it  into  the  desired  reference  frame  [Ref.  10]. 


^  qdOn  ^  abcn 

"  (2^ 

I  (  Ik'] 

cos#  cos  # - 

v  3 ; 

1  V  3) 

2 

.  .  (  In' 

(  2  k' 

K  =  — 

sin#  sm  #-  — 

sin  #  +  — 

,  3 

‘V  3  J 

V  3) 

1  1  1 


(Eq.  5-1) 


(Eq.  5-2) 


2  2  2 


v 


qdOn 


V, * 

/On 


(Eq.  5-3) 
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(Eq.  5-4) 


6  -  jcodt  +  00 

o 

For  the  stationary  reference  frame,  co  is  zero  and  6 Q  is  selected  to  also  equal  zero  for 
convenience. 

Assuming  that  the  load  is  balanced  the  zero  sequence  voltage,  V0n,  will  be  zero. 
The  voltages  Vqn  and  can  be  plotted  along  the  stationary  q-  and  d-axes,  respectively. 


Figure  5-1  -  q-  and  d-axis  Voltages 


Figure  5-2  illustrates  the  vector  addition  of  the  two  instantaneous  voltages  directed  along 
the  q-  and  d-axes.  The  resultant  V*  provides  a  vector  representation  for  the  three  phase 
voltages,  as  desired. 


Figure  5-2  -  Single  Vector  Representation 
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Since  the  phase  voltages  are  sinusoids  and  thus  time  varying,  the  transformation 
into  the  stationary  reference  ffame(u>  =  0  and  0O  =  0)  will  produce  a  time-varying 
resultant  vector  as  well.  If  the  q-axis  is  viewed  as  the  imaginary  axis  and  the  d-axis  the 
real  axis,  this  resultant  voltage  can  be  thought  of  as  a  phasor  representation  of  the  phase 
voltages.  By  specifying  the  magnitude  and  phase  angle  of  a  reference  phasor,  any  desired 
set  of  balanced  three-phase  voltages  can  be  generated  and  represented  conveniently.  The 
space  vector  method  of  control  uses  this  reference  vector  in  conjunction  with  breaking  the 
plane  into  six  sectors  to  establish  the  rationale  for  switching  the  three-phase  inverter  used 
to  synthesize  the  voltages. 

In  a  six-switch  inverter  there  are  eight  possible  switch  combinations  as  shown  in 
Figure  5-3.  Note  that  electrically  State  7  and  State  8  are  identical,  both  producing  zero 
volts  across  each  phase. 


Figure  5-3  -  Eight  Possible  Switch  States  for  Six-Switch  Inverter 
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If  it  is  assumed  that  the  inverter  is  connected  to  a  balanced  three-phase  wye- 
connected  load,  the  line-to-neutral  voltages  corresponding  to  each  of  the  above  states  are 
summarized  in  Table  5-1. 

Table  5-1  -  Phase  Voltages 


State  # 

Van 

Vbn 

Vcn 

1 

-V^ 

2 

1 

Vdc/3 

-2  Vdc/3 

3 

■  Vdc/3 

2Vdc/3 

jaaWai 

4 

-2Vdc/3 

Vdc/3 

5 

-  Vdc/3 

-  Vdc/3 

2  Vdc/3 

6 

Vdc/3 

-2  Vdc/3 

Vdc/3 

7 

0 

0 

0 

8 

0 

0 

0 

Transforming  the  phase  voltages  corresponding  to  each  state  into  the  stationary 


q-d-0  frame  results  in  the  mapping  illustrated  in  Figure  5-4.  Note  that  the  magnitude  of 
.  2 

each  of  the  vectors  is  -v  . 

3  00 


Figure  5-4  -  Inverter  Voltages  in  q-d  Frame 
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These  states  divide  the  plane  into  six  sectors  each  spanning  60  degrees,  as  shown 


in  Figure  5-5. 


Next,  a  ‘mean’  space  vector,  V ,  is  defined  as  the  weighted  average  of  the  line-to- 
neutral  voltages  used  during  a  switching  period,  expressed  as  a  single  vector  in  the 
stationary  q-d-0  frame.  A  switching  period  covers  the  time  period  it  takes  to  generate  the 
average  voltage  vector.  If  four  separate  switch  states  are  used  to  create  V ,  then  the 
period  would  cover  the  time  required  to  switch  to  all  four  of  these  states.  In  SVC 

normally  only  three  separate  line-to-neutral  voltage  combinations  are  used  to  generate  V . 

These  line-to  neutral  voltages  are  controlled  so  that  V  tracks  the  desired  voltage  vector 
V*.  The  error  for  one  switching  period  is  the  mean  voltage  vector  subtracted  from  the 
desired  voltage  vector  and  integrated  over  one  switching  period. 
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'/ 

error  -  J(v*-V)df  (Eq.  5-5) 

'o 

where, 

V*  =  desired  Voltage  Vector 
V  =  actual  Voltage  Vector 
t0  =  beginning  of  switching  period 
tf  =  end  of  switching  period 

This  error  can  be  maintained  small  if  the  cycle  time  is  as  short  as  possible  and  the 
ratio  of  switching  frequency  to  output  frequency  is  large.  Both  of  the  previous 
techniques,  maintaining  a  short  cycle  time  and  a  large  switching  to  output  frequency 
ratio,  operate  to  keep  V*  approximately  constant  over  a  switching  period.  Since  V*  is 
rotating  in  the  stationary  q-d-0  plane  at  the  electrical  output  frequency,  e.g.  60Hz,  if  the 
switching  frequency  is  sufficiently  large,  somewhere  in  the  kilohertz  range,  then  V*  can 
be  assumed  constant  over  one  switching  period.  Keeping  V*  constant  during  a  switching 
period  is  desirable  because  the  generation  of  V  will  not  be  attempting  to  hit  a  “moving 
target.” 

To  keep  the  cycle  time  short,  the  three  switching  states  adjacent  to  the  reference 
vector  are  only  used,  and  these  states  are  used  only  once  each  to  generate  an  average 
voltage  vector  equal  to  the  reference  vector.  For  example,  if  V*  lies  in  Sector  I  then  only 

the  states  adjacent  to  Sector  I  (States  1,2, 7,8)  will  be  used  to  generate  V .  In  addition  the 
switching  frequency  is  made  as  large  as  practical  so  that  the  switching  period,  Tz  =  — - — 

f switch 
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the  sum  of  the  times  at  States  1,2,7  and  8,  will  be  small.  If  the  switching  frequency  were 


selected  to  be  1 0kHz  then  the  cycle  time  would  be  T  = - =  0.000 1  sec . 

z  10000 Hz 

The  switching  sequences  are  arranged  in  such  a  way  that  a  transition  from  one 
state  to  the  next  is  performed  by  switching  only  one  inverter  leg;  this  minimizes  the 
switching  frequency  for  each  inverter  leg.  It  can  be  seen  in  the  switching  diagram  that 
these  conditions  are  met  if  the  inverter  legs  are  switched  in  a  sequence  starting  at  one  zero 
state  and  ending  at  the  other  zero  state.  Thus,  both  zero  states  must  be  used.  For 
example,  if  the  reference  vector  lies  in  Sector  I  then  the  switching  sequence  would  be 
. .  .8, 1,2, 7, 2, 1,8,1, 2,7,2  . . .  Starting  at  State  8,  phase  A  is  switched  high  to  arrive  at  State 
1 .  From  State  1 ,  phase  C  is  switched  high  to  arrive  at  State  2.  To  arrive  at  State  7  phase 
C  is  switched  high.  Note  that  the  switching  sequence  is  8, 1,2, 7...  and  not  8, 2, 1,7... .  This 
is  to  minimize  the  number  of  switchings  required.  To  go  from  State  8  to  State  2  both 
phase  A  and  phase  B  must  be  switched;  however,  to  go  from  State  8  to  State  1  only  phase 
A  must  be  switched.  The  times  spent  at  each  state  are  variables  to  be  calculated  and  are 
used  to  force  the  average  voltage  equal  to  the  desired  voltage.  This  process  repeats  until 
V*  rotates  into  a  new  sector.  The  transition  to  the  next  sector  is  accomplished  when  the 
sequence  is  at  a  zero  state,  since  the  zero  states  are  used  in  every  sector.  Table  5-2  lists 
the  switching  sequence  to  be  used  for  each  sector. 

The  average  voltage  generated  for  one  cycle  of  switching  from  states  1-2-7  is  then 


(Eq.  5-6) 
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Table  5-2  -  Switching  Sequence 


Sector 

Sequence 

I  ! 

...8, 1,2, 7, 2,1, 8... 

II 

...8, 3,2, 7, 2, 3, 8... 

III 

...8, 3,4, 7, 4, 3, 8... 

IV 

...8, 5,4, 7, 4, 5, 8... 

V 

...8, 5,6, 7, 6,5, 8... 

VI 

...8, 1,6, 7, 6, 1,8... 

Upon  determining  tslalel ,  half  of  the  time  will  used  for  State  7  and  half  of  the  time  will  be 
used  for  State  8.  Since  these  two  states  are  electrically  equivalent,  the  average  voltage 


will  not  be  affected.  This  average  voltage  value  is  equal  to  V ,  which  is  approximately 
the  desired  voltage  vector,  V*. 


Defining  variables  as  :  Tz  = 


1 


f , 'switch 

V  =  V 

y  1  Y  state  \ 


V  =  V 

T  2  state  2 


^1  ^  state  1 


^2  ^  state  2 


^ state  1  ^  state!  ^ state  1 


it  is  clear  from  Equation  5-6  that 


VA+V2t2=VT; 


(Eq.  5-7) 


Describing  the  space  vectors  of  Figure  5-6  in  rectangular  coordinates, 
equation  becomes: 


,the  above 


t  —  V 

*1  DC 


+  h \v°c 


co: 

sin) 


s(6o°y 

(60-)_ 


=  T, -rx  a\ 


cos(y) 

sin(y) 


(Eq.  5-8) 


where 


y  =  cot  and  a  = 


|V| 


DC 
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Figure  5-6  -  Determination  of  Switching  Times 


Equation  5-8  is  solved  for  t{  and  t2,  as  shown  in  Equations  5-9  and  5-10.  Knowing  t,  and 
t2,  the  zero  state  time,  t7,  can  be  easily  found  from  Equation  5-11. 


sm 


Ta  v  ’ 

1  sin(60°) 

(Eq.  5-9) 

_r,  sinM 

1  sin(60°) 

(Eq.  5-10) 

=  71  -  (r,  +  r2 ) 

(Eq.  5-11) 

Since  the  voltage  is  zero  at  both  State  7  and  State  8,  it  is  found  to  be  beneficial  in 
terms  of  symmetry  to  allocate  half  of  the  zero  voltage  time  to  State  7  and  half  to  State  8. 


By  redefining  t7  as  T0  (time  at  zero  state)  and  dividing  by  two,  the  desired  zero  state  times 

T 

are  generated.  Thus,  states  seven  and  eight  are  both  at  zero  for  -2-  seconds.  This  leads  to 

2 

the  pulse  pattern  shown  in  Figure  5-7.  Switching  in  other  sectors  follows  the  same 
general  rules  as  in  Sector  one. 

A  situation  may  arise  where  the  reference  vector  is  so  large  that  it  cannot  be 
generated  in  one  cycle  i.e.  V*  >  V,  +  V2  (Figure  5-6).  This  problem  can  be  handled  by 
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Figure  5-7  -  Voltage  Pattern,  Sector  One 


limiting  the  control  vector  to  2/3  of  the  DC  line  voltage,  if  a  value  larger  than  this  is 
required  the  q-  and  d-components  of  the  control  must  be  scaled  to  achieve  a  smaller 
magnitude  while  maintaining  the  same  relative  proportions.  This  is  accomplished  by 

|V|  o 

limiting  the  value  of  a  = - to  -~VDC  .  Additionally,  the  time  to  switch  from  a  zero 

V/x  3 

state  to  an  intermediate  state  must  be  considered.  A  minimum  switching  time  exists,  this 
may  vary  depending  on  the  type  of  inverter  in  use.  Thus  the  largest  voltage  which  can  be 
generated  in  one  cycle  becomes 

|v|  =-Vnr  ~27^.  (Eq.  5-12) 

I'  Max  3  00  Tz 

B.  SIMULATION 

Simulation  of  a  six-switch  inverter  controlled  by  space  vector  modulation  was 
accomplished  using  Advanced  Continuous  Simulation  Language  (ACSL).  An  initial 
open-loop  simulation  without  voltage  or  current  feedback  shows  the  feasibility  of  the 
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control  scheme.  The  circuit  used  in  the  simulation  is  shown  in  Figure  5-8.  The  program 
(listed  in  Appendix  E)  implements  the  LC  output  filter  and  the  resistive  load  shown  in  the 
circuit  diagram.  Values  used  in  the  simulation  are  reported  in  Table  5-3,  and  the 
flowchart  for  the  simulation  is  shown  in  Figure  5-9. 

Table  5-3  -  Simulation  Parameters 


Vdc 

500V 

Lf 

lOmH 

cf 

lmF 

■^load 

25  Ohms 

f 

^Switch 

6kHz 

v,.n 

140V 

fou. 

60Hz 

Figure  5-8  -  Circuit  Diagram  for  Simulation 


The  initial  section  of  the  ACSL  code  sets  up  the  integration  algorithm  and 
constants  for  use  by  the  rest  of  the  program.  The  main  body  of  the  program,  the 
derivative  section,  begins  by  generating  a  desired  three-phase  voltage  waveform  and 
converts  it  into  the  stationary  q-d-0  frame.  Next  the  sector  of  operation  is  determined 
from  the  q-d-0  voltage  values.  Switching  times  are  calculated  using  angles  from  the 
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Figure  5-9  -  Simulation  Flowchart 


q-d-0  values,  a  maximum  limit  is  placed  on  the  switching  times  by  using  an  if  statement 
to  limit  the  magnitude  of  V  (Vbar)  to  VDC.  A  minimum  switching  time  is  incorporated 
by  the  use  of  a  bound  statement  in  the  generation  of  the  switching  times.  Schedule 
statements  perform  the  switching  by  calling  appropriate  subroutines  located  in  discrete 
blocks  at  the  end  of  the  derivative  section.  These  subroutines  generate  the  appropriate 
output  voltages  for  the  switch  configuration.  The  remainder  of  the  derivative  section 
transforms  the  output  voltages  into  q-d-0  quantities,  calculates  state  variables,  and  then 
transforms  quantities  back  into  the  a-b-c  reference  frame. 

Figure  5-10  shows  the  output  voltage  across  a  single  phase  of  the  load  resistance 
and  the  current  output  from  a  single  phase  of  the  converter.  As  expected  the  current  and 
voltage  are  out  of  phase  and  are  both  sinusoidal. 
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Figure  5-10  -  Output  Voltage  and  Current  -  Phase  A 

The  voltage  pattern  for  two  complete  switching  cycles  is  shown  in  Figure  5-11. 
As  expected  each  phase  changes  polarity  only  once  per  cycle. 


Figure  5-11  -  Voltage  Pattern  -  Sector  Six 


C.  ADVANTAGES  OF  SPACE  VECTOR  CONTROL 

With  the  advent  of  the  microprocessor  Space  Vector  Control  (SVC)  has  become 


very  practical  control  method  to  implement.  A  simple  look-up  table  based  on  the 


magnitude,  phase  angle  and  sector  location  of  a  reference  phasor  can  be  used  by  a 
microprocessor  to  generate  an  SVC  switching  sequence.  Traditional  sine-triangle 
modulation,  where  a  sinewave  is  superimposed  on  a  high-frequency  triangle  wave,  is 
much  simpler  to  implement  using  analog  methods  than  SVC.  As  most  modem 
controllers  are  microprocessor  based,  in  which  case  both  sine-triangle  and  SVC  are 
implemented  using  look-up  tables,  this  analog  disadvantage  of  SVC  is  lost. 

To  compare  the  maximum  output  amplitudes  possible  with  sine-triangle  and  SVC 
the  modulation  indices  must  be  studied.  The  modulation  index  for  sine-triangle 
modulation,  mST,  is  defined  as  the  fraction  of  possible  maximum  pulse  width.  The  line- 
to-neutral  voltage  for  sine-triangle  modulation  in  terms  of  the  modulation  index  is  given 


as: 

(Eq.5-13) 

For  SVC  the  modulation  index,  msvc,  can  be  obtained  by  referring  to  Figure  5-6  as: 


V,.> 


(Eq.  5-14) 


In  sine-triangle  modulation  when  the  modulation  index  becomes  larger  than  one  the 


overmodulation  region  is  entered.  This  is  normally  undesirable  since  the  output 


waveform  becomes  more  distorted  and  lower  order  harmonics  are  generated. 

Comparing  Equations  5-13  and  5-14,  a  relationship  between  mST  and  msyc  can  be 

uncovered. 


4 

mST  —  2  mSVC 


(Eq.  5-15) 
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Further  investigation  shows  that  the  largest  line-to-line  voltage  possible  without 
distortion  for  SVC  is  the  DC  source  voltage  V^.  Using  Equation  5-14  this  gives  the 
maximum  modulation  index  for  nondistorted  waveform  from  SVC  as: 


^-/max=^=V3F/_„=V3/WiVC|Fi 


DC 


(Eq.  5-16) 


Relating  this  back  to  sine-triangle  modulation  gives  a  modulation  index  of 


4  4  -s/3 

mST  —  ^  mSVC  —  2  2 


-7=  =  1.155. 

s 


(Eq.  5-17) 


Comparing  Equations  5-16  and  5-17  it  can  be  seen  that  the  peak  amplitude  non-distorted 
waveform  which  can  be  generated  by  SVC  is  approximately  15%  higher  than  the 
maximum  sine-triangle  generated  waveform. 

When  comparing  the  switching  frequency  of  individual  devices  using  SVC  with 
that  of  sine-triangle  modulation,  it  is  noted  that  for  one  cycle  of  SVC  each  of  the  three 
legs  of  a  six-switch  converter  change  state  once.  So  each  switch  is  operated  once. 
However  for  sine-triangle  modulation  each  leg  of  the  converter  changes  state  twice 
during  one  cycle,  so  each  switch  is  operated  twice.  Therefore  when  comparing 
waveforms  generated  by  each  of  the  modulation  techniques  this  relationship.  Equation  5- 
18,  must  be  taken  into  account. 


f  =  —  f 

J  switch  sine-triangle  ^  J  switch  SVC 


(Eq.  5-18) 
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The  previous  chapters  have  dealt  entirely  with  circuit  description  and  system 
simulation.  It  is  desired  to  test  the  previous  findings  on  an  actual  circuit  model.  This  will 
validate  the  findings  and  give  greater  confidence  to  the  simulation  results.  The  next 
chapter  focuses  on  comparing  simulated  and  actual  circuit  performance  under  various 
conditions. 
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VI.  IMPLEMENTATION  AND  TESTING  OF  CONTROL  ALGORITHMS 


In  order  to  test  the  control  algorithms  a  hard-switched  inverter  was  constructed 
using  lab  components.  Control  algorithms  were  realized  with  a  real-time  interface  system 
consisting  of  a  486  computer  and  a  Digital  Signal  Processor  (DSP).  This  system  allowed 
the  control  algorithm  to  be  programmed  in  SIMULINK,  automatically  translated  into  C- 
code  and  then  downloaded  into  the  DSP.  This  chapter  includes  a  brief  description  of  the 
system  used  for  testing,  a  discussion  of  the  process  of  implementing  control  schemes,  a 
presentation  of  test  results,  and  a  correlation  of  the  test  results  with  theory. 

A.  d-SPACE  SYSTEM  DESCRIPTION 

The  dSPACE  (digital  signal  processing  and  control  engineering)  system  is  a 
technology  which  allows  for  real-time  processing,  analysis  and  control  of  digital  signals. 
The  specific  hardware  used  is  the  dSPACE-DSl  102  controller  board,  which  is  built 
around  a  TMS320C31  floating  point  DSP.  The  DS1 102  contains  two  DSP’s.  The  first  is 
the  main  processor,  a  TMS320C3 1 ,  which  executes  the  main  user  application.  The 
second  is  the  TMS320P14  slave-DSP  used  for  many  different  digital  I/O  functions.  [18] 
SIMULINK  is  normally  used  to  generate  the  main  user  application.  However,  for  the 
DSP  to  be  able  to  read  the  program  it  must  be  in  C  code.  To  solve  this  problem  the 
system  incorporates  an  automatic  C  code  generator  which  converts  the  MATLAB  code 
generated  from  SIMULINK  into  C  and  downloads  it  into  the  main  processor.  Figure  6-1 
illustrates  a  flowchart  for  the  dSPACE  system. 
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Figure  6-1  -  dSPACE  System  Flowchart 


B.  SYSTEM  USED  FOR  TESTING 

A  hard-switched  inverter  was  constructed  using  a  Transistor  Voltage  Source 
Inverter  manufactured  by  INVERPOWER  Controls  Ltd.  It  consists  of  six 
transistor/diode  pairs  arranged  for  use  in  three-phase  applications  as  illustrated  in  Figure 
6-2.  Darlington  power  transistors  are  used  for  the  switches.  Gating  signals  for  the 
transistors  are  supplied  by  the  output  of  the  d-SPACE  board,  which  is  passed  through  a 
pulse  amplifier  module  (also  built  by  INVERPOWER)  as  the  d-SPACE  output  is  not 
appropriate  to  directly  drive  the  six  transistors. 

Source  power  for  the  circuit  was  provided  by  a  three-phase  rectifier  connected  to  a 
large  capacitor  bank  (a  total  of  30mF).  This  provided  a  stable  DC  voltage.  A  varic  was 
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connected  between  the  three-phase  source  and  the  rectifier.  This  arrangement  allowed  the 
DC  voltage  to  be  varied  form  zero  to  220  volts. 

A  three-phase  RL-load  was  built  with  standard  laboratory  equipment.  The  resistor 
banks  used  are  also  manufactured  by  INVERPOWER  and  consist  of  six  116  Q/500  Watt 
resistors  which  can  be  easily  arranged  in  series  or  parallel  as  required.  The  inductors 
used  are  43.5  mH/10  Amp  reactors,  built  by  INVERPOWER,  and  have  taps  at  5% 
increments  from  zero  to  full  inductance.  The  load  used  for  testing  was  arranged  as  a 
three-phase  Y-connected  load  with  an  open  neutral.  Hall  effect  sensors  were  used  to 
provide  current  feedback.  These  sensors  provide  an  output  current  which  is  1/1000  th  of 
the  measured  current.  By  using  a  load  resistor  this  current  was  converted  to  a  voltage 
signal,  which  is  more  easily  handled  by  the  dSPACE  system.  The  hall  effect  sensors 
require  an  external  grounded  voltage  supply.  This  fact  must  be  considered  when 
connecting  to  the  feedback  lines  to  avoid  ground  loops. 
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C.  SIMULINK  IMPLEMENTATION  OF  CONTROL  ALGORITHM 


The  stationary  reference  frame  current  control  algorithm  discussed  in  Chapter  IV 
was  programmed  in  SIMULINK.  Figure  6-3  shows  the  block  diagram  for  the  code. 


The  Current  Sensing  block  contains  the  A/D  conversion  of  the  feedback  signals 
from  the  Hall  Effect  sensors.  Signals  from  -10  to  10  volts  are  accepted  by  the  A/D 
converter,  outside  this  range  the  unit  saturates.  These  -10  to  10  volt  analog  values  are 
scaled  to  a  digital  signal  in  the  range  -1  to  1.  As  stated  before  care  must  be  taken  to 
ensure  that  a  ground  loop  does  not  exist  between  the  feedback  signals  and  the  power  for 
the  DS1 102  board.  To  protect  the  A/D  converter  from  overvoltage  two  10  volt  zener 
diodes  were  used  to  limit  the  feedback  signal  to  +10  volts.  The  Preprocessing  block 
scales  the  current  signals  so  that  the  digital  signal  matches  the  analog  value  of  current 
seen  by  the  circuit  and  removes  any  bias  in  the  signal.  Digital  filtering  can  be  performed 
in  this  block,  if  required,  by  using  any  one  of  the  built-in  SIMULINK  filters.  One  must 
be  careful  when  introducing  a  filter  in  that  the  filter  will  introduce  an  amplitude  scale  and 
phase  shift.  For  this  testing  a  filter  was  not  necessary.  In  the  Desired  Currents  block  the 
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desired  sinusoidal  a-b-c  phase  currents  are  generated  and  a  ramp-up  signal  is  generated 
for  use  in  start-up.  The  ramp  causes  a  linear  increase  in  the  desired  signal  from  0  to 
100%.  For  testing,  the  duration  of  the  ramp-up  was  10  seconds.  The  Ks  and  Ks-1  blocks 
perform  the  algebraic  transformation  into  and  from  the  stationary  reference  frame.  Both 
the  q-  and  d-axis  controller  blocks  contain  proportional-integral  controllers.  Limited 
integrators  are  used  to  prevent  wind-up.  The  Scale  to  +1  block  limits  the  input  to  the 
PWM  to  +1  and  also  scales  the  incoming  signal  so  that  a  maximum  commanded  value 
corresponds  to  +1  input  to  the  PWM  block.  The  PWM  unit  requires  a  +1  signal.  A  -1 
corresponds  to  a  duty  cycle  of  zero  and  +1  to  a  duty  cycle  of  100%.  The  unit  uses  the 
slave  DSP  chip  to  generate  a  modulating  signal  with  a  duty  cycle  specified  by  the  input. 
The  frequency  of  this  modulating  signal  can  be  set  to  one  of  four  values;  1.5kHz,  6kHz, 
25kHz  or  100kHz.  The  output  is  directed  through  channels  45, 46  and  41.  These 
channels  are  electrical  connections  within  the  DS1 102  board.  By  use  of  a  ribbon  cable 
connected  to  the  board  these  signals  are  brought  to  an  external  interface  where  they  can 
be  used  by  the  converter.  At  this  point  the  possibility  of  a  ground  loop  between  the 
external  circuit  and  the  DS1 102  board  exists  also.  Some  type  of  isolation  amplifier  is 
normally  used  to  protect  the  DSP  board.  The  details  of  each  block  are  included  in 
Appendix  F. 

Once  the  modulation  signals  were  passed  out  of  the  DS1 102  board  a  time  delay 
was  introduced  to  allow  one  switch  to  turn  off  before  the  opposite  switch  was  turned  on. 
This  delay  is  necessary  to  prevent  both  switches  being  closed  simultaneously  causing  a 
short  across  the  source.  This  was  implemented  in  hardware  using  dual  input  AND  and 
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NOR  gates  (with  an  RC  network  on  one  of  the  inputs)  and  LM3 1 1  comparators.  The 
value  of  the  time  delay  was  approximately  8psec,  this  seemed  to  work  well  with  the 
INVERPOWER  switches. 

D.  TESTING  OF  STATIONARY  REFERENCE  FRAME  CONTROLLER 

Values  used  for  testing  are  shown  in  Table  6-1 .  The  value  of  switching  frequency 
was  selected  based  upon  the  limitations  of  the  switches  used  in  the  INVERPOWER  unit. 

Table  6-1  -  Testing  Parameters 


Lf 

42.5mH 

^load 

29  Ohms 

Ki 

30 

Kp 

2 

«>e 

377  rad/sec 

f 

^switch 

1.5  kHz 

vDC 

100  V 

From  Chapter  IV,  it  is  known  that  a  steady-state  error  will  exist  between  the 
commanded  and  actual  currents  due  to  the  control  being  performed  in  the  stationary 
reference  frame.  To  predict  the  error  a  simulation  using  the  ACSL  code  for  the  current 
controller  was  performed.  Figure  6-4  shows  one  of  the  phase  currents  from  the 
simulation  with  a  commanded  current  amplitude  of  one  amp  peak  (two  amps  peak-to- 
peak). 

As  can  be  seen  in  Figure  6-4  a  substantial  steady-state  error  exists.  A  phase 
current  plot  from  the  hardware,  Figure  6-5,  verifies  this  error. 

Figures  6-4  and  6-5  show  a  very  good  correlation.  The  simulation  shows  a 
current  value  of  +.5  Amps  and  the  circuit  produced  a  value  of  about  +.55  Amps.  The 
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Figure  6-4  -  ACSL  Output  for  1  Amp  Commanded  Current 


Figure  6-5  -  Circuit  Output  for  1  Amp  Commanded  Current 
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difference  in  peak  magnitudes  can  be  attributed  to  the  delays  in  the  sensing  circuitry  and 
other  non-idealities  which  were  not  modeled  by  the  ACSL  simulation.  The  1.5  kHz 
harmonic  superimposed  on  the  60  Hz  fundamental  is  attributed  to  the  switching  and  is 
observed  on  both  waveforms. 

A  commanded  current  amplitude  of  1  amp  corresponds  to  synchronous  reference 
frame  commanded  values  of  i‘*s  =  1.0  and  i£  =  0.0 .  Equations  4-30  through  4-32 
(repeated  here  as  Equations  6-1  through  6-3)  can  be  used  to  predict  the  steady-state  error. 

(Eq.  6-1) 
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(Eq.  6-2) 
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(Eq.  6-3) 


The  values  K'p  and  in  Equations  6-2  and  6-3  include  the  forward  path  gains  of  the 

system  including  the  PWM  gain.  For  the  controller  as  implemented  the  forward  path  gain 
is  1/3.5,  since  this  was  the  value  used  in  the  scale  to  +1  block,  and  the  PWM  gain  is  the 
DC  voltage  divided  by  two.  This  gives  values  of  28.6  and  428.6  for  K'p  and  £/ , 
respectively.  The  steady-state  values  predicted  by  Equation  6-1  are  found  to  be 


ie  =  0.4604  and  =  0.1018 .  These  currents  can  be  used  directly  to  calculate  the  stator 


current  phasor  as  0.472Z12.5’ .  Thus,  the  actual  current,  as  predicted  by  Equations  6-1 


96 


through  6-3,  will  be  47%  of  the  magnitude  of  the  commanded  current  and  will  be  out  of 
phase  with  the  commanded  current  by  12.5  degrees. 

The  dSPACE  system  has  a  built  in  program  called  TRACE31 W  which  allows  for 
monitoring  of  internal  control  signals.  Using  the  TRACE3 1 W  program  the  values  of 
commanded  current  and  sensed  current  were  captured  and  shown  in  Figure  6-6.  The 
measured  current  is  the  signal  at  the  output  of  the  sensing  block.  The  important 
information  from  the  plot  is  that  a  phase  difference  of  about  10  degrees  exists  between 
the  desired  and  sensed  current.  The  presence  of  this  phase  offset  corresponds  to  the  value 
of  12.5  degrees  calculated  previously,  confirming  our  prediction.  The  amplitudes  of  the 
two  signals  cannot  be  compared  as  the  sensed  current  is  measured  before  any  processing. 
The  dc  offset  in  the  measured  signal  is  due  to  the  bias  of  the  hall  effect  sensor. 


The  commanded  current  was  raised  to  2  amps  and  the  output  from  the  simulation 
and  circuit  were  plotted  in  Figures  6-7  and  6-8  respectively.  Again  a  steady-state  error  of 
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0.5  Amps/Di  v 


Figure  6-8  -  Circuit  Output  for  2  Amp  Commanded  Current 
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approximately  50%  exists.  A  1.5  kHz  switching  harmonic  signal  can  be  seen  on  each  of 
the  plots. 

E.  TESTING  OF  SYNCHRONOUS  REFERENCE  FRAME  CONTROLLER 

The  synchronous  reference  frame  controller  was  then  implemented  in  SIMULINK 
and  merged  into  the  dSPACE  system.  Values  used  for  testing  the  synchronous  controller 
are  the  same  as  for  the  stationary  controller  except  that  the  DC  source  voltage  was  raised 
to  150V  to  maintain  operation  in  the  linear  range.  The  SIMULINK  block  diagram  for  the 
controller  is  shown  in  Figure  6-9.  The  details  of  the  blocks  are  included  in  Appendix  F. 


Figure  6-9  -  SIMULINK  Flowchart  for  Synchronous  Current  Controller 


It  was  shown  in  Chapter  4  that  a  controller  in  the  synchronous  reference  frame  can 
achieve  a  zero  steady-state  error.  An  ACSL  simulation  was  performed  with  the 
parameters  used  in  testing.  The  output  of  the  ACSL  code  (Figure  6-10)  shows  a  current 
waveform  of  one  amp  peak  value  for  a  one  amp  commanded  current  and  thus  a  zero 
steady-state  error.  Actual  circuit  output  for  one  amp  commanded  current  (Figure  6-1 1)  is 
shown  to  be  in  perfect  agreement  with  the  ACSL  simulation.  Figures  6-12  and  6-13 
illustrate  the  ACSL  and  circuit  output  for  a  commanded  current  of  two  amps. 
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Figure  6-10  -  ACSL  Output  for  1  Amp  Commanded  Current 


Figure  6-11  -  Circuit  Output  for  1  Amp  Commanded  Current 
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I.OAmps/Div  ^  |  |  Amps 


Figure  6-12  -  ACSL  Output  for  2  Amps  Commanded  Current 
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In  comparing  the  stationary  and  synchronous  reference  frame  controllers  it  is 
obvious  that  the  synchronous  controller  performs  much  better  with  regards  to  steady-state 
error.  Although  the  steady-state  error  of  the  stationary  controller  can  be  calculated,  and 
therefore  corrected  for,  the  error  varies  with  the  load.  Unless  the  load  is  known  precisely 
the  correction  may  not  be  accurate.  For  this  reason  the  synchronous  controller  is  the 
preferred  choice. 

One  drawback  of  the  synchronous  controller  is  that  it  is  more  computationally 
intensive  that  the  stationary  controller.  When  compiling  the  control  algorithm  and 
generating  the  C-code  to  be  used  by  the  DSP  chip  a  time  step  and  integration  algorithm 
must  be  specified.  These  are  specified  as  in  the  SIMULINK  environment.  The  dSPACE 
system  has  a  limit  on  how  many  calculations  can  be  performed  per  second.  As  the  time 
step  specified  becomes  small  this  limits  the  number  of  computations  which  can  occur 
during  each  step.  Complicated  algorithms  require  more  computations  than  simple  ones. 
So  as  the  algorithms  become  more  involved  the  time  step  at  which  they  can  be  run 
becomes  larger.  When  implementing  these  two  control  algorithms  in  the  dSPACE  using 
the  Euler  algorithm  it  was  found  that  the  minimum  time  step  for  the  stationary  algorithm 
was  .0002  seconds  corresponding  to  a  frequency  of  5  kHz,  and  for  the  synchronous 
algorithm  was  .0003  sec  corresponding  to  3.33  kHz.  This  difference  was  not  significant 
for  the  test  setup  used.  In  addition,  if  other  microprocessor  resources  are  used  to 
implement  these  algorithms  they  may  be  much  faster  and  capable  of  minimizing  the 
implementation  issues  between  the  two  current  regulator  algorithms. 
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To  compare  the  dynamic  responses  of  the  controllers  a  step  change  in  the 
commanded  current  from  one  to  two  amps  was  performed.  Using  the  TRACE3 1 W 
program  the  commanded  and  actual  q-axis  synchronous  current  were  captured.  Figures 
6-14  and  6-15  show  the  stationary  and  synchronous  controller  responses.  The  “noisy” 
signal  in  the  plots  represents  the  actual  q-axis  values  and  the  other  signal  is  the 
commanded  values.  The  response  times  of  both  controllers  are  comparable;  however  the 
stationary  controller  has  a  large  steady-state  error  as  shown  previously.  Overall  the 
synchronous  controller  outperformed  the  stationary  controller. 

The  measured  waveforms  in  Figures  6-14  and  6-15  appear  jagged.  This  is  a 
function  of  the  noise  present  in  the  signals  and  the  limitations  of  the  TRACE31W 
program. 
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Figure  6-15  -  Synchronous  Controller  Step  Response 


F.  SPACE  VECTOR  CONTROL  IMPLEMENTATION 

The  Space  Vector  Control  modulation  technique  was  not  tested,  however  this 
section  includes  a  brief  discussion  on  how  SVC  could  be  implemented  using  the  dSPACE 
system. 

A  S  VC-based  control  algorithm  with  feedback  could  be  implemented  in 
SIMULINK.  The  same  basic  set-up  of  sensing  and  output  channels  which  were  used  to 
test  the  current  regulation  algorithms  could  be  utilized  in  a  SVC  program.  The  PWM  unit 
would  be  replaced  by  a  Digital  to  Analog  converter,  as  SVC  does  not  utilize  sine-triangle 
PWM.  The  basic  steps  in  the  algorithm  would  include  generating  a  desired  voltage  signal 
(assuming  a  voltage  regulator  is  implemented),  generating  an  actual  voltage  signal  in  the 
correct  reference  frame,  comparing  the  two  voltage  signals  to  arrive  at  a  commanded 
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voltage  signal,  and  finally  converting  this  commanded  signal  into  the  desired  switching 


sequence. 

The  slowest  part  of  the  algorithm  will  be  in  converting  the  commanded  signal  into 
the  desired  output.  Since  none  of  the  slave  functions  of  the  DS1 102  board  can  be  utilized 
for  this  purpose  all  of  the  computations  must  be  performed  on  the  master  DSP.  This  will 
greatly  limit  the  frequency  at  which  a  SVC  algorithm  could  be  run  on  dSPACE. 

An  algorithm  based  on  a  look-up  table  would  be  faster,  but  it  is  not  readily 
implemented  using  the  dSPACE  system.  If  an  alternate  DSP  were  available  for 
programming,  a  reasonably  fast  algorithm  could  be  produced. 

This  chapter  has  undertaken  the  hardware  implementation  and  testing  of  various 
algorithms  on  a  hard-switched  converter.  The  following  chapter  investigates  the 
operation  of  an  actual  ARCP  converter. 
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VII.  ARCP  PROTOTYPE  UNIT 


A  prototype  ARCP  converter  was  constructed  by  the  Applied  Research 
Laboratory  /  Penn  State  University  (ARL/PSU).  This  converter  was  delivered  to  the 
Naval  Post  Graduate  School  for  incorporation  into  the  Power  Electronic  Building  Block 
Network  Evaluation  Simulator  (PEBBNES)  which  is  currently  under  construction. 
Funding  for  both  the  ARL/PSU  converter  and  the  PEBBNES  is  provided  from  the  Naval 
Surface  Warfare  Center.  This  chapter  discusses  the  prototype  ARCP  unit  parameters, 
operational  characteristics,  the  proposed  controller  to  be  used  with  the  prototype  ARCP 
unit,  and  recommendations  for  future  ARCP  evolution. 

A.  PHYSICAL  DESCRIPTION  OF  UNIT 

The  converter  consists  of  a  User  Interface  (UI)  board  and  three  identical  Pole 

boards  (basically  a  leg  of  the  inverter  as  described  in  Chapter  II)  mounted  on  a  heat 

sink/base  plate  made  of  3/8”  aluminum.  Overall  dimensions  of  the  unit  are  14”  wide,  24” 

deep  and  approximately  2  54”  high[19].  Restrictions  were  placed  on  the  size  to  allow 

incorporation  into  the  PEBBNES.  Complete  specifications  for  the  unit  are  listed  in 

Reference  19.  Figure  7-1  shows  a  photograph  of  the  entire  ARCP  prototype  unit. 

The  UI  board  contains  a  master  control  section  and  an  output  current  detection 

and  feedback  section  for  each  of  the  three  phases.  A  24  volt  DC  control  voltage  is 

required  by  the  UI  board,  as  well  as  by  each  of  the  Pole  boards.  Each  Pole  board  requires 

a  DC  bus  voltage,  with  a  minimum  output  capacitance  of  1 100  pF,  for  power  conversion. 

Optical  signals  are  required  by  the  Pole  boards  to  drive  the  main  switches,  one  control 
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Figure  7-1  -  ARCP  Prototype  Unit 

signal  for  each  switch.  Figure  7-2  shows  a  photograph  of  a  single  phase  of  the  ARCP 
including  all  of  the  major  connections  required.  The  main  and  auxiliary  switches  are 
located  below  the  board,  mounted  directly  onto  the  heat  sink,  with  the  electrical 
connections  on  the  board. 
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Figure  7-2  -  Single  Phase  of  ARCP  Prototype 

The  controller  interface  is  designed  to  be  electrically  isolated  from  the  ARCP  unit 
(by  use  of  optical  connections)  except  when  the  current  feedback  circuit  is  connected. 


For  this  reason  the  current  feedback  circuit  is  not  grounded  to  any  part  of  the  ARCP  unit. 
Thus,  the  ground  in  the  current  feedback  circuit  will  be  in  reference  to  the  controller 
ground. 

The  main  switches  used  in  the  prototype  unit  are  IRCPC50U  UltraFast  IGBT’s 
made  by  International  Rectifier.  They  have  rise  and  fall  times  rated  at  28  and  45  nsec 
respectively. [20]  With  switching  times  so  fast  the  IGBT’s  allow  for  very  high  frequency 
operation.  The  circuitry  necessary  to  turn  on  and  off  the  main  and  auxiliary  switches  is 
contained  entirely  on  the  Pole  board.  All  that  is  required  of  the  user  is  to  provide  24  volt 
power  and  an  optical  signal  to  the  main  switches.  The  boost  current  level  is  set  at 
approximately  one  amp  for  no  load  and  automatically  increases  to  approximately  16  amps 
at  full  load.  The  threshold  current  level  is  12  amps. 

B.  OPERATION 

Operation  of  the  prototype  ARCP  closely  mimics  the  ideal  circuit  operation 
presented  in  Chapter  II.  One  of  the  main  areas  of  concern  is  the  circuit  behavior  during 
commutation  between  states.  Figure  7-3  illustrates  the  transition  from  low  to  high 
voltage  state  with  positive  current,  this  corresponds  to  a  Case  1  condition  (diode 
conducting)  from  Chapter  II.  The  top  waveform  is  the  current  through  the  auxiliary 
circuit  and  the  bottom  waveform  is  output  voltage  of  the  ARCP  measured  across  the 
lower  main  switch. 

Figure  7-3  can  be  compared  to  the  ideal  commutation  of  Figure  2-10.  Load 
current  for  the  conditions  shown  in  Figure7-3  is  approximately  one  amp,  the  DC  supply 
voltage  is  25  volts. 
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Figure  7-3  -  ARCP  Commutation  from  Conducting  Diode 


Following  the  trace  of  the  current  waveform,  as  the  value  reaches  about  one  amp 
(the  load  current  value)  the  lower  main  switch  is  opened  and  the  circuit  begins  to 
resonate.  As  the  current  reaches  its  peak  value  (2.6  amps)  the  voltage  is  approximately  at 
half  of  the  final  value.  When  the  voltage  reaches  25  volts  the  upper  main  IGBT  is  turned 
on,  at  this  point  the  current  in  the  auxiliary  circuit  is  about  one  amp.  This  is  the  same 
value  as  the  load  current,  since  the  load  current  is  being  supplied  by  the  auxiliary  circuit 
at  this  point.  As  time  progresses  the  upper  IGBT  picks  up  the  load  and  the  auxiliary 
current  falls  to  zero,  at  which  point  the  auxiliary  switch  is  turned  off,  completing  the 
commutation.  The  major  discrepancies  between  the  real  and  ideal  case  is  the  small 
amount  of  ringing  seen  at  the  turn-on  of  the  main  IGBT  and  the  turn-off  of  the  auxiliary 
switches. 
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Commutation  with  the  switch  initially  conducting  is  illustrated  in  Figure  7-4,  this 
corresponds  to  a  Case  2  commutation.  A  Case  2  commutation  can  be  found  by  observing 
the  circuit  for  a  transition  from  upper  to  lower  voltage  states  with  positive  load  current, 
referenced  out  of  the  converter.  In  the  previous  condition  the  upper  switch  would  be 
conducting.  Load  conditions  are  the  same  as  before,  DC  voltage  25  volts  and  load 
current  one  amp.  The  ideal  case  for  the  conditions  represented  in  Figure  7-4  is  covered  in 
Chapter  II,  Figure  2-16.  Case  3  commutation  could  be  examined  by  increasing  the  load 
current  above  the  threshold  level,  and  monitoring  the  commutation  from  a  conducting 
switch. 


Figure  7-4  -  ARCP  Commutation  from  Conducting  Switch 
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Single-phase  testing  of  the  unit  was  performed  by  utilizing  the  dSPACE  system  to 
generate  modulating  signals  and  using  Hewlett  Packard  HFBR-1521  optical  transmitters 
and  HFBR-PNS005  fiber  optic  cable  to  interface  with  the  ARCP  prototype  unit.  A  delay 
circuit  feeding  LM3 1 1  comparators  was  used  to  drive  the  optical  transmitters.  The  delay 
circuit  was  required  to  achieve  a  lpsec  delay  between  switching  of  the  main  switches  (a 
specification  required  by  the  prototype  unit.) 

Although  three-phase  testing  was  not  completed,  due  primarily  to  delivery  delays, 
extension  of  the  testing  to  include  all  three  phases  is  straightforward.  In  particular,  it  will 
consist  of  building  two  additional  delay  circuits  identical  to  the  one  used  for  one-phase 
testing,  connecting  two  additional  outputs  from  the  dSPACE  system  to  drive  the 
additional  phases,  and  connecting  supply  and  control  power  to  the  additional  phases. 

The  effects  on  the  current  waveforms  of  changing  the  ARCP  switching  frequency 
and  of  changing  the  inductance  in  the  RL  load  may  be  observed  in  the  following  studies. 
With  a  switching  frequency  of  1.5kHz,  a  DC  supply  voltage  of  50  volts,  a  resistive  load 
of  16  Q,  and  a  filter  inductance  of  1  ImH,  the  current,  illustrated  in  Figure  7-5,  follows  a 
generally  sinusoidal  shape,  but  carries  quite  a  bit  of  switching  harmonics.  Increasing  the 
filter  inductance  to  values  of  43.5mH  and  86mH  while  maintaining  the  switching 
frequency,  supply  voltage  and  load  resistance  constant  is  illustrated  in  Figures  7-6  and  7- 
7.  The  larger  filter  values  reduce  the  harmonics  of  the  output  voltage,  but  the  amplitude 
of  the  voltage  is  reduced  due  to  the  larger  impedance  seen  by  the  ARCP  output. 
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Figure  7-5  -  ARCP  Open-Loop  Resistor  Voltage,  fswitch=1.5kHz,  L=1  IraH 
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Figure  7-7  -  ARCP  Open-Loop  Resistor  Voltage,  fswitch=1.5kHz,  L=86mH 

Figures  7-8  and  7-9  show  the  effect  of  increasing  the  switching  frequency  to  6kHz 
and  25kHz,  respectively.  The  DC  input  voltage  is  50  volts  and  the  load  resistance  is  16Q, 
as  before.  The  value  of  filter  inductance  is  1  lmH.  Comparing  the  voltage  waveforms 
with  Figure  7-5,  it  is  apparent  that  the  increased  switching  frequency  reduces  the 
harmonic  content  of  the  output  voltage.  However,  unlike  increasing  the  filter  inductance 
value,  higher  switching  frequencies  do  not  cause  a  lower  amplitude  response. 

Figures  7-5  through  7-9  are  generated  open-loop,  closed-loop  control  algorithms 
(such  as  those  presented  in  Chapter  IV)  can  also  be  used  to  reduce  the  distortion  seen  in 
the  output  voltage  waveform. 
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Figure  7-8  -  ARCP  Open-Loop  Resistor  Voltage,  fswitch=6kHz,  L=1  ImH 


Figure  7-9  -  ARCP  Open-Loop  Resistor  Voltage,  fswitch=25kHz,  L=1  ImH 
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C.  PROPOSED  CONTROLLER  FOR  ARCP 

The  Naval  Surface  Warfare  Center  has  constructed  a  controller  specifically 
designed  to  interface  with  the  prototype  ARCP  unit.  These  control  units  are  limited  in 
number  and  were  not  available  for  integration  with  the  ARCP  hardware  unit. 

The  controllers  are  built  around  Texas  Instrument’s  TMS320C30  DSP’s.  On  and 
off  chip  RAM  are  made  available  for  the  DSP’s  to  use  in  controlling  the  ARCP.  The 
capability  to  download  various  control  algorithms  into  the  DSP  is  possible.  Currently  the 
process  of  changing  control  algorithms  involves  reprogramming  EPROM’s,  but  it  is 
possible  through  further  development  to  download  a  new  algorithm  directly  from  a 
computer  into  the  DSP.  The  controller  can  become  very  useful  in  testing  various  control 
algorithms  as  the  capability  to  change  control  schemes  improves. 

Optical  output  is  provided  to  interface  with  the  ARCP,  four  optical  outputs  per 
phase  are  provided.  This  allows  the  controller  to  command  both  the  main  and  auxiliary 
switches  if  desired.  Electrical  connections  for  use  in  current  feedback,  voltage  feedback 
or  a  combination  of  current  and  voltage  feedback  control  schemes  are  provided.  A  total 
of  10  feedback  connections  are  possible.  Each  controller  is  capable  of  receiving  feedback 
from  and  providing  control  to  an  operating  three-phase  ARCP. 

Capability  to  optically  interconnect  controllers  is  provided,  this  may  lead  to 
parallel  ARCP  operation  or  connection  of  the  controller  to  a  master  control  unit.  The 
interconnection  of  controllers  will  allow  communication  between  units.  Thus,  the 
controllers  can  be  part  of  an  intelligence  network  built  into  the  power  system. 
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D.  RECOMMENDATIONS  FOR  FUTURE  ARCP  DEVELOPMENTS 

Currently  the  controller  used  to  drive  the  ARCP  must  generate  signals  for  the 
upper  and  lower  main  switches  and  ensure  that  the  proper  delay  between  signals  is 
achieved.  This  places  an  unnecessary  burden  on  the  controller.  Since  the  only  states 
which  the  ARCP  can  be  in  is  with  the  upper  main  switch  closed,  or  the  lower  main  switch 
closed,  or  in  a  transition  between  the  two  previous  states;  only  one  signal  is  needed  per- 
phase  for  the  ARCP.  With  one  signal,  a  logical  one  could  represent  the  upper  switch 
closed  and  a  logical  zero  the  lower  switch  closed.  The  necessary  delays  could  be  built 
into  the  ARCP  unit.  This  would  guarantee  a  more  efficient  commutation  process  for  the 
ARCP  since  the  switching  points  could  be  located  very  close  to  the  optimum  times.  In 
addition  the  timing  problem  between  the  ARCP  and  the  controller  would  be  eliminated, 
and  computation  time  used  by  the  controller  to  generate  two  separate  signals  with  delays 
could  be  used  to  perform  other  tasks.  The  effect  of  having  one  control  signal  per  phase 
would  have  the  effect  of  decentralizing  the  intelligence  in  the  power  network,  allowing 
the  system  to  operate  more  independently  and  efficiently. 

The  ability  to  change  the  boost  and  threshold  current  levels  would  allow  the 
ARCP’s  to  be  ‘tuned’  to  a  specific  load  condition.  This  should  be  made  an  option  so  that 
the  ARCP  could  operate  in  a  default  level  to  cover  all  operational  conditions. 

A  physical  improvement  which  would  improve  the  operability  of  the  ARCP 
would  be  to  provide  optical  feedback  to  the  controller  vice  electrical  feedback.  Electrical 
connections  typically  pick  up  noise,  especially  in  an  high  EMI  environment  such  as 
would  be  encountered  by  a  power  converter.  Optical  feedback  lines  would  allow  the 
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controller  to  be  moved  farther  from  the  ARCP  unit  away  from  EMI,  and  thus  improving 
performance.  In  addition  having  no  common  electrical  connections  between  the 
controller  and  the  ARCP  would  eliminate  the  possibility  of  a  fault  at  the  ARCP  affecting 
the  controller,  or  vice  versa. 

Providing  some  type  of  LC  filter  at  the  output  of  the  ARCP  would  allow  for 
‘clean’  power  to  be  provided  to  all  loads.  With  a  high  switching  frequency,  the  filter 
components  would  not  have  to  be  large,  and  may  not  be  needed  at  all.  However, 
providing  them  at  the  ARCP  output  may  eliminate  the  need  for  each  load  connected  to 
the  ARCP  to  provide  its  own  filtering. 

If  a  filter  were  connected  to  the  ARCP,  additional  feedback  signals  for  inductor 
current  and  output  capacitor  voltage  should  be  provided  to  the  controller  using  optical 
hardware.  This  would  allow  the  voltage  feedback  with  inner  current  control  loop 
algorithm  discussed  in  Chapter  IV  to  be  easily  implemented. 
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VIII.  CONCLUSIONS 


The  Auxiliary  Resonant  Commutated  Pole  Converter  provides  a  practical 
architecture  for  use  as  a  power  converter  in  shipboard  applications.  The  savings  provided 
by  increased  efficiency  outweigh  the  cost  of  the  added  complexity  of  ARCP  converters 
when  compared  to  hard-switched  converters.  This  savings  becomes  more  pronounced 
when  considering  very  high-frequency  switching  converters,  since  the  switching  losses  of 
conventional  hard-switched  converters  increase  proportionally  with  switching  frequency. 

It  was  shown,  through  detailed  simulation,  that  for  evaluation  of  control 
algorithms  the  ARCP  could  be  modeled  as  a  hard-switched  converter.  Although  losses  in 
a  hard  switched  inverter  are  higher,  the  behaviors  of  the  ARCP  and  hard-switched 
converter  are  very  similar  on  a  macroscopic  scale.  This  allowed  the  evaluation  and 
design  of  control  algorithms  using  simplified  hard-switched  converter  models,  reducing 
the  development  time. 

Several  control  algorithms  were  investigated  through  simulation  and  reduced- 
scale  circuit  modeling.  The  major  conclusion  from  the  study  of  these  algorithms  is  that 
when  controlling  time-varying  quantities,  such  as  sinusoidal  voltages  and  currents,  the 
steady-state  response  of  the  controllers  have  a  zero  error  when  the  control  is  performed  in 
a  reference  frame  where  the  controlled  quantities  are  DC  values.  For  the  cases  studied, 
the  synchronous  reference  frame  provided  DC  quantities  to  be  controlled.  Although 
slightly  more  computationally  intensive,  control  in  the  synchronous  reference  frame  is  the 
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preferred  method  when  dealing  with  control  algorithms  implemented  using  Digital  Signal 
Processors. 

In  Chapter  V  a  new  type  of  pulse-width-modulation  technique,  Space  Vector 
Control  (SVC),  was  investigated  and  simulated.  Several  advantages  of  SVC  over  sine- 
triangle  modulation  were  pointed  out;  most  significant  is  the  capability  of  the  SVC 
algorithm  to  realize  a  waveform  15%  larger  than  the  sine-triangle  control  without  going 
into  overmodulation.  No  hardware  testing  was  performed  using  this  modulation 
technique.  It  is  recommended  that  a  study  be  performed,  using  either  an  ARCP  or  hard- 
switched  converter,  which  compares  traditional  sine-triangle  modulation  with  SVC 
modulation.  An  efficient  implementation  of  the  SVC  algorithm  may  prove  to  be  the  one 
drawback  of  Space  Vector  Control  which  makes  it  challenging  for  use  in  a  DSP-based 
controller. 

The  ARCP  prototype  unit  discussed  in  Chapter  VII  provides  an  excellent  tool  for 
evaluating  the  ARCP  interaction  in  a  shipboard  DC  power  distribution  network.  With 
some  modifications  for  power  connection  and  control  input  this  unit  can  be  incorporated 
into  the  Power  Electronics  Building  Block  Network  Evaluation  Simulator  (PEBBNES) 
currently  under  construction  at  the  Naval  Postgraduate  School. 

Recommendations  for  improvements  to  the  ARCP  prototype  unit,  given  in 
Chapter  VII,  focus  on  simplifying  the  control  interface  (recommending  all  fiber  optic 
connections  for  control  and  feedback)  and  freeing  computation  time  on  the  controller 
(recommending  only  one  control  signal  for  each  phase.) 
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APPENDIX  A.  LAPLACE  TRANSFORM  SOLUTION  OF  ARCP  RESONANT 

STATE  DYNAMICS 

This  appendix  develops  the  equations  for  the  resonant  phase  of  commutation  for 
the  ARCP  converter.  LaPlace  transforms  are  used  to  arrive  at  a  solution.  Figure  A-l 
shows  the  equivalent  circuit  for  the  analysis.  The  filter  inductor  is  considered  a  constant 
current  source  during  commutation  and  the  DC  capacitors  are  considered  constant  voltage 
sources. 


Figure  A-l  -  Equivalent  Circuit  during  Resonant  Phase 


Loop  Equations 

Loopl:  T-F«-L.(|-t)=0 

L0°P2:  lL-L'&-f')-F«=0  (E*A'2) 

Defining 

*'r=*Wi  (Eq.A-3) 
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Loop  1  becomes 


Loop  2  becomes 


%-La-Lr^0 

2  C1  r  dt 


—■  +  L  -j1-  -VC2-0 
2  r  dt  C2 


Subtracting  Equation  A-5  from  A-4 


dir 

Va-VC2-2Lr-^  =  0 


Defining 


V  =  V 

y  C2  v  out 


Using  the  relationship 


Vci  becomes 


V  +V  -V 

VC\^VC2  ydc 


V  =  V  -V  =V  - V 

yC\  ydc  yC 2  ydc  y  out 


And  Equation  A- 10  results 


Vdc-Fo„,-Foo,-2Ir^  =  0 


V  ,  -  — +  L  —  =  0 

°u,  2  r  dt 


Node  Equation 

Summing  currents  at  the  node  between  VC1  and  VC2 


.  .  .  Cr  dVa  Cr  dVC2 

,  -i  -l,  +— - — - : —  =  0 

2  1  L  2  dt  2  dt 


(Eq.  A-4) 

(Eq.  A-5) 

(Eq.  A-6) 

(Eq.  A-7) 

(Eq.  A-8) 

(Eq.  A-9) 


(Eq.  A- 10) 


(Eq.  A-ll) 
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Substituting  the  expression  for  ir  from  Equation  A-3 


c AdVn 

2  \  dt 


dVc£\ 

dt  ) 


+  1, 


iL=0 


Simplifying  gives  Equation  A- 13 

Cjd^-Vj  d\C 

2  dt  dt  j 


(Eq.  A- 12) 


'Cr-^L  +  ir_iL=0  (Eq.  A- 13) 

LaPlace  Transform  Solution 

The  initial  Conditions  assuming  commutation  from  low  to  high  voltage  with  diode  D2 


conducting  are  given  by 

n*(0>  =  0  (Eq.  A- 14) 

i,(0)  =  iL  +1^  (Eq.  A-15) 

Taking  the  LaPlace  transform  of  Equation  A- 10  yields 

n.,(j)-^-  +  [2L,(sI,(s)-it(0))]  =  0  (Eq.  A-16) 

Substituting  in  the  initial  conditions  results  in 

ym  W-^  +  [2Lt(sl,(s)-it  =  0  (Eq.  A-17) 

Taking  the  LaPlace  transform  of  Equation  A- 13  gives 

-  C,[sE„„  (s)-V„  (0)]  +  l,(s)- ^  =  0  (Eq.  A-18) 

Substituting  in  the  initial  conditions  results  in 

-C,sF„„(j)  +  Ir(s)-!J-  =  0  (Eq.A-19) 
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Solving  for  Vout 

Solving  Equation  A- 19  for  Ir(s) 

irW  =  — +c,sr„(s) 

s 

Substituting  into  Equation  A- 17 


(Eq.  A-20) 


(Eq.  A-21) 


Simplifying  Equation  A-21 


LW-^+l,C,s!V.(S)-Lrito,  =0 


^,(s)(l  +  L,C,s2)  =  ^-Lri 


boost 


dc 


boost 


(\  +  LrCrs2)ls  l+LrCrs2 


1  K 


dc 


1 


,  _  L  C  2  L  C 

VouXs)=  7-p 


r  boost ) 


A 


L  C  +S  J 

'  r  r  ' 


1 


LrCr 


■  +  s 


(Eq.  A-22) 


Taking  the  inverse  LaPlace  transform 


K.  (<)■=■ y 


1-cos 


'  i  ' 

t 


AT7 


:Lri  boost  sinl 


VLA  rboost  IVla 


(Eq.  A-23) 


Defining  the  resonant  frequency  asco0  =  ,  = . ,  Equation  A-23  can  be  rewritten  as 

V^r  Cr 


(Eq.  A-24) 
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Solving  for  I, 


Solving  Equation  A- 17  for  V0 


Kut(s) =  ^‘-Lr(sIr(s)-iL  -iboost) 


(Eq.  A-25) 


Substituting  Equation  A-25  into  Equation  A- 19 


-C,J  -^-L,(sI,(s)-iL-iloo,,)  +I,(s)-7  =  0 


(Eq.  A-26) 


Simplifying  Equation  A-26 


-  c  L,slr(s)  +  LriL  +  L,i_,  +  I,(s) =  0 


1 

-  +  LrCrs2Ir  (s)  -  LrCrs(iL  +  i^, )  +  It  (s)  -  ^  =  0 


Ir(s)(l  +  LrCrs2)-  ^  +  dc2  r  +  LrCrs(iL  +iboost) 


j  /  \  _  _ j_L _  ^dc^-r  LrCrs(iL  +  iboost) 

r  s(l  +  LrCrs2)  2(l  +  LrCrs2)  (l  +  LrCrs2) 


Ir(s)  = 


LrCr  FdcCr  + 


2rr+s' 


(Eq.  A-27) 


Taking  the  inverse  LaPlace  transform  yields  an  expression  for  the  resonant  current 


1  _  .  1 

(  /  ,  Y\  /T  r  ^dc^rS,n  rr~r  t  (  \ 

ir(t)  =  iL  1-cos  —j=^=t  + — — - - - —  +0l  +  'boost)cos  -/f==t  (Eq‘ A’28) 

V  WL  A,  JJ  1  WLrCr  J 
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Substituting  coQ 


,  Equation  A-28  may  be  rewritten  as 


v  Ic 

ir(t)  =  hi1  ~  cos6>0t)  +  fl^sin*>0t  +  (iL  +  ib0( 


Combining  terms  results  in  Equation  A-30 


)cos<u0t  (Eq.  A-29) 


(Eq.  A-30) 
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APPENDIX B.  ACSL  SIMULATION  PROGRAMS 


A.  SINGLE-PHASE  HARD  SWITCHED  INVERTER  SIMULATION 

PROGRAM 

INITIAL 


" - This  Program  simulates  operation  of  a  hard  switched  inverter  with 

"  sine-triangle  PWM.  Ideal  switching  and  components  are  assumed. 
"  Load  consists  of  an  LC  filter  with  a  resistive  load. 

"  DC  Voltage  =  500  Volts 
"  Switching  Frequency  =  6KHz 
"  Output  Frequency  =  60Hz 


MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint=1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.03 

CONSTANT  Lf=10.1e-3 
CONSTANT  Cf=2.0e-3 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  ilic=15.0 
CONSTANT  VLoadic=250.0 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1 .0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*Vtripk*ftriang 
we=60.0*2.0*3. 14159265 
Vout=0.0 

fl 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 


"maximum  integration  step  size  " 

"min  time  step:  for  var  step  algorithm" 
"data  communication  interval 
"integration  algorithm— R.K.  4th  " 

"stop  point  for  integration  " 

"filter  (output)  inductance  " 

"filter  (output)  capacitance  " 

"DC  source  voltage  " 

"load  resistance  " 

"initial  load  current  " 

"initial  load  voltage  " 

"reference  triangle  frequency  " 

"reference  triangle  peak  value 
"reference  triangle  period  " 

"reference  triangle  half  period 
"reference  triangle  slope 
"radian  frequency  of  control  sinusiod" 
"initial  output  voltage 
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DERIVATIVE 


" — sine-triangle  modulation" 
tt=mod(t,Ttri) 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0 . 0+slope  *  (tt-T  o  ver2) 

END  IF 

Vsin  =  5.0*cos(we*t) 

SCHEDULE  hi.XN.(Vtri  -  Vsin) 
SCHEDULE  lo.XP.(Vtri  -  Vsin) 

"-—calculate  derivatives  of  state  variables" 
iLoad=VLoad/Rload 
PVLoad=(iL-iLoad)/Cf 
piL=(V  out- VLoad)/Lf 

" — calculate  state  variables" 

VLoad=integ(pVLoad,VLoadic) 

iL=integ(piL,iLic) 

END  !  "of  derivative" 

DISCRETE  hi 
Vout=Vdc 

END  !  "of  discrete" 

DISCRETE  lo 
Vout=0.0 

END  !  "of  discrete" 

END  !  "of  dynamic" 

END  !  "of  program" 


!  "form  ref  triangle 


!  "control  sinusiod 
! "switch  to  HI  state 
!  "switch  to  LO  state 
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B.  SINGLE-PHASE  ARCP  SIMULATION 


PROGRAM 

INITIAL 


" - This  Program  simulates  operation  of  the  Auxiliary  Resonant  " 

"  Commutated  Pole  (ARCP)  converter.  Ideal  switching  and  components" 
"  are  assumed.  Load  consists  of  an  LC  filter  with  a  resistive  load.  " 
"  Resonant  Frequency  =  250KHz  " 

"  DC  Voltage  =  500  Volts 

"  Switching  Frequency  =  6KHz  " 

"  Output  Frequency  =  60  Hz  " 


MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.03 


!  "maximum  integration  step  size  " 

!"min  time  step:for  var  step  algorithm" 
!  "data  communication  interval  " 

!  "integration  algorithm-R.K.  4th  " 

!  "stop  point  for  integration  " 


LOGICAL  SI 

Sl=.false. 

S2=.true. 

Al=.false. 

A2=.false. 

state  1=.  false. 

state2=.false. 

state3=.false. 

state4=.false. 

state5=.false. 

state6=.false. 


,S2,A1  ,A2,statel,state2,state3,state4,state5,state6 

!  "status  of  upper  main  switch,  switch  #1 
!"  status  of  lower  main  switch,  switch  #2 
!  "status  of  auxiliary  switch  #1 
!"  status  of  auxiliary  switch  #2 
!  "transition  Lo  to  Hi,  IL  pos 
!  "transition  Hi  to  Lo,  IL  pos  and  >  Ithres 
!  "transition  Hi  to  Lo,  IL  pos  and  <  Ithres 
!  "transition  Lo  to  Hi,  IL  neg  and  <  Ithres 
!  "transition  Hi  to  Lo,  IL  neg 
!  "transition  Lo  to  Hi,  IL  neg  and  >  Ithres 


ti 


ft 


CONSTANT  Lr=2.9e-6 
CONSTANT  Cr=.3e-6 
CONSTANT  Lf=10.1e-3 
CONSTANT  Cf=2.0e-3 
CONSTANT  RLoad=25.0 
CONSTANT  Vdc=500.0 
CONSTANT  iLic=15.0 
CONSTANT  Voutic=0.0 
CONSTANT  VLoadic=250.0 
CONSTANT  iric=0.0 
CONSTANT  iboostswitch=4.0 


!  "resonant  inductance 

!  "resonant  cap  (total  both  switches) 

["filter  (output)  inductance 

["filter  (output)  capacitance 

["load  resistance 

!"DC  source  voltage 

["initial  load  current 

["initial  output  voltage 

["initial  load  voltage 

["initial  auxiliary  circuit  current 

["boost  current  reqd  before  switching 
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CONSTANT  ithreshold=18 
wo=l  .0/(sqrt(Lr*Cr)) 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1.0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*Vtripk*fitriang 
we=60.0*2.0*3. 141 59265 


("threshold  current 
!  "resonant  radian  frequency 
("reference  triangle  frequency  " 

!  "reference  triangle  peak  value  " 

("reference  triangle  period 
("reference  triangle  half  period 
("reference  triangle  slope 
("radian  frequency  of  control  sinusiod" 


END  ("of initial" 


DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — sine-triangle  modulation" 

tt=mod(t,Ttri)  ("form  ref  triangle 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0.0+slope*  (tt-T  over2) 


ENDIF 

Vsin  =  5.0*cos(we*t) 

("control  sinusiod 

it 

SCHEDULE  hi.XN.(Vtri  -  Vsin) 

("switch  to  HI  state 

ft 

SCHEDULE  lo.XP.(Vtri  -  Vsin) 

("switch  to  LO  state 

it 

" — schedule  switching  points  of  specific  components" 

SCHEDULE  S 1  on.XP.  V  out-V  dc 

("SI  on 

It 

SCHEDULE  S 1  offS2.XN.ir+iboostswitch 

("SI  off  State2 

tt 

SCHEDULE  S 1  offS5.XN.(iboostswitch-iL+ir) 

("SI  off  State5 

ti 

SCHEDULE  S2on.XN.Vout 

!"S2  on 

II 

SCHEDULE  S2offS  1  ,XN.(iboostswitch+iL-ir) 

!"S2  off  State  1 

It 

SCHEDULE  S2offS4.XP.ir-iboostswitch 

!"S2  off  State4 

it 

SCHEDULE  Aloff.XP.ir 

!"A1  off 

ti 

SCHEDULE  A2off.XN.ir 

!"A2  off 

It 

" — determine  characteristic  equations  for  switch  configuration" 
IF(.not.(Sl.or.S2).and.(Al  ,or.A2))  THEN 
p V  out=(ir-iL)/ Cr 
ELSE  IF  (,not.(Sl.or.S2))  THEN 
pVout=-iL/Cr 
ELSE 
pVout=0.0 
ENDIF 


4- 
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IF  (Al.or.A2)  THEN 
pir=(V  dc-2.0*Vout)/(2.0*Lr) 

ELSE 

pir=0.0 

ENDIF 

" — calculate  diode  and  switch  currents" 
IF(Sl.and.((iL-ir).GT.O.O))  THEN 
iSl=iL-ir 
ELSE 

iDl=-(iL-ir) 

ENDIF 

IF(S2.and.((iL-ir).LT.0.0))  THEN 
iS2=-(iL-ir) 

ELSE 

iD2=iL-ir 

ENDIF 

" — calculate  state  variables" 

Vout=integ(pVout,Voutic) 

ir=integ(pir,iric) 

iLoad=VLoad/RLoad 

pVLoad=(iL-iLoad)/Cf 

VLoad=integ(pVLoad,VLoadic) 

piL=(Vout-VLoad)/Lf 

iL=integ(piL,iLic) 

END  !  "of  derivative" 

" — determine  state  of  converter  during  switching  interval" 
DISCRETE  lo 
IF(iL.LT.O.O)  THEN 
state5=.true. 

Al=.  true. 

ELSE  IF (iL  .LT.ithreshold)  THEN 
state2=.true. 

Al=.true. 

ELSE 

state3=.true. 

Sl=.false. 

ENDIF 

END  !  "of  discrete" 
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DISCRETE  hi 
IF(iL.GT.O.O)  THEN 
state  l=.true. 

A2=.true. 

ELSE  IF ( AB S (iL) .LT.ithreshold)  THEN 
state4=.true. 

A2=.true. 

ELSE 

state6=.true. 

S2=.false. 

ENDIF 

END  !  "of  discrete" 

" — determine  which  switches  are  closed" 
DISCRETE  Sion 
Sl=.true. 
pVout=0.0 
Vout=Vdc 
state  l=.false. 
state4=.false. 
state6=.false. 

END  !  "of  discrete" 

DISCRETE  SloffS2 
IF(state2)  THEN 
Sl=.false. 

ENDIF 

END  !  "of  discrete" 

DISCRETE  SloffS5 
IF(state5)  THEN 
Sl=.  false. 

END  IF 

END  !  "of  discrete" 

DISCRETE  S2on 
S2=.  true. 
pVout=0.0 
Vout=0.0 
state2=.false. 
state3=.false. 
state5=.false. 

END  !  "of  discrete" 
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DISCRETE  S2ofifSl 
IF(statel)  THEN 
S2=.false. 

END  IF 

END  !  "of  discrete" 

DISCRETE  S2offS4 
IF(state4)  THEN 
S2=.false. 

END  IF 

END  !  "of  discrete" 

DISCRETE  Aloff 
Al=.  false. 
pir=0.0 
ir=0.0 

END  !  "of  discrete" 

DISCRETE  A2off 
A2=.false. 
pir=0.0 
ir=0.0 

END  !  "of  discrete" 
END  !  "of  dynamic" 


END  !  "of  program' 


C.  THREE-PHASE  HARD-SWITCHED  INVERTER  SIMULATION 

PROGRAM 

INITIAL 

" - This  Program  simulates  operation  of  a  3-Phase  hard  switched  inverter" 

"  with  sine-triangle  PWM.  Ideal  switching  and  components  are  assumed" 
"  Load  consists  of  an  LC  filter  with  a  resistive  load. 

"  DC  Voltage  500  Volts 
"  Switching  Frequency  =  6kHz 
"  Output  Frequency  =  60Hz 

MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.03 
CONSTANT  Lfr=10.1e-3 
CONSTANT  CF=1.0e-4 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1.0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*Vtripk*ftriang 
we=60.0*2.0*3. 14159265 
pi2o3=2.0*3. 14159265/3.0 
CONSTANT  iqssic=8.0 
CONSTANT  idssic=3.0 
CONSTATN  A=9.0 
iLa=iqssic 

ilb=-0. 5  *  iqssic-idssic*  sqrt(3 . 0)/2 .0 
ilc=-0. 5  *  iqssic+idssic*  sqrt(3 . 0)/2 . 0 
CONSTANT  VqssLic=200.0 
CONSTANT  VdssLic=0.0 
VLa=VqssL 

VLb=-0.5*VqssLic-(sqrt(3.0)/2.0)*VdssLic 
VLc— 0 . 5  *  VqssLic+(sqrt(3 . 0)/2 .0)*  V  dssLic 
Vag=0.0 
Vbg=0.0 
Vcg=0.0 

END  !  "of  initial" 


M 

M 

tl 

!  "maximum  integration  step  size 
!"min  time  step:for  var  step  algorithm" 

!  "data  communication  interval 
! "  integration  algorithm-R.K.  4th 

! "  stop  point  for  integration  " 

!  "filter  (output)  inductance 

!  "filter  (output)  capacitance 

!  "DC  source  voltage 

!  "load  resistance 

!  "reference  triangle  frequency 

!  "reference  triangle  peak  value 

!  "reference  triangle  period 

!  "reference  triangle  half  period 

!  "reference  triangle  slope 

!  "radian  frequency  of  control  sinusiod" 

! "  1 20  degrees  in  radians 

["initial  q  axis  current 

!  "initial  d  axis  current 

! "  control  voltage  magnitude 
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DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — sine-triangle  modulation" 

tt=mod(t,Ttri)  ["form  ref  triangle  " 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0.0+slope*  (tt-Tover2) 

END  IF 

!  "control  sinusiod  " 
!  "switch  A  to  HI  state  " 
!  "switch  A  to  LO  state" 

!  "control  sinusoid  " 
!  "switch  B  to  HI  state  " 
!  "switch  B  to  LO  state" 

!  "control  sinusiod  " 
!"  switch  C  to  HI  state  " 
!  "switch  C  to  LO  state" 

" — calculate  line  voltages" 

V  as=(2 . 0/3 . 0)  *  Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 

Vbs=(2.0/3.0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 

Vcs=(2.0/3.0)*Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

" — transform  to  stationary  reference  frame" 
Vqss=(2.0/3.0)*(Vas-0.5*Vbs-0.5*Vcs) 

Vdss=(sqrt(3 .0)/3.0)*(-Vbs+Vcs) 

" — find  current  derivatives  in  stationary  frame" 
piqss=(V  qss-VqssL)/Lf 
pidss=(V  dss- V  dssL)/Lf 

" — find  load  voltage  derivatives  in  stationary  frame" 
pVqssL=(iqss-(V  qssL/Rload))/Cf 
pVdssL=(idss-(V  dssL/Rload))/Cf 

"——calculate  stationary  currents" 
iqss=integ(piqss,iqssic) 
idss=integ(pidss,idssic) 
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Vsina  =  A*cos(we*t) 

SCHEDULE  hia.XN.(Vtri  -  Vsina) 
SCHEDULE  loa.XP.(Vtri  -  Vsina) 

Vsinb  =  A*cos(we*t-pi2o3) 
SCHEDULE  hib.XN.(Vtri  -  Vsinb) 
SCHEDULE  lob.XP.(Vtri  -  Vsinb) 

Vsinc  =  A*cos(we*t+pi2o3) 
SCHEDULE  hic.XN.(Vtri  -  Vsinc) 
SCHEDULE  loc.XP.(Vtri  -  Vsinc) 


" — calculate  stationary  voltages" 

V  qssL=integ(pV  qssL,  V  qssLic) 

V  dssL=integ(pV  dssL,  VdssLic) 

"--transform  back  to  ABC  quantaties" 
iLa=iqss 

iLb=-0.5*iqss-(sqrt(3.0)/2.0)*idss 
iLc=-0.5  *  iqss+(sqrt(3 .0)/2.0)*idss 

VLa=VqssL 

VLb=-0. 5  *  VqssL-(sqrt(3 . 0)/2.0)  *  VdssL 
VLc=-0. 5  *  V  qssL+(sqrt(3 ,0)/2.0)*  VdssL 

iLoada=VLa/RLoad 

iLoadb=VLb/RLoad 

iLoadc=VLc/Rload 

END  !  "of  derivative" 

DISCRETE  hia 
Vag=Vdc 

END  !  "of  discrete" 

DISCRETE  loa 
Vag=0.0 

END  !  "of  discrete" 

DISCRETE  hib 
Vbg=Vdc 

END  !  "of  discrete" 

DISCRETE  lob 
Vbg=0.0 

END  !  "of  discrete" 

DISCRETE  hie 
Vcg=Vdc 

END  !  "of  discrete" 

DISCRETE  loc 
Vcg=0.0 

END  !  "of  discrete" 


END  !  "of  dynamic 
END  !  "of  program 
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D.  TRHEE-PHASE  ARCP  SIMULA  TION 

PROGRAM 

INITIAL 


"  This  Program  simulates  operation  of  the  Three  Phase  Auxiliary 
"  Resonant  Comutated  Pole  (ARCP)  converter. 

Ideal  switching  and  components  are  assumed 
"  Resonant  Frequency  =  250KHz 
"  DC  Voltage  =  500  Volts 
"  Switching  Frequency  =  6KHz 
Output  Frequency  =  60Hz 


MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.00003 


!  "maximum  integration  step  size  " 

!"min  time  step:  for  var  step  algorithm" 
.'"data  communication  interval  " 

!  "integration  algorithm-R.K.  4th  " 

! "stop  point  for  integration  " 


LOGICAL  S 1  a,S2a,Ala,A2a,state  1  a,state2a,state3a,state4a,state5a,state6a 


Sla=.  false. 

S2a=.true. 

Ala=.false. 

A2a=.false. 

state  la=.false. 

state2a=.false. 

state3a=.false. 

state4a=.false. 

state5a=.false. 

state6a=.false. 


!  "status  of  upper  main  switch,  switch  #  1 
!  "status  of  lower  main  switch,  switch  #2  " 

!  "status  of  auxiliary  switch  #1  " 

!  "status  of  auxiliary  switch  #2  " 

!  "transition  Lo  to  Hi,  IL  pos 
!  "transition  Hi  to  Lo,  IL  pos  and  >  Ithres 
!  "transition  Hi  to  Lo,  IL  pos  and  <  Ithres 
!  "transition  Lo  to  Hi,  IL  neg  and  <  Ithres  " 
!  "transition  Hi  to  Lo,  IL  neg 
!  "transition  Lo  to  Hi,  IL  neg  and  >  Ithres  " 


LOGICAL  S 1  b,S2b, A 1  b,A2b,state  1  b,state2b,state3b,state4b,state5b,state6b 


Sib- false. 
S2b=.  true. 
Alb=.false. 
A2b=.false. 


!  "status  of  upper  main  switch,  switch  #1 
!  "status  of  lower  main  switch,  switch  #2 
!  "status  of  auxiliary  switch  #1 
.'"status  of  auxiliary  switch  #2 


state  lb=.false. 
state2b=.false. 
state3b=.false. 


!  "transition  Lo  to  Hi,  IL  pos  " 

!  "transition  Hi  to  Lo,  IL  pos  and  >  Ithres  " 
!  "transition  Hi  to  Lo,  IL  pos  and  <  Ithres 
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state4b=.false. 

state5b=.false. 

state6b=.false. 


!  "transition  Lo  to  Hi,  IL  neg  and  <  Ithres 
!  "transition  Hi  to  Lo,  IL  neg 
!  "transition  Lo  to  Hi,  IL  neg  and  >  Ithres 


LOGICAL  S 1  c,S2c, A1  c,A2c,state  1  c,state2c,state3  c,state4c,state5c,state6c 


["status  of  upper  main  switch,  switch  #1 
["status  of  lower  main  switch,  switch  #2 
["status  of  auxiliary  switch  #1 
["status  of  auxiliary  switch  #2 
["transition  Lo  to  Hi,  IL  pos 
["transition  Hi  to  Lo,  IL  pos  and  >  Ithres 
["transition  Hi  to  Lo,  IL  pos  and  <  Ithres 
["transition  Lo  to  Hi,  IL  neg  and  <  Ithres 
["transition  Hi  to  Lo,  IL  neg 
["transition  Lo  to  Hi,  IL  neg  and  >  Ithres 


Sic*,  false. 

S2c=.  true. 

Alc=.  false. 

A2c=.false. 
state  lc=.false. 
state2c=. false. 
state3c=.false. 
state4c=.false. 
state5c=.false. 
state6c=.false. 

CONSTANT  Lr=2.9e-6 
CONSTANT  Cr=.3e-6 
CONSTANT  Lf=10.1e-3 
CONSTANT  Cf=T.0e-4 
CONSTANT  RLoad=25.0 
CONSTANT  Vdc=500.0 
CONSTANT  iqssic=5.0 
CONSTANT  idssic=3.0 
CONSTANT  Vagic=0.0 
CONSTANT  Vbgic=0.0 
CONSTANT  Vcgic=0.0 
CONSTANT  iraic=0.0 
CONSTANT  irbic=0.0 
CONSTANT  ircic=0.0 
CONSTANT  iboostswitch=5.0 
CONSTANT  ithreshold=18 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1.0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*  Vtripk*  ftriang 
we=60.0*2.0*3. 141 59265 
pi2o3=2.0*3.141 59265/3.0 
CONST ATN  A=5.0 
iLa=iqssic 

ilb=-0 . 5  *  iqssic-idssic*  sqrt(3 .0)/2 .0 
ilc=-0 . 5  *  iqssic+idssic*  sqrt(3 .0)/2 .0 
CONSTANT  VqssLic=200.0 


["resonant  inductance 

["resonant  cap  (total  both  switches) 

["filter  (output)  inductance 

["filter  (output)  capacitace 

["load  resistance 

!"DC  source  voltage 

["initial  q  axis  current 

["initial  d  axis  current 

[  "initial  output  current 

["initial  output  current 

[  "initial  output  current 

["initial  auxiliary  circuit  current 

[  "initial  auxiliary  circuit  current 

["initial  auxiliary  circuit  current 

["boost  current  reqd  before  switching" 

["threshold  current 

["reference  triangle  frequency 

["reference  triangle  peak  value 

["reference  triangle  period 

["reference  triangle  half  period 

["reference  triangle  slope 

["radian  frequency  of  control  sinusiod' 

["120  degrees  in  radians 

["control  voltage  magnitude 
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CONSTANT  VdssLic=0.0 
VLa=VqssL 

VLb=-0.5*VqssLic-(sqrt(3.0)/2.0)*VdssLic 

VLc=-0.5*VqssLic+(sqrt(3.0)/2.0)*VdssLic 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 
DERIVATIVE 


" — sine-triangle  modulation" 

tt=mod(t,Ttri)  !  "form  ref  triangle  " 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*  slope 
ELSE 

Vtri=-1 0.0+slope*(tt-Tover2) 

END  IF 


Vsina  =  A*cos(we*t) 

SCHEDULE  hia.XN.(Vtri  -  Vsina) 
SCHEDULE  loa.XP.(Vtri  -  Vsina) 


!  "control  sinusioid  " 
["switch  A  to  HI  state  " 
["switch  A  to  LO  state" 


Vsinb  =  A*cos(we*t-pi2o3) 
SCHEDULE  hib.XN.(Vtri  -  Vsinb) 
SCHEDULE  lob.XP.(Vtri  -  Vsinb) 


["control  sinusoid  " 
["switch  B  to  HI  state  " 
["switch  B  to  LO  state" 


Vsinc  =  A*cos(we*t+pi2o3) 
SCHEDULE  hic.XN.(Vtri  -  Vsinc) 
SCHEDULE  loc.XP.(Vtri  -  Vsinc) 


["control  sinusiod  " 
["switch  C  to  HI  state  " 
["switch  C  to  LO  state" 


" — schedule  switching  points  of  specific  components  for  Phase  A  " 


SCHEDULE  S 1  aon.XP.  Vag-Vdc  !  "S 1  on 

SCHEDULE  SI aoff.XN.switcha  ! "SI  off 

SCHEDULE  S2aon.XN.Vag  !"S2  on 

SCHEDULE  S2aoff.XN.switcha  !"S2ofF 

SCHEDULE  Alaoff.XP.ira  !"A1  off 

SCHEDULE  A2aoff.XN.ira  !"A2off 


139 


" _ schedule  switching  points  of  specific  components  for  Phase  B  " 


SCHEDULE  Slbon.XP.Vbg-Vdc  !"S1  on 

SCHEDULE  SlboffXN.switchb  !"S1  off 

SCHEDULE  S2bon.XN.Vbg  '■  "S2  on 

SCHEDULE  S2boff.XN.switchb  !  "S2  off 

SCHEDULE  Alboff.XP.irb  !"A1  off 

SCHEDULE  A2boff.XN.irb  ! "  A2  off 

•'.—schedule  switching  points  of  specific  components  for  Phase  C 
SCHEDULE  Slcon.XP.Vcg-Vdc  !"S1  on 

SCHEDULE  Slcoff.XN.switchc  !"S1  off 

SCHEDULE  S2con.XN.Vcg  !"S2  on 

SCHEDULE  S2coff.XN .  switchc  !"S2  off 

SCHEDULE  Alcoff.XP.irc  !"A1  off 

SCHEDULE  A2coff.XN.irc  '■ "  A2  off 


"-—generate  control  for  main  switches 

switcha=(ira+iboostswitch)*state2a  +  (iboostswitch-iLa+ira)*state5a  +  & 
(iboostswitch+iLa-ira) *  state  1  a  -  (ira-iboostswitch)*state4a 
switchb=(irb+iboostswitch)*state2b  +  (iboostswitch-iLb+irb)*state5b  +  & 
(iboostswitch+iLb-irb)*  state  1  b  -  (irb-iboostswitch)*state4b 
switchc=(irc+iboostswitch)*state2c  +  (iboostswitch-iLc+irc)*state5c  +  & 
(iboostswitch+iLc-irc)  *  state  1  c  -  (irc-iboostswitch)*state4c 

" _ determine  characteristic  equations  for  switch  configuration  phase  A" 

IF(.not.(S  1  a.or.S2a).and.((Al  a.or.  A2a)))  THEN 
p  V  ag=(ira-iLa)/ Cr 
ELSE  IF  (.not.(Sla.or.S2a))  THEN 
pVag=-iLa/Cr 
ELSE 
pVag=0.0 
ENDIF 

IF  ((Ala.or.A2a))  THEN 
pira=(Vdc-2.0*Vag)/(2.0*Lr) 

ELSE 

pira=0.0 

ENDIF 

" — calculate  diode  and  switch  currents" 

IF(S  1  a.and.((iLa-ira).GT.O.O))  THEN 
iSla=iLa-ira 
ELSE 

iDla=-(iLa-ira) 

ENDIF 
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IF(S2a.and.((iLa-ira).LT.0.0))  THEN 
iS2a=-(iLa-ira) 

ELSE 

iD2a=iLa-ira 

ENDIF 

" — determine  characteristic  equations  for  switch  configuration  phase  B" 
IF(.not.(Slb.or.S2b).and.((Alb.or.A2b)))  THEN 
pVbg=(irb-iLb)/Cr 
ELSE  IF  (.not.(Slb.or.S2b))  THEN 
PVbg=-iLb/Cr 
ELSE 
pVbg=0.0 
ENDIF 

IF  ((Alb.or.A2b))  THEN 
pirb=(Vdc-2.0*Vbg)/(2.0*Lr) 

ELSE 

pirb=0.0 

ENDIF 

" — calculate  diode  and  switch  currents" 

IF(Slb.and.((iLb-irb).GT.O.O»  THEN 
iSlb=iLb-irb 
ELSE 

iDlb=-(iLb-irb) 

ENDIF 

IF  (S2b  .and.  ((iLb-irb)  .LT.0 . 0»  THEN 
iS2b=-(iLb-irb) 

ELSE 

iD2b=iLb-irb 

ENDIF 

"-—determine  characteristic  equations  for  switch  configuration  phase  C" 
IF(.not.(Slc.or.S2c).and.((Alc.or.A2c)))  THEN 
pV  cg=(irc-iLc)/ Cr 
ELSE  IF  (.not.(SIc.or.S2c))  THEN 
pVcg=-iLc/Cr 
ELSE 
pVcg=0.0 
ENDIF 
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IF  ((Alc.or.A2c))  THEN 
pirc=(V  dc-2.0*  Vcg)/(2.0*Lr) 

ELSE 

pirc=0.0 

ENDIF 

" — calculate  diode  and  switch  currents" 

IF(S  1  c.and.((iLc-irc).GT.O.O))  THEN 
iSlc=iLc-irc 
ELSE 

iDlc=-(iLc-irc) 

ENDIF 

IF(S2c.and.((iLc-irc).LT.0.0))  THEN 
iS2c=-(iLc-irc) 

ELSE 

iD2c=iLc-irc 

ENDIF 

" — calculate  state  variables  phase  A" 
Vag=integ(pV  ag,  V  agic) 
ira=integ(pira,iraic) 

" — -calculate  state  variables  phase  B" 

Vbg=integ(pVbg,Vbgic) 

irb=integ(pirb,irbic) 

" — calculate  state  variables  phase  C" 

Vcg=integ(pVcg,Vcgic) 

irc=integ(pirc,ircic) 

" — calculate  line  voltages" 

Vas=(2.0/3.0)*Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 
Vbs=(2.0/3.0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 
Vcs=(2. 0/3.0)* Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

" — transform  to  stationary  reference  frame" 

Vqss=(2.0/3.0)*(Vas-0.5*Vbs-0.5*Vcs) 

Vdss=(sqrt(3.0)/3.0)*(-Vbs+Vcs) 

" — find  current  derivatives  in  stationary  frame" 
piqss=(V  qss- V  qssL)/Lf 
pidss=(V  dss-V  dssL)/Lf 


* 
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"-—find  load  voltage  derivatives  in  stationary  frame" 
p  V  qssL=(iqss-(V  qssL/Rload))/Cf 
pV  dssL=(idss-(V  dssL/Rload))/Cf 

"-—calculate  stationary  currents" 

iqss=integ(piqss,iqssic) 

idss=integ(pidss,idssic) 

"——calculate  stationary  voltages" 

V  qssL=integ(pV  qssL,  V  qssLic) 

V  dssL=integ(pV  dssL,  V  dssLic) 

" — transform  back  to  ABC  quantaties" 
iLa=iqss 

iLb— 0 .5  *  iqss-(sqrt(3 . 0)/2 .0)*  idss 
iLc=-0. 5  *  iqss+(sqrt(3 .0)12 .0)  *  idss 

VLa=VqssL 

VLb=-0.5*VqssL-(sqrt(3.0)/2.0)*VdssL 

VLc=-0.5*VqssL+(sqrt(3.0)/2.0)*VdssL 

"-—calculate  load  currents" 
iLoada=VLa/RLoad 
iLoadb=VLb/RLoad 
iLoadc=VLc/RLoad 
END  !  "of  derivative" 

"—determine  state  of  phase  A  during  switching  interval" 
DISCRETE  loa 
IF(iLa.LT.0.0)  THEN 
state5a=.true. 

Ala=.  true. 

ELSE  IF(iLa.LT.ithreshold)  THEN 
state2a=.true. 

Ala=.  true. 

ELSE 

state3a=.true. 

Sla=.false. 

ENDIF 

END  !  "of  discrete" 
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DISCRETE  hia 
IF(iLa.GT.O.O)  THEN 
state  la=.true. 

A2a=.true. 

ELSE  IF (AB  S(iLa).LT.  ithreshold)  THEN 
state4a=.true. 

A2a=.true. 

ELSE 

state6a=.true. 

S2a=.  false. 

ENDIF 

END  !  "of  discrete" 

" — determine  state  of  phase  B  during  switching  interval" 
DISCRETE  lob 
IF(iLb.LT.O.O)  THEN 
state5b=.true. 

Alb=.true. 

ELSE  IF(iLb.LT.ithreshold)  THEN 
state2b=.  true. 

Alb=.true. 

ELSE 

state3b=.true. 

Slb=.false. 

ENDIF 

END  !  "of  discrete" 

DISCRETE  hib 
IF(iLb.GT.O.O)  THEN 
state  lb=.true. 

A2b=.true. 

ELSE  IF(ABS(iLb).LT.ithreshold)  THEN 
state4b=.true. 

A2b=.true. 

ELSE 

state6b=.true. 

S2b=.false. 

ENDIF 

END  !  "of  discrete" 
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"-—determine  state  of  phase  C  during  switching  interval" 
DISCRETE  loc 
IF(iLc.LT.O.O)  THEN 
state5c=.true. 

Alc=.true. 

ELSE  IF (iLc .LT.ithreshold)  THEN 
state2c=.true. 

Alc=.true. 

ELSE 

state3c=.true. 

Slc=.false. 

ENDIF 

END  !  "of  discrete" 

DISCRETE  hie 
IF(iLc.GT.O.O)  THEN 
statelc=.true. 

A2c=.true. 

ELSE  IF ( AB S(iLc) .LT.ithreshold)  THEN 
state4c=.true. 

A2c=.true. 

ELSE 

state6c=true. 

S2c=  false. 

ENDIF 

END  !  "of  discrete" 

" — determine  which  switches  are  closed  Phase  A" 
DISCRETE  Slaon 
Sla=.true. 
pVag=0.0 
Vag=Vdc 
statela=.  false. 
state4a=.false. 
state6a=.false. 

END  !  "of  discrete" 

DISCRETE  Slaoff 
Sla=.false. 

END  !  "of  discrete" 
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DISCRETE  S2aon 
S2a=.true. 
pVag=0.0 
Vag=0.0 
state2a=.false. 
state3a=.false. 
state5a=. false. 

END  !  "of  discrete" 

DISCRETE  S2aoff 

S2a=.false. 

END  !  "of  discrete" 

DISCRETE  Alaoff 
Ala=.  false. 
pira=0.0 
ira=0.0 

END  !  "of  discrete" 

DISCRETE  A2aoff 
A2a=.false. 
pira=0.0 
ira=0.0 

END  !  "of  discrete" 

" — determine  which  switches  are  closed  Phase  B" 

DISCRETE  Slbon 
Slb=.true. 
pVbg=0.0 
Vbg-Vdc 
state  lb=.false. 
state4b=.false. 
state6b=.false. 

END  !  "of  discrete" 

DISCRETE  Slboff 

Slb=.false. 

END  !  "of  discrete" 

DISCRETE  S2bon 
S2b=.true. 
pVbg=0.0 
Vbg=0.0 
state2b=.false. 
state3b=.false. 
state5b=.false. 

END  !  "of  discrete" 
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DISCRETE  S2boff 
S2b=.false. 

END  !  "of  discrete' 


DISCRETE  Alboff 
Alb=.false. 
pirb=0.0 
irb=0.0 

END  !  "of  discrete" 

DISCRETE  A2boff 
A2b=.false. 
pirb=0.0 
irb=0.0 

END  !  "of  discrete" 

"-—determine  which  switches  are  closed  Phase  C" 

DISCRETE  SI  con 
Slc=.true. 
pVcg=0.0 
Vcg=Vdc 
state  lc=.false. 
state4c=.false. 
state6c=.  false. 

END  !  "of  discrete" 

DISCRETE  Slcoff 

Slc=.false. 

END  !  "of  discrete" 

DISCRETE  S2con 
S2c=.true. 
pVcg=0.0 
Vcg=0.0 
state2c=.false. 
state3c=.false. 
state5c=.false. 

END  !  "of  discrete" 

DISCRETE  S2coff 

S2c=.false. 

END  !  "of  discrete" 
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DISCRETE  Alcoff 
Alc=.false. 
pirc=0.0 
irc=0.0 

END  !  "of  discrete" 

DISCRETE  A2coff 
A2c=.  false. 
pirc=0.0 
irc=0.0 

END  !  "of  discrete" 
END  !  "of  dynamic" 


END  !  "of  program' 


APPENDIX  C.  TOTAL  HARMONIC  DISTORTION  CALCULATION 


When  evaluating  total  harmonic  distortion  it  is  assumed  that  the  fundamental 
frequency  is  known.  Although  it  is  not  absolutely  necessary  to  know  the  fundamental 
frequency  a  priori,  it  must  at  least  be  estimated  to  perform  the  calculation.  Since  the 
estimation  can  become  extremely  complicated  our  assumption  simplifies  the  task. 

Total  harmonic  distortion  is  defined  as 


v 2  +v2  +v2  +• 

THD%  =  100 x  si - *3.  *4 

V  K] 


=  100x 


y 

ifTAK, 


(Eq.  C-l) 


where  V  indicates  an  RMS  voltage,  and  h  indicates  the  order  of  the  harmonic. 

To  evaluate  the  harmonic  components  the  principle  of  orthogonal  functions  will 
be  used.  The  following  discussion  shows  the  development  of  the  procedure. 

Initially  the  quantity  in  question,  in  this  case  voltage,  is  defined  in  Equation  C-2. 


v,=  2>„sin(n<V  +  0„) 

n=l,3,5,... 

Using  a  trigonometric  expansion  Equation  C-2  becomes 


(Eq.  C-2) 


oo 

vs  =  ^F„[sin(nfi)0r)cos^n +cos(nfi)0/l)sin^„]  (Eq.  C-3) 


n=l,3,5,... 


Multiplying  vs  by  sin(w6i0t)  and  integrating  over  a  fundamental  period,  the  function  C 
is  obtained 


2 

c  =  -  jvs  sin (ma0t)dt 
1  A 


(Eq.  C-4) 
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Substituting  Equation  C-3  into  Equation  C-4,  Equation  C-5  is  obtained  where  C,  and  C2 
are  given  by  Equations  C-6  and  C-7. 

c-|  j>X(C,+C2)  (Eq.  C-5) 

1  n=  1,3,5,... 

T 

C,  =  Jsin(w®0t)cos#„  sin(mco0t)dt  =  0  for  m^n  (Eq.  C-6) 

o 

T 

C2  -  jcos(ncoot)cos0n  sin (mco0t)dt  =  0  for  all  m  and  n  (Eq.  C-7) 

o 

Reducing  Equation  C-5,  Equation  C-8  is  obtained 

2  t 

C  =  —  Jcos^„  sin 2(nco0t)dt  (Eq.  C-8) 

T  o 

Integrating  and  solving,  Equation  C-8  becomes 

C  =  V„  cos  0 n  (Eq.  C-9) 

The  same  procedure  is  repeated  again,  except  vs  is  initially  multiplied  by 
cos(ma>0t) ,  this  will  produce  a  sine  function  result.  The  results  are  summarized  as 


2  T 

s  =  - Jv,  cos(mco0t)dt 
*  0 

(Eq.  C-10) 

s  =  I  i>„(S,+S2) 

(Eq.  C-ll) 

T 

Sl  =  Jsin(n6)or)cos0„  cos(ma>Qt)dt  =  0  for  all  m  and  n 

0 

(Eq.  C-12) 

T 

S2  =  Jcos(tt£y0r)sin#B cos{mo)0t)dt  =  0  for  m*n 

(Eq.  C-13) 

0 
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(Eq.  C-14) 


2  1 

s  =  -v„  Jsin0„  cos 2(nco0t)dt 
1  0 

S  =  Vn  sin#„  (Eq.C-15) 

With  both  the  sine  and  cosine  functions  available,  the  magnitude  of  the  nth  harmonic,  Vn , 
and  the  respective  phase  angle,  9 n  ,  can  be  evaluated  using  Equations  C-l 6  and  C-17. 

Vn=4s2  +C2  (Eq.C-16) 

^=tan’‘(f)  (Eq.  C-17) 


The  procedure  is  shown  graphically  in  Figure  C-l. 


Figure  C-l  -  Harmonic  Analyzer 
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APPENDIX  D.  CONTROLLER  SIMULATION  PROGRAMS 


A.  MATLAB  SIMULATION  FOR  VOLTAGE  FEEDBACK  CONTROLLER 

%  This  program  simulates  a  stationary  voltage  feedback  controller  in  the  synchronous 
%  reference  frame.  A  step  response  is  generated  and  the  poles  and  zeros  of  the  transfer 
%  function  along  with  the  steady  state  matrix  are  output. 

Kp=7  %  Proportional  Gain 

Ki=50  %  Integral  Gain 

U=50;  %  Step  Response  Magnitude 

t=linspace(0, . 03,20000); 

%%  Simulation  Parameters  %% 

CfHOOOl; 

Lf=  01; 

Rload=25; 
we=2*pi*60; 

%%  System  Matrices  %% 

A=[-l/(Cf*Rload)  -we  1/Cf 

we  -l/(Cf*Rload)  0 

-(Kp+1)/Lf  0  0 

0  -(Kp+1)/Lf  we 

-10  0 

0-10 
B=[l  0  Kp/Lf  0  0  0;0  1  0  Kp/Lf  0  0]'; 

C=U*[1  0  0  0  0  0]; 

D=[0  0]; 

%%  Poles  and  Zeros  of  Transfer  Function  %% 

[num,den]=ss2tf(A,B,C,D,  1 ); 

R=roots(den) 

P=roots(num) 

%%  Steady  State  Matrix  %% 

SS=-[1  0  0  0  0  0;0  1  0  0  0  0]*inv(A)*B 

%%  Step  Response  Output  %% 
step(A,B,C,D,l,t),title('MATLAB  Step  Response') 


0  0  0 

1/Cf  0  0 

-we  Ki/Lf  0 

0  0  Ki/Lf 

0  0  -we 

0  we  0]; 
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B.  ACSL  SIMULATION  FOR  VOLTAGE  FEEDBACK  CONTROLLER 


PROGRAM 

INITIAL 

"  This  Program  simulates  operation  of  a  3-Phase  hard  switched  inverter 
"  with  sine-triangle  PWM.  Ideal  switching  and  components  are  assumed 
"  Load  consists  of  an  LC  filter  with  a  resistive  load. 


"  DC  Voltage  500  Volts 
"  Switching  Frequency  =  6kHz 
MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.03 
CONSTANT  Lf=10.1e-3 
CONSTANT  CF=1.0e-4 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1.0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*Vtripk*ftriang 
we=60.0*2.0*3. 14159265 
pi2o3=2.0*3. 141 59265/3.0 
CONSTANT  iqssic=8.0 
CONSTANT  idssic=3.0 


Output  Frequency  =  60Hz 
["maximum  integration  step  size 
!"min  time  step:for  var  step  algorithm 
["data  communication  interval 
["integration  algorithm— R.K.  4th 

["stop  point  for  integration 

["filter  (output)  inductance 

["filter  (output)  capacitance 

!"DC  source  voltage 

["load  resistance 

["reference  triangle  frequency 

["reference  triangle  peak  value 

["reference  triangle  period 

["reference  triangle  half  period 

["reference  triangle  slope 

["radian  frequency  of  control  sinusiod 

["120  degrees  in  radians 

["initial  q  axis  current 

["initial  d  axis  current 


iLa=iqssic 


iLb=-0.5  *  iqssic-idssic*  sqrt(3 .0)/2.0 
iLc=-0.5*iqssic+idssic*sqrt(3.0)/2.0 
CONSTANT  VqssLic=250.0 
CONSTANT  VdssLic=0.0 
VLa=VqssL 

VLb=-0.5*VqssLic-(sqrt(3.0)/2.0)*VdssLic 

VLc=-0.5*VqssLic+(sqrt(3.0)/2.0)*VdssLic 

Vag=0.0 

Vbg=0.0 

Vcg=0.0 


CONSTANT  Kl=10.0 
CONSTANT  K2=0.0 
CONSTANT  A=250.0 
Vmax=(V  dc/20.0) 
END  !  "of  initial" 


"proportional  gain 
"integral  gain 

"desired  voltage  magnitude 
"linear  pulse  width  modulation  gain 


DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — sine-triangle  modulation" 

tt=mod(t,Ttri)  !"  form  ref  triangle  " 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0.0+slope*(tt-Tover2) 

END  IF 

" — calculate  synchronous  reference  frame  angle" 
thetae=mod(integ(we,0.0),twopi) 

" — calculate  desired  voltages" 

Vdesa  =  A*cos(we*t)  ! "desired  phase  A  voltage  " 

Vdesb  =  A*cos(we*t-pi2o3)  !  "desired  phase  B  voltage  " 

Vdesc  =  A*cos(we*t+pi2o3)  !  "desired  phase  C  voltage  " 

" — transform  desired  voltages  to  qdO  quantities" 

Vqsdes=(2.0/3.0)*(Vdesa-0.5*Vdesb-0.5*Vdesc)  ["desired  q-axis  voltage" 
Vdsdes=(sqrt(3.0)/3.0)*(-Vdesb+Vdesc)  ["desired  d-axis  voltage" 

"-—determine  error  voltages" 

V  qserr= V  qsdes-VqssL 

V  dserr=V  dsdes- V  dssL 

"——generate  commanded  voltages  (PI  Controller)" 

Vqscom=(Kl*Vqserr  +  K2*integ(Vqserr,0.0)) 

Vdscom=(Kl  *  Vdserr  +  K2*integ(V  dserr,0.0)) 

" — transform  commanded  quantaties  to  ABC  frame" 

Vacom=Vqscom 

Vbcom=-0.5*Vqscom-(sqrt(3.0)/2.0)*Vdscom 

V  ccom=-0 . 5  * V  qscom+(sqrt(3 .0)/2.0)*Vdscom 

" — normalize  voltages  for  use  in  switching  control" 

V  acont=V  aconvV  max 
Vbcont=Vbcom/V  max 

V  ccont=V  ccomA^  max 

"—transform  desired  Voltages  to  synchronous  qdO  quantities" 
Vqsedes=(2.0/3 .0)*(V  desa*cos(thetae)+Vdesb*cos(thetae-pi2o3)+& 
Vdesc*cos(thetae+pi2o3)) 
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Vdsedes=(2.0/3.0)*(Vdesa*sin(thetae)+Vdesb*sin(thetae-pi2o3)+& 

Vdesc  *  sin(thetae+pi2o3 )) 

" — transform  actual  currents  to  synchronous  qdO  quantities" 
Vqse=(2.0/3.0)*(VLa*cos(thetae)+VLb*cos(thetaepi2o3)+& 

VLc*  cos(thetae+pi2o3 )) 

Vdse=(2.0/3.0)*(VLa*sin(thetae)+VLb*sin(thetae-i2o3)+VLc*sin(thetae+pi2o3)) 

" — schedule  inverter  switching" 

SCHEDULE  hia.XN.(Vtri  -  Vacont) 

SCHEDULE  loa.XP.(Vtri  -  Vacont) 

SCHEDULE  hib.XN.(Vtri  -  Vbcont) 

SCHEDULE  lob.XP.(Vtri  -  Vbcont) 

SCHEDULE  hic.XN.(Vtri  -  Vccont) 

SCHEDULE  loc.XP.(Vtri  -  Vccont) 

" — calculate  line  voltages" 

Vas=(2.0/3.0)*Vag  -  (1 .0/3.0)*(Vbg  +  Vcg) 

Vbs=(2.0/3.0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 

Vcs=(2.0/3.0)*Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

" — transform  to  stationary  reference  frame" 
Vqss=(2.0/3.0)*(Vas-0.5*Vbs-0.5*Vcs) 

V  dss=(sqrt(3 . 0)/3 .0)*(-Vbs+Vcs) 

" — find  current  derivatives  in  stationary  frame" 
piqss=(V  qss-VqssL)/Lf 
pidss=(V  dss-VdssL)/Lf 

" — find  load  voltage  derivatives  in  stationary  frame" 
pV  qssL=(iqss-(V  qssL/Rload))/ Cf 
pVdssL=(idss-(V  dssL/Rload))/Cf 

" — calculate  stationary  currents" 
iqss=integ(piqss,iqssic) 
idss=integ(pidss,idssic) 

" — calculate  stationary  voltages" 

V  qssL=integ(p  V  qssL,V  qssLic) 

V  dssL=integ(pV  dssL,  V  dssLic) 


!"  switch  A  to  HI  state" 

! "switch  A  to  LO  state" 

["switch  B  to  HI  state  " 
["switch  B  to  LO  state" 

["switch  C  to  HI  state  " 
["switch  C  to  LO  state" 
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" — transform  back  to  ABC  quantaties" 
iLa=iqss 

iLb=-0.5*iqss-(sqrt(3.0)/2.0)*idss 

iLc=-0.5*iqss+(sqrt(3.0)/2.0)*idss 

VLa=VqssL 

VLb=-0.5*VqssL-(sqrt(3.0)/2.0)*VdssL 

VLc=-0.5*VqssL+(sqrt(3.0)/2.0)*VdssL 

iLoada=VLa/RLoad 

iLoadb=VLb/RLoad 

iLoadc=VLc/RLoad 

END  !  "of  derivative" 

DISCRETE  hia 
Vag=Vdc 

END  !  "of  discrete" 

DISCRETE  loa 
Vag=0.0 

END  !  "of  discrete" 

DISCRETE  hib 
Vbg=Vdc 

END  !  "of  discrete" 

DISCRETE  lob 
Vbg=0.0 

END  !  "of  discrete" 

DISCRETE  hie 
Vcg=Vdc 

END  !  "of  discrete" 

DISCRETE  loc 
Vcg=0.0 

END  !  "of  discrete" 

END  !  "of  dynamic" 

END  !  "of  program" 
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C.  ACSL  SIMULATION  FOR  CURRENT  FEEDBACK  CONTROLLER  IN  THE 
STATIONARY  REFERENCE  FRAME 

PROGRAM 


INITIAL 

"  This  Program  simulates  operation  of  a  3-Phase  hard  switched  inverter " 

"  with  sine-triangle  PWM.  Ideal  switching  and  components  are  assumed." 
"  Load  consists  of  an  LC  filter  with  a  resistive  load. 

"  A  Controller  operating  in  the  stationary  reference  frame  is  included 
"  DC  Voltage  =  500  Volts  " 

"  Switching  Frequency  =  6KHz  Output  Frequency  =  60Hz 


MAXTERVAL  maxt=  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.032 
CONSTANT  Lf=10.1e-3 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1 .0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4 . 0  *  V  tripk*  ftriang 
we=60.0*2.0*3. 14159265 
pi2o3=2.0*3. 14159265/3.0 
CONSTANT  twopi  =  3.14259265 
CONSTANT  iqssic=5.0 
CONSTANT  idssic=3.0 
iLa=iqssic 

ilb=-0 . 5  *  iqssic-idssic*  sqrt(3 .0)/2 .0 

ilc=-0.5*iqssic+idssic*sqrt(3.0)/2.0 

Vag=0.0 

Vbg-0.0 

Vcg=0.0 

CONSTANT  Kl=0.10 
CONSTANT  K2=0.001 
CONSTANT  A=9.0 
V  max=(V  dc/20 .0) 

END  !  "of  initial" 


["maximum  integration  step  size 
!"min  time  step:for  var  step  algorithm  " 
!  "data  communication  interval 
["integration  algorithm— R.K.  4th 

!  "stop  point  for  integration 

["filter  (output)  inductance 

!"DC  source  voltage 

["load  resistance 

["reference  triangle  frequency 

!  "reference  triangle  peak  value 

["reference  triangle  period 

!  "reference  triangle  half  period 

["reference  triangle  slope 

["radian  frequency  of  contrlo  sinusiod" 

["120  degrees  in  radins 

!  "360  degrees  in  radians 

!  "initial  q  axis  current 

!  "initial  d  axis  current 


["proportional  gain 
["integral  gain 
["desired  current  magnitude 
["pulse  width  modulation  gain 
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DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — sine-triangle  modulation" 

tt=mod(t,Ttri)  ["form  ref  triangle  " 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0.0+slope*(tt-Tover2) 

END  IF 

"-—calculate  synchronous  reference  frame  angle" 
thetae=integ(we,0.0) 

"-—calculate  desired  currents" 

Idesa  =  A*cos(thetae)  ! "desired  phase  A  current  " 

Idesb  =  A*cos(thetae-pi2o3)  ! "desired  phase  B  current  " 

Idesc  =  A*cos(thetae+pi2o3)  ! "desired  phase  C  current  " 

" — transform  desired  currents  to  stationary  qdO  quantities" 
Iqsdes=(2.0/3.0)*(Idesa-0.5*Idesb-0.5*Idesc) 

Idsdes=(sqrt(3 .0)/3 .0)*  (-Idesb+Idesc) 

"—transform  desired  currents  to  synchronous  qdO  quantities" 
Iqsedes=(2.0/3.0)*(Idesa*cos(thetae)+Idesb*cos(thetae-pi2o3)& 
+Idesc*cos(thetae+pi2o3)) 

Idsedes=(2.0/3.0)*(Idesa*sin(thetae)+Idesb*sin(thetae-pi2o3)& 

+Idesc*sin(thetae+pi2o3)) 

" — transform  actual  currents  to  synchronous  qdO  quantities" 

Iqse=(2.0/3.0)*(ILa*cos(thetae)+ILb*cos(thetae-pi2o3)+ILc*cos(thetae+pi2o3)) 

Idse=(2.0/3.0)*(ILa*sin(thetae)+ILb*sin(thetae-pi2o3)+ILc*sin(thetae+pi2o3)) 

" — determine  error  currents  in  stationary  reference  frame" 

Iqserr=Iqsdes-iqss 

Idserr=Idsdes-idss 

"-—generate  commanded  currents  (PI  Controller)" 

Iqscom=(Kl*Iqserr  +  K2*integ(Iqserr,0.0)) 

Idscom=(Kl*Idserr  +  K2*integ(Idserr,0.0)) 
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"-—transform  commanded  quantaties  to  ABC  frame" 
Iacom=Iqscom 

Ibcom=-0.5*Iqscom  -  (sqrt(3.0)/2.0)*Idscom 
Iccom=-0.5*Iqscom  +  (sqrt(3.0)/2.0)*Idscom 


"——normalize  voltages  for  use  in  switching  control" 


V  acont=Iacom*  Rload/V  max 
Vbcont=Ibcom*  Rload/V  max 

V  ccont=Iccom*  Rload/V  max 

"—schedule  inverter  switching" 
SCHEDULE  hia.XN.(Vtri  -  Vacont) 
SCHEDULE  loa.XP.(Vtri  -  Vacont) 

! "switch  A  to  HI  state" 
!  "switch  A  to  LO  state' 

SCHEDULE  hib.XN.(Vtri  -  Vbcont) 
SCHEDULE  lob.XP.(Vtri  -  Vbcont) 

!  "switch  B  to  HI  state  ' 
! "switch  B  to  LO  state' 

SCHEDULE  hic.XN.(Vtri  -  Vccont) 
SCHEDULE  loc.XP.(Vtri  -  Vccont) 

!  "switch  C  to  HI  state  ' 
! "switch  C  to  LO  state' 

" — calculate  line  voltages" 

Vas=(2.0/3.0)*Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 
Vbs=(2.0/3.0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 
Vcs=(2.0/3.0)*Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 


" — transform  to  stationary  reference  frame" 
Vqss=(2.0/3.0)*(Vas-0.5*Vbs-0.5*Vcs) 

V  dss=(sqrt(3 . 0)/3 . 0)  *  (- Vbs+ V  cs) 

"——find  current  derivatives  in  stationary  frame" 
piqss=(V  qss-RLoad*iqss)/Lf 
pidss=(V  dss-Rload*  idss)/Lf 

" — calculate  stationary  currents" 

iqss=integ(piqss,iqssic) 

idss=integ(pidss,idssic) 

" — transform  back  to  ABC  quantaties" 
iLa=iqss 

iLb=-0.5*iqss-(sqrt(3.0)/2.0)*idss 

iLc=-0.5*iqss+(sqrt(3.0)/2.0)*idss 

END  !  "of  derivative" 
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DISCRETE  hia 
Vag=Vdc 

END  !  "of  discrete" 

DISCRETE  loa 
Vag=0.0 

END  !  "of  discrete" 

DISCRETE  hib 
Vbg=Vdc 

END  !  "of  discrete" 

DISCRETE  lob 
Vbg=0.0 

END  !  "of  discrete" 

DISCRETE  hie 
Vcg=Vdc 

END  !  "of  discrete" 

DISCRETE  loc 
Vcg=0.0 

END  !  "of  discrete" 
END  !  "of  dynamic" 


END  !  "of  program' 


D.  ACSL  SIMULATION  FOR  CURRENT  FEEDBACK  CONTROLLER  IN  THE 
SYNCHRONOUS  REFERENCE  FRAME 

PROGRAM 


INITIAL 

" - This  Program  simulates  operation  of  a  3-Phase  hard  switched  inverter " 

"  with  sine-triangle  PWM.  Ideal  switching  and  components  are  assumed." 
"  Load  consists  of  an  LC  filter  with  a  resistive  load.  " 

"  A  Controller  operating  in  the  stationary  reference  frame  is  included  " 
"  DC  Voltage  =  500  Volts 


"  Switching  Frequency  =  6KHz 

MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.032 
CONSTANT  Lf=10.1e-3 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1 .0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*Vtripk*fitriang 
we=60.0*2.0*3. 14159265 
pi2o3=2.0*3. 14159265/3.0 
CONSTANT  twopi  =  3.14259265 
CONSTANT  iqssic=5.0 
CONSTANT  idssic=3.0 
iLa=iqssic 

ilb=-0.5*iqssic-idssic*sqrt(3.0)/2.0 

ilc=-0.5  *iqssic+idssic*sqrt(3 .0)/2.0 

Vag=0.0 

Vbg=0.0 

Vcg=0.0 

CONSTANT  K1=0. 10 
CONSTANT  K2=0.001 
CONSTANT  A=9.0 
Vmax=(Vdc/20.0) 

END  !  "of  initial" 


Output  Frequency  =  60Hz 

!  "maximum  integration  step  size 
!"min  time  step:for  var  step  algorithm " 
!"data  communication  interval 
! "integration  algorithm— R.K.  4th 

!"stop  point  for  integration  " 

!  "filter  (output)  inductance 

!"DC  source  voltage  " 

!  "load  resistance  " 

["reference  triangle  frequency  " 

!  "reference  triangle  peak  value  " 

["reference  triangle  period 
["reference  triangle  half  period  " 

["reference  triangle  slope  " 

["radian  frequency  of  contrlo  sinusiod" 
["120  degrees  in  radins 
["360  degrees  in  radians  ' ' 

["initial  q  axis  current 
["initial  d  axis  current  " 


"proportional  gain 
"integral  gain 
"desired  current  magnitude 
"pulse  width  modulation  gain 
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DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — sine-triangle  modulation" 

tt=mod(t,Ttri)  !  "form  ref  triangle  " 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0.0+slope*(tt-Tover2) 

END  IF 

" — calculate  synchronous  reference  frame  angle" 
thetae=integ(we,0.0) 

" — calculate  desired  currents" 

Idesa  =  A*cos(thetae)  ! "desired  phase  A  current  " 

Idesb  =  A*cos(thetae-pi2o3)  ["desired  phase  B  current  " 

Idesc  =  A*cos(thetae+pi2o3)  ["desired  phase  C  current  " 

" — transform  desired  currents  to  synchronous  qdO  quantities" 
Iqsdes=(2.0/3.0)*(Idesa*cos(thetae)+Idesb*cos(thetae-pi2o3)+& 
Idesc*cos(thetae+pi2o3)) 

Idsdes=(2.0/3.0)*(Idesa*sin(thetae)+Idesb*sin(thetae-pi2o3)+& 

Idesc*  sin(thetae+pi2o3)) 

"-—transform  actual  currents  to  synchronous  qdO  quantities" 

Iqse=(2.0/3.0)*(ILa*cos(thetae)+ILb*cos(thetae-pi2o3)+ILc*cos(thetae+pi2o3)) 

Idse=(2.0/3.0)*(ILa*sin(thetae)+ILb*sin(thetae-pi2o3)+ILc*sin(thetae+pi2o3)) 

" — determine  error  currents" 

Iqserr=Iqsdes-iqse 

Idserr=Idsdes-idse 

"-—generate  commanded  currents  (PI  Controller)" 

Iqscom=(Kl  *Iqserr  +  K2*integ(Iqserr,0.0)) 

Idscom=(Kl  *Idserr  +  K2*integ(Idserr,0.0)) 

"-—transform  commanded  quantaties  to  ABC  frame" 

Iacom=Iqscom*  cos(thetae)+  Idscom*  sin(thetae) 
Ibcom=Iqscom*cos(thetae-pi2o3)  +  Idscom*sin(thetae-pi2o3) 
Iccom=Iqscom*cos(thetae+pi2o3)  +  Idscom*sin(thetae+pi2o3) 
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" — normalize  voltages  for  use  in  switching  control" 
Vacont=Iacom*Rload/V  max 
Vbcont=Ibcom*  Rload/V  max 
V  ccont=Iccom*  Rload/V  max 


"——schedule  inverter  switching" 
SCHEDULE  hia.XN.(Vtri  -  Vacont) 
SCHEDULE  loa.XP.(Vtri  -  Vacont) 

SCHEDULE  hib.XN.(Vtri  -  Vbcont) 
SCHEDULE  lob.XP.(Vtri  -  Vbcont) 

SCHEDULE  hic.XN.(Vtri  -  Vccont) 
SCHEDULE  loc.XP.(Vtri  -  Vccont) 

" — calculate  line  voltages" 
Vas=(2.0/3.0)*Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 
Vbs=(2.0/3.0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 
Vcs=(2.0/3.0)*Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

" — transform  to  stationary  reference  frame" 

V  qss=(2 .0/3 . 0)  *  (Vas-0. 5  *  Vbs-0. 5  *  Vcs) 

V  dss=(sqrt(3 . 0)/3 . 0)*  (- Vbs+V  cs) 


["switch  A  to  HI  state  " 
["switch  A  to  LO  state" 

["switch  B  to  HI  state  " 
["switch  B  to  LO  state" 

["switch  C  to  HI  state  " 
["switch  C  to  LO  state" 


" — find  current  derivatives  in  stationary  frame" 
piqss=(V  qss-RLoad*iqss)/Lf 
pidss=(Vdss-Rload*  idss)/Lf 

" — calculate  stationary  currents" 

iqss=integ(piqss,iqssic) 

idss=integ(pidss,idssic) 

"——transform  back  to  ABC  quantaties" 
iLa=iqss 

iLb=-0. 5  *  iqss-(sqrt(3 .0)/2 .0)  *  idss 
iLc=-0 . 5  *  iqss+(sqrt(3 . 0)/2 . 0)  *  idss 

END  !  "of  derivative" 


DISCRETE  hia 
Vag=Vdc 

END  !  "of  discrete" 
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DISCRETE  loa 
Vag=0.0 

END  !  "of  discrete" 

DISCRETE  hib 
Vbg=Vdc 

END  !  "of  discrete" 

DISCRETE  lob 
Vbg=0.0 

END  !  "of  discrete" 

DISCRETE  hie 
Vcg=Vdc 

END  !  "of  discrete" 

DISCRETE  loc 
Vcg=0.0 

END  !  "of  discrete" 
END  !  "of  dynamic" 


END  !  "of  program' 


E.  MATLAB  SIMULATION  FOR  VLOTAGE  CONTROLLER  WITH  INNER 
CURRENT  CONTROL  LOOP 

%  Voltage  Controller  with  Inner  Current  Control  Loop 
%  State  Space  Model 

%  States 
% 

%  [xvqe  xvde  xiqe  xide  Vloadqe  Vloadde  iloadqe  iloadde  ilqe  ilde]' 

% 


%%  Gains  %% 

Kpv  =  0.25  %  Voltage  Proportional  Gain 

Kiv  =1.0  %  Voltage  Integral  Gain 

Kpc  =10.0  %  Current  Proportional  Gain 

Kic  =1.0  %  Current  Integral  Gain 

Kpwm=1.0  %  PWM  Gain 
U=50.0  %  Magnitude  of  Step  Change 

%%  Simulation  Parameters  %% 

we=  2*pi*60;  %  Electrical  Frequency 
Lf=.0101;  %  F ilter Inductance 

Cf=.0001;  %  Filter  Capacitance 

Rload  =  25;  %  Load  Resistance 

Lload  =  .0001;%  Load  Inductance 

%%  System  Matrices  %% 

A=[0, 0,  0,  0,  -1,  0, 0, 0,  0, 0; 

0, 0,  0,  0, 0,-l,0, 0,  0,  0; 

Kiv,  0,  0,  0,  -Kpv,  0, 1, 0,  -1, 0; 

0,  Kiv,  0,  0,  0,  -Kpv,  0,  1,  0,  -1; 

0, 0,  0,  0,  0,  -we,  -1/Cf,  0, 1/Cf,  0; 

0,  0,  0,  0,  we,  0,  0,  -1/Cf,  0, 1/Cf; 

0,  0,  0,  0, 1/Lload,  0,  -Rload/Lload,  -we,  0,  0; 

0,  0,  0,  0,  0, 1/Lload,  we,  -Rload/Lload,  0, 0; 

Kpc*Kiv*Kpwm/Lf,  we*Kic*Kpwm,  Kic*Kpwm/Lf,  0,  ((l-Kpc*Kpv)*Kpwm-l)/Lf, 
-Kpv*we*Kpwm,  Kpc*Kpwm/Lf,  we*Kpwm,  -Kpc*Kpwm/Lf,  -we; 

-we*Kiv*Kpwm;  Kpc*Kiv*Kpwm/Lf,  0,  Kic*Kpwm/Lf;  Kpv*we*Kpwm, 
((l-Kpc*Kpv)*Kpwm-l)/Lf,  -we*Kpwm,  Kpc*Kpwm/Lf,  we,  -Kpc*Kpwm/Lf;]; 
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B=[l,  0,  Kpv  -we*Cf,  0, 0, 0, 0,  Kpwm*(Kpc*Kpv/Lf-(we*we*Cf)), 
-Kpwm*  (Kpc*  we*  Cf/Lf+Kpv*  we)] ' ; 

C=U*[0  000  1  00  0  00]; 


D=[0]; 


%%  Determine  Poles  and  Zeros  %% 

[num,den]=ss2tf(A,B5C,D); 

R=roots(den) 

P=roots(num) 

%%  Display  Step  Response  %% 

figure(l) 

step(A,B,C,D,l) 

title(['Kpv-,num2str(Kpv),'  Kpc-,num2str(Kpc),'  Kiv=',num2str(Kiv), 
'  Kic-,num2str(Kic)]) 
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F.  ACSL  SIMULATION  FOR  VOLTAGE  CONTROLLER  WITH  INNER 
CURRENT  CONTROL  LOOP 

PROGRAM 


INITIAL 

" - This  Program  simulates  operation  of  a  3-Phase  hard  switched  inverter  " 

"  with  sine-triangle  PWM.  Ideal  switching  and  components  are  assumed." 
"  Load  consists  of  an  LC  filter  with  an  LR  load.  A  current  controlled 
"  voltage  regulated  pulse  width  modulation  (VRPWM)  control  scheme  " 

"  is  implimented  in  the  synchronous  reference  frame.  " 

"  DC  Voltage  500  Volts  VRPWM  Control 

"  Switching  Frequency  =  6kHz  Output  Frequency  =  60Hz  " 


MAXTERVAL  maxt  -  1 ,0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint  =  1.0e-7 
ALGORITHM  ialg  =  5 
NSTEPS  nstp  =  1 
CONSTANT  tstop  =  0.03 

CONSTANT  Lf=l  0.1  e-3 
CONSTANT  CF=1.0e-4 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  Lload=0.1e-3 
CONSTANT  ftriang  =  6000.0 
CONSTANT  Vtripk  =  10.0 
Ttri  =  1 ,0/ftriang 
Tover2  =  Ttri/2.0 
slope  =  4.0*Vtripk*ftriang 
we=60.0*2.0*3. 14159265 
pi2o3=2.0*3 .141 59265/3 .0 
twopi=2.0*3. 14159265 
CONSTANT  iqssic=0.0!5.5 
CONSTANT  idssic=0.0!6.0 
CONSTANT  iqssLic=0.0!6.0 
CONSTANT  idssLic=0.0 
CONSTANT  Kpv=  1.0 
CONSTANT  Kiv=  1.0 
CONSTANT  Kpc=  1.0 
CONSTANT  Kic=  1.0 
iLa=iqssic 

ilb=-0.5  *  iqssic-idssic*  sqrt(3 . 0)/2. 0 


["maximum  integration  step  size  " 

!"min  time  step:for  var  step  algorithm  " 
! "  data  communication  interval  " 

! "integration  algorithm-R.K.  4th 

! "  stop  point  for  integration  ' ' 

!"  filter  (output)  inductance  " 

["filter  (output)  capacitance 
!"DC  source  voltage  " 

["load  resistance  " 

["load  inductance  " 

["reference  triangle  frequency 
!  "reference  triangle  peak  value 
["reference  triangle  period 
!  "reference  triangle  half  period 
["reference  triangle  slope 
["radian  frequency  of  control  sinusiod" 
["120  degrees  in  radians  " 

["360  degrees  in  radians  " 

["initial  q  axis  converter  current  " 

["initial  d  axis  converter  current 
["initial  q  axis  load  current 
["initial  d  axis  load  current 
["voltage  loop  proportional  gain  " 

["voltage  loop  integral  gain 
!  "current  loop  proportional  gain 
["current  loop  integral  gain 


$ 
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ilc=-0. 5  *  iqssic+idssic*  sqrt(3 . 0)/2.0 

CONSTANT  VqssLic=150.0 

CONSTANT  VdssLic=0.0 

Vag=0.0 

Vbg=0.0 

Vcg=0.0 

CONSTANT  Vcom=150.0 
Vmax=(Vdc/(2.0*Vtripk)) 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — sine-triangle  modulation" 

tt=mod(t,Ttri)  !  "form  ref  triangle 

IF(tt.LT.Tover2)  THEN 
Vtri=10.0  -  tt*slope 
ELSE 

Vtri=- 1 0.0+slope*  (tt-Tover2) 

END  IF 

" — calculate  synchronous  reference  frame  angle" 
thetae=mod(integ(we,0.0),twopi) 

"---calculate  feedforward  capacitor  current" 
icqcom=(Vcom-VqssL)*Kpv  +  integ(Vcom-VqssL,0.0)*Kiv 
icdcom=-VdssL*Kpv  -  integ(VdssL,0.0)*Kiv  -  Vcom*we*Cf 

" — calculate  feedforward  current" 
ifqcom=icqcom  +  iqssL 
ifdcom=icdcom  +  idssL 

"——calculate  feedback  current  error" 
iqerr=ifqcom  -  iqss 
iderr=ifdcom  -  idss 

"-—calculate  voltage  control  signals" 

Vgq=iqerr*Kpc  +  integ(iqerr,0.0)*Kic  +  we*Lf*ifdcom  +  VqssL 
Vgd=iderr*Kpc  +  integ(iderr,0.0)*Kic  -  we*Lf^ifqcom  +  VdssL 

" — transform  control  quantaties  to  ABC  frame" 
Vacom=Vgq*cos(thetae)  +  Vgd*sin(thetae) 
Vbcom=Vgq*cos(thetae-pi2o3)  +  Vgd*sin(thetae-pi2o3) 
Vccom=Vgq*cos(thetae+pi2o3)  +  Vgd*sin(thetae+pi2o3) 
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" — normalize  voltages  for  use  in  switching  control" 

V  acont=V  acom/V  max 
Vbcont=Vbcom/V  max 

V  ccont=V  ccom/V  max 

" — schedule  inverter  switching" 

SCHEDULE  hia.XN.(Vtri  -  Vacont) 

SCHEDULE  loa.XP.(Vtri  -  Vacont) 

SCHEDULE  hib.XN.(Vtri  -  Vbcont) 

SCHEDULE  lob.XP.(Vtri  -  Vbcont) 

SCHEDULE  hic.XN.(Vtri  -  Vccont) 

SCHEDULE  loc.XP.(Vtri  -  Vccont) 

"——calculate  line  voltages" 

Vas=(2.0/3.0)*Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 

Vbs=(2.0/3.0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 

Vcs=(2.0/3.0)*Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

"——transform  to  synchronous  reference  frame" 

V  qss=(2 .  0/3 .0)*  (Vas*  cos(thetae)+Vbs*  cos(thetae-i2o3 )+ V  cs*  cos(thetae+pi2o3 )) 

V  dss=(2. 0/3 .0)*  ( Vas*  sin(thetae)+Vbs*  sin(thetae-pi2o3 )+ V  cs*  sin(thetae+pi2o3 )) 

" — find  converter  current  derivatives  in  synchronous  frame" 
piqss=-we*idss  +  (Vqss-VqssL)/Lf 
pidss=  we*iqss  +  (Vdss-VdssL)/Lf 

" — find  load  current  derivatives  in  synchronous  frame" 
piqssL=-we*idssL  +  (VqssL-Rload*iqssL)/Lload 
pidssL=  we*iqssL  +  (VdssL-Rload*idssL)/Lload 

" — find  load  voltage  derivatives  in  synchronous  frame" 
pVqssL=-we*VdssL  +  (iqss-iqssL)/Cf 
pVdssL=  we*VqssL  +  (idss-idssL)/Cf 

"—calculate  synchronous  converter  currents" 

iqss=integ(piqss,iqssic) 

idss=integ(pidss,idssic) 

" — calculate  synchronous  load  currents" 

iqssL=integ(piqssL,iqssLic) 

idssL=integ(pidssL,idssLic) 


!"  switch  A  to  HI  state  " 
!  "switch  A  to  LO  state" 

!  "switch  B  to  HI  state  " 
! "switch  B  to  LO  state" 

!  "switch  C  to  HI  state  " 
!" switch  C  to  LO  state" 
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" — calculate  synchronous  load  voltages" 

V  qssL=integ(p  V  qssL,  V  qssLic) 
VdssL=integ(pVdssL,VdssLic) 

" — transform  back  to  ABC  quantaties" 
iLa=iqss*cos(thetae)+  idss*  sin(thetae) 
iLb=iqss*cos(thetae-pi2o3)  +  idss*sin(thetae-pi2o3) 
iLc=iqss*cos(thetae+pi2o3)  +  idss*sin(thetae+pi2o3) 

VLa=V  qssL*cos(thetae)+  VdssL*  sin(thetae) 

VLb=V  qssL  *  cos(thetae-pi2o3 )  +  VdssL*sin(thetae-pi2o3) 
VLc=VqssL*cos(thetae+pi2o3)  +  VdssL*  sin(thetae+pi2o3) 

iLoada=iqssL*cos(thetae)+  idssL*sin(thetae) 
iLoadb=iqssL*cos(thetae-pi2o3)  +  idssL*sin(thetae-pi2o3) 
iLoadc=iqssL*cos(thetae+pi2o3)  +  idssL*sin(thetae+pi2o3) 

END  !  "of  derivative" 

DISCRETE  hia 
Vag=Vdc 

END  !  "of  discrete" 

DISCRETE  loa 
Vag=0.0 

END  !  "of  discrete" 

DISCRETE  hib 
Vbg=Vdc 

END  !  "of  discrete" 

DISCRETE  lob 
Vbg=0.0 

END  !  "of  discrete" 

DISCRETE  hie 
Vcg=Vdc 

END  !  "of  discrete" 

DISCRETE  loc 
Vcg=0.0 

END  !  "of  discrete" 


END  !  "of  dynamic' 
END  !  "of  program' 
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APPENDIX  E.  SPACE  VECTOR  CONTROL  SIMULATION  PROGRAM 


PROGRAM 

INITIAL 


This  Program  simulates  operation  of  a  3-Phase  hard  switched  inverter 
w/  Space  Vector  PWM.  Ideal  switching  and  components  are  assumed. 


MAXTERVAL  maxt  =  1.0e-8 
MINTERVAL  mint  =  1.0e-12 
CINTERVAL  cint=1.0e-7 
ALGORITHM  ialg  =  5 
CONSTANT  tstop  =  0.03 

CONSTANT  Lf=10.1e-3 
CONSTANT  CF=1.0e-4 
CONSTANT  Vdc=500.0 
CONSTANT  Rload=25.0 
CONSTANT  fswitch  =  6000.0 
LOGICAL  sequence 
sequence=.true. 
state=7 

Tsw  =  1.0/fswitch 
wosin=60.0*2.0*3. 14159265 
rl50=5.0*3. 14159265/6.0 
rl20=2.0*3. 14159265/3.0 
r90=3. 141 59265/2.0 
r60=3. 14159265/3.0 
r30=3. 14159265/6.0 
rm30=-3. 14159265/6.0 
rm90=-3. 141 59265/2.0 
rml50=-5.0*3. 141 59265/6.0 
CONSTANT  iqssic=8.0 
CONSTANT  idssic=3.0 
CONSTANT  AA=100.0 
iLa=iqssic 

iLb=-0.5*iqssic-idssic*sqrt(3.0)/2.0 
iLc=-0.5*iqssic+idssic*sqrt(3.0)/2.0 
CONSTANT  VqssLic=200.0 
CONSTANT  VdssLic=0.0 


"maximum  integration  step  size 
"min  time  step  for  var  step  algorithm 
"data  communication  interval 
"integration  algorithm~R.K.  4th  ' 

"stop  point  for  integration  1 

"filter  (output)  inductance  ’ 

"filter  (output)  capacitance 
"DC  source  voltage  1 

"load  resistance  ’ 

"reference  switch  frequency  ’ 

"true  =  7  to  8  ’ 

"initialization  ' 

"starting  switch  configuration  ' 

"reference  switching  period  ' 

"radian  frequency  of  desired  sinusiod 
"150  degrees  in  radians  ' 

"120  degrees  in  radians  ' 

"90  degrees  in  radians  ' 

"60  degrees  in  radians  ' 

"30  degrees  in  radians  ’ 

"-30  degrees  in  radians 
"-90  degrees  in  radians  ' 

"-150  degrees  in  radians  ' 

"initial  q  axis  current  ' 

"initial  d  axis  current  ' 

"commanded  voltage  magnitude 
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VLa=VqssL 

VLb=-0.5*VqssLic-(sqrt(3.0)/2.0)*VdssLic 

VLc=-0.5*VqssLic+(sqrt(3.0)/2.0)*VdssLic 

Vag=0.0 

Vbg=0.0 

Vcg=0.0 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

" — space  vector  modulation" 

Vcosa  =  AA*cos(wosin*t)  !"des  phase  A  Voltage" 

Vcosb  =  AA*cos(wosin*t-rl20)  !"des  phase  B  Voltage" 

Vcosc  =  AA*cos(wosin*t+rl20)  !"des  phase  C  Voltage" 

Vq=(2.0/3.0)*(Vcosa-0.5*Vcosb-0.5*Vcosc)  !"des  q-axis  voltage  " 

Vd=(sqrt(3.0)/3.0)*(-Vcosb+Vcosc)  !"des  d-axis  voltage  " 

"——determine  sextant  of  operation" 
beta=atan2(V  d,Vq) 

IF  ((beta.  GT.r90).  AND.  (beta.LT.r  1 5  0))  THEN 
sextant=l 

ELSE  IF  ((beta. GT.r  150). OR.(beta.LT.rm  1 5 0))  THEN 
sextant=2 

ELSE  IF  ((beta.GT.rml50).AND.(beta.LT.rm90))  THEN 
sextant=3 

ELSE  IF  ((beta.GT.rm90).AND.(beta.LT.rm30))  THEN 
sextant=4 

ELSE  IF  ((beta.GT.rm30).AND.(beta.LT.r30))  THEN 
sextant=5 
ELSE 
sextant=6 
END  IF 

" — calculate  switching  times" 

Vbar=sqrt(Vq**2+Vd**2)  ! "desired  magnitude  " 

IF  (Vbar.gt.Vdc)  THEN  !"limit  magnitude  " 

Vbar=Vdc 
END  IF 
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IF  (beta.GT.0.0)  THEN 

gamma=mod(atan2(Vd,Vq),r60)  !  "desired  angle 

ELSE 

gamma=mod(atan2(V  d,Vq),r60)+r60 
END  IF 

a=Vbar/Vdc 

K=Tsw*a/sin(r60) 

IF  (sequence)  THEN 
1 1  =bound(5e-6,T  sw,K*  sin(r60-gamma)) 
t2=bound(5e-6,T  sw,K*  sin(gamma)) 

ELSE 

1 1  =bound(5  e-6,T  sw,K*  sin(gamma)) 
t2=bound(5  e-6,Tsw,K*  sin(r60-gamma)) 

END  IF 

tzero=(T  sw-(t  1  +t2))/2 . 0 

"-—generate  switching  of  inverter" 
tt=mod(t,Tsw) 

SCHEDULE  switchtl.xp.tt-tzero 
SCHEDULE  switcht2.xp.tt-(tl+tzero) 

SCHEDULE  switcht3  .xp.tt-(t  1  +t2+tzero) 

" — calculate  line  voltages" 

V  as=(2 . 0/3 . 0)  *  Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 

Vbs=(2 . 0/3 . 0)*  Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 

V  cs=(2 .  0/3 . 0)  *  Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

" — transform  to  stationary  reference  frame" 

V  qss=(2.0/3 .0)*  (Vas-0. 5  *  Vbs-0. 5  *  Vcs) 

V  dss=(sqrt(3 . 0)/3 . 0)*  (- Vbs+V  cs) 

"-—find  current  derivatives  in  stationary  frame" 
piqss=(V  qss-VqssL)/Lf 
pidss=(Vdss-VdssL)/Lf 

"-—find  load  voltage  derivatives  in  stationary  frame" 
p  V  qssL=(iqss-(V  qssL/Rload))/Cf 
pV  dssL=(idss-(V  dssL/Rload))/Cf 
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" — calculate  stationary  currents" 

iqss=integ(piqss,iqssic) 

idss=integ(pidss,idssic) 

" — calculate  stationary  voltages" 

V  qssL=integ(pV  qssL,  V  qssLic) 

V  dssL=integ(pV  dssL,  V  dssLic) 

" — transform  back  to  ABC  quantaties" 
iLa=iqss 

iLb=-0. 5  *  iqss-(sqrt(3 . 0)/2.0)  *  idss 
iLc=-0.5*iqss+(sqrt(3.0)/2.0)*idss 

VLa=VqssL 

VLb--0.5*VqssL-(sqrt(3.0)/2.0)*VdssL 
VLc=-0 . 5  *  VqssL+(sqrt(3 .0)/2 . 0)*  VdssL 

iLoada=VLa/RLoad 

iLoadb=VLb/RLoad 

iLoadc=VLc/RLoad 

END  !  "of  derivative" 

DISCRETE  switchtl 
IF  (sequence)  THEN 

IF  ((sextant.EQ.  l).OR.(sextant.EQ.2))  THEN 
Vcg=0.0 
state=2 

ELSE  IF  ((sextant.EQ.3).OR.(sextant.EQ.4))  THEN 
Vag=0.0 
state=4 

ELSE 

Vbg=0.0 

state=6 

END  IF 
ELSE 

IF  ((sextant.EQ.  l).OR.(sextant.EQ.6))  THEN 
Vag=Vdc 
state=l 

ELSE  IF  ((sextant.EQ.2).OR.(sextant.EQ.3))  THEN 
Vbg=Vdc 
state=3 

ELSE 

Ycg=Vdc 

state=5 
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END  IF 
END  IF 

END  !  "of  discrete" 

DISCRETE  switcht2 
IF  (.notsequence)  THEN 
IF  ((sextant.EQ.  l).OR.(sextant.EQ.2))  THEN 
if  (state.EQ.l)  then 
Vbg=Vdc 
else 

Vag=Vdc 
end  if 
state=2 

ELSE  IF  ((sextant.EQ.3).OR.(sextant.EQ.4»  THEN 
if  (state.EQ.3)  then 
Vcg=Vdc 
else 

Vbg=Vdc 
end  if 
state=4 
ELSE 

if  (state.EQ.5)  then 
Vag=Vdc 
else 

Vcg=Vdc 
end  if 
state=6 
END  IF 
ELSE 

IF  ((sextant.EQ.  l).OR.(sextant.EQ.6))  THEN 
if  (state.EQ.2)  then 
Vbg=0.0 
else 

Vcg=0.0 
end  if 
state=l 

ELSE  IF  ((sextant.EQ.2).OR.(sextant.EQ.3»  THEN 
if  (state.EQ.2)  then 
Vag=0.0 
else 

Vcg=0.0 
end  if 
state=3 
ELSE 

if  (state.EQ.4)  then 
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Vbg=0.0 

else 

Vag=0.0 
end  if 
state=5 
END  IF 
END  IF 

END  !  "of  discrete” 

DISCRETE  switcht3 
IF  (state.EQ.l)  THEN 
Vag=0.0 
state=8 

sequence=.false. 

ELSE  IF  (state.EQ.2)  THEN 
Vcg=Vdc 
state=7 

sequence=.true. 

ELSE  IF  (state.EQ.3)  THEN 
Vbg=0.0 
state=8 

sequence=.false. 

ELSE  IF  (state.EQ.4)  THEN 
Vag=Vdc 
state=7 

sequence=.true. 

ELSE  IF  (state.EQ.5)  THEN 
Vcg=0.0 
state=8 

sequence=.false. 

ELSE 

Vbg=Vdc 

state=7 

sequence=.true. 

END  IF 

END  !  "of  discrete" 

END  !  "of  dynamic" 

END  !  "of  program" 
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APPENDIX  F.  SIMULINK  PROGRAMS 


A.  STATIONARY  REFERENCE  FRAME  CONTROLLER 


Desired  Currents  Ks  <3*aKis 

Controller 


Figure  F-l  -  Stationary  Controller 


Figure  F-2  -  Current  Sensing  Block 


Figure  F-3  -  Controller  Block 
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Figure  F-4  -  Pre-Processing  Block 


Figure  F-5  -  Desired  Currents  Block 


Ks23 


Figure  F-6  -  Ks  Block 


Ks-1(32) 


Figure  F-7  -  Ks  Inverse  Block 


Figure  F-8  -  Scale  to  +  1  Block 
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B.  SYNCHRONOUS  REFERENCE  FRAME  CONTROLLER 


The  Current  Sensing,  Pre-Processing,  Controller  and  Scale  to  +  1  blocks  for  the 
synchronous  controller  are  the  same  as  in  the  stationary  controller  above.  Blocks  Ke, 
Desired  Currents  and  Ke  Inverse  are  shown  here. 


Figure  F-9  -  Synchronous  Controller 


Figure  F-10  -  Ke  Block 
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Figure  F-ll  -  Desired  Currents  Block 


Figure  F-12  -  Ke  Inverse  Block 
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